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Because of the exchange interaction among the moment 
bearing magnetic atoms (or ions) in random magnetic alloys, 
various kinds of magnetic ordering take place depending on 
the tempo I'atuie and the concentration, c, of the magnetic 
co-nstituent(s) . In this thesis we have investigated the 
nature of the transition from paramagnet-to-spin glass state 
and the behaviour of various physical properties across the 
transition. 


Recently there have been several attempts to regard 
a spin glass as the prototype for many more complicated and 
less understood physical systems. Unfortunately, even the 
theoretical understanding of the physics of spin glasses 
is far from complete. There have been two distinct and 
apparently irreconcilable approaches - the cooperative 
phase transition model and the model of gradual blocking 
of superparamagnetic clusters. Both of these have had 
limited success. All the successes and limitations of 
these models have been critically reviewed in part A of 



this thesis 



AiJLJL; 


Part B highlights some further inadequacies of the 
cooperati-ve phase transition model of Edwards- Anderson type 
with single relaxation time. Part C proposes an alter- 
natiTre theoretical jframework whose predictions are, at least 
qu^itatively, closer to the experimental results. The 
latter pa-ves the way for an unified approach to the theoii.es 
of spin glasses by incorporating cooperative effects on the 
one hand, and features of a gradual freezing process on the 
other. 

Part B begins with Chapter 2 where we investigate 
the physical conditions under which both lorentzian and 
Gaussian local field distributions are good approximations 
for I sing spin passes with BKKY exchange interaction. In 
this chapter we also solve Sherrington-Kirkpatiick equations 
for two canonical spin glasses numerically and self-consis- 
tently so as to find the variation of the order parameters 
with temperature and external magnetic field. 

Mookerjee’s theory of magneto re si stance has been 
applied to two canonical spin glasses in Chapter 

Comparing the magnetoresistance thus calculated with 
the corresponding experimental data we calculate the strength 
of the s-d exchange interaction, 

The results of our investigation of relaxation spin 
dynamics in spin glasses using Master equation approach has 
been applied in Chapter 4 to calculate the resistivity by 
Boltzmann tiansport equation method* The resistivity thus 



calculated, using single relaxation time, shows cusp at T , 

o 

contrary to the smooth beha-viour observed experimentally. 

This is shown to be a consequence of the inadequacy of 
cooperative phase transition model with single relaxation 
time . 

Chapter 5 is an extension of Chapter 4 to vector spin 

glasses in the presence of external magnetic field. In this 

Chapter we improve Mookerjee’s theory of magneto resistance 

incorporating effects of spin dynamics. 

Part C of this thesis begins with Chapter 6 where we 

develop the percolation model of paramagne t-to-spin glass 

transition and calculate the transition temperature T as a 

a 

function of concentration c taking the effects of finite mean 

free path and fiustration properly into account. In this 

Chapter we derive Yogel- Pule her law and explain the variation 

of with time and length scales of probes. Besides, we 
a 

calculate the low- field low- frequency a.c. susceptibility and 
show that linear response theory holds. 

In Chapter 7 the mathematical relation between the 
percolation model and Wohlfarth’s model is established and 
interpreted physically. Comparing our theoretical results 
with experimental data we calculate the distribution of 
blocking temperatures. 

The logarithmic relaxation of magnetization is derived 
analytically in Chapter 8 so as to reveal a deeper connection 
with l/f noise in various other physical systems. 



ft are extended to calculate 
Ideas developed xn Chapter 8 

the resistivity of spin glasses m Ghaptt^r 9. 

. + r.-e random walk we calculate 

Generalizing the concept of rana 

^ - 1 . 4 -- v> rs-f ^i-oin clusters in Chapter 10 
the magnetization distribution o P 

for both classical as well as quantum sp' 

. . e^xistence of a 'hole* in 

In Chapter 11 we prove the exis 

^ ni ^ntlv low temperatures, 

the local field distribution at suffioi^a y 

model with some 

In Chapter 12 we compare our percoiatio 

^ ,,^e that "a spin glass 

other models heuristically and c 

is a glass”. 
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IITmODUGTIOU 


It is usual pmctice to begin a Ph.D. thesis with a 
brief review of the existing literatux^ on the subject of 
investigation followed by a set of emphatic statements that 
try to convince the reader that the investigation stemmed 
from a logical analysis, followed a rational course of 
reasoning and reached a firm conclusion. Besides, no work 
in theoretical physics can be relevant without comparison 
with experimental results. Moreover, a detailed critical 
review (instead of a bxief account) of the existing theore- 
tical models often helps in better assessment of their 
merits and at the same time throws light on their limita- 
tions. Therefore, I shall begin with a critical review of 
the theoretical models, a brief summaiy of the experimental 
results together with a comparison of the predictions of the 
former with the latter in part A. 

To achieve maximum logical coherence, details of 
our original contribution will be presented in two separate 
parts — B and 0| the order being far from chronological. 
Briefly speaking, part B of this thesis highlights some 
further inadequacies of the homogeneous mean field theoiy 
(with single relaxation time) of cooperative phase transition 
model of - Edwards- Anders on (EA) type for the paramagnet-to- 
spin glass transition. In part G we propose an alternative 
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thaoretical framework whose predictions are, at least 
qualitatively, closer to the experimental results. The 
latter also paves the way for an unified approach to the 
theories of spin glasses hy incorporating cooperative 
effects on the one hand and features of a gradual freezing 
process on the other. Some of the details of calculations 
have been given in Appendices 1 and 2. 
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PifflT A 

I®)AU FIELD THEORIES OF SPIU GLASSES; 
A CRITICAL REVIEW 



CHAPOIaR 1* 

mm MEIP THBOBIES OF SPIN GLASSES; A CBITIO-AL BBVIBW 


1 . 1 What is a Spin Glass ? 

The name "spin glass" (hereafter referred SG) was 
first introduced by Coles [ 1,2] to denote the "entire class 
of (random) magnetic alloys of moderate dilution in which the 
magnetic structure no longer resembles that of the pure metal 
but close enough so that their exchange interactions dominate 
other energies such as the Kondo effect and other free elect- 
ron interactions". Such an alloy, "after cooling to a very 
low temperature in zero applied field has the solute moments 
"frozen" in local molecular fields. These fields having a 
distribution of magnitudes and directions such that the net 
magnetization of any region containing a few teas of solute 
atoms is zero". So far as the nature of the exchange inter- 
action is concerned, according to Anderson [ 1 ], "the 
physical model that has come to be accepted for the substan- 
ces is that the exchange interaction follows the EKKI 
theory,. ..the very important property of EKII Yibloh gives 
SG the unique properties is tho oscillation". 


A review article based on this Chapter has been submitted 
for publication in Physics Reports. 
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Tills was the picture of SG- in 1972. During the last 
decade many o feher SG systems have bean discovered vyhsre the 
physical origin of the randomiiess and the nature of exchange 
int r-racnioiis are different from the ones considered above. 
Thsrefore, Binder [3] defined SG as "magnetic systems where 
the interactions between the spins are ’in conflict' v/ith 
each other, due to some 'quenched (i.e. frozsn-iii) disorder' 
in the system". The terms 'quenched disorder' and 'in 
conflict' will be discussed in detail later. So the essen- 
tial ingo^dient of SG is the random nature of the exchange 
interaction which does not allow the system to exhibit long 
range order (ferromagnetic or anti ferromagnetic) but imposes 
a new type of order where the spins are thought to be more or 
less 'frozen' in random directions. 

What is the justification for calling these systems 
'glass'? In glasses spatial degrees of freedom freeze 
without any long range order in space. Similarly, the frozen 
spin degrecss of freedom in SG exhibit no long range order, 
Notice that we have not conclusively stated anything about 
thd dynamical processes of this freezing. Later on vre shall 
see that some particular class of theories draw the analogy 
between different properties of glasses and SG and argue that 
SG aro also me tas table states like glasses. Various aspects 
of spin glasses have been reviewed in our recent article [ 4 ] 
and in many other theoretical [ 5-29] and experimental [ 30-43] 
review articles. 
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1 *2 Spin Q-Iass Systems : 

Mydosh [30], and later, Kumar [31] and Pord [35] 
enlisted and systematically classified the various types of 
SG-. We shall mention. ‘ some typical examples of each of 
these classes and then add some more recently discovered 
ones to the list. 

(a) Kohle metal- transition metal alloys; These are called 
canonical or archetypal SG-, typical examples being 
AuPe , QuTin, igIVin etc. The concentration of the transi- 
tion metal varies betv/een 0.1 to 10 at 

(b) Transition metal- transition metal alloys; The simplest 
member is PdJVIn [44] | others include MoMn, Rhivin, VPe etc. 

(c) Rare earth SO alloys; Some typical examples are 

(d) Transition metal compounds; An extensively studied 

member of this category is (^4. 0^ [45-50] . 

(e) Amorphous Metallic SO; This class [ 51] includes alloys 

of the type ^16®6 '^3 

and s’= (Mn, li, Or, Mo), Besides amorphous films 
of GdAl, laGd etc. also exhibit SG behaviour [ 52,53]. 

(f) Amorphous Semiconducting SG, e.g. amorphous Si, 
amorphous Ge [54] and (Sb 2 S^) 2 j.(SbI^)y.I'e 2 [55] • 

(g) Insulating and semiconducting SG; Simplest and widely 

studied systems are Bu^j. Sr^__^ S[ 5^ , Gd^ S[57] 

and Olg [58]. Besides 0d^_^ Mn^ Te and 
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Gd^_^ m^Se [ 59,6q O^^&MFeB'g 

[ 62 ] aad PbMaPeFr^ [63] also fall in this broad category, 
liany insulating amorphous aluniino silicate glasses 
[ 64-67 j also show SG- behaviour, 

(h) Intercalated graphite; Intercalating i^Cl^ into 
highly ordered graphite gives rise to SG properties 
of the system [ 68 ] . 

(i) 9 also an induced moment SG [69]. 

In this thesis we shall be mostly concerned with 
canonical SG. 


1 .3 Brief Review of Exnerimental Results ; 

The experiHBnts performed on SG can be broadly 
divided into two classes - one of which exhibits a sharp 
anomaly at a temperature T while the others show a smeared 

O 

behaviour (a peak or hump) not necessarily around T . 

s 


Sharp Anomaly 

Magne ti c Sus c e p ti bi li ty 

Remanence 

Mossbauer 


Smeared Behaviour 
Specific Heat 
Resistivity 
Thermopower 
Feutron Scattering 


Muon Spin Relaxation 
Anomalous Hall Effect 

We shall describe briefly experimental results of 
only the following properties - Mossbauer effect, magneti- 
zation and susceptibili-ty, specific heat, Resistivity 
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and magnetoresistance j and neutron scattering experiments. 

Our discussion will be mostly concerned with the experimental 
results on canonical spin glasses. 


1.5.1 Mo s s baue r E xpe ^ - 

The Mo'ssbauer experiment measures the local magnetic 
field at or near a particular site where a radioactive isotope 
emits or absorbs -rays. 

Hyper fine splitting, apparently below a sharp 
temperature, was observed by Violet and Borg [ 70] and by 
Window and coworkers [ 71,72]. 

Violet and Borg [70] divided the spectra into three 
classes, depending on their appearance (Figure 1 ) as follows: 

(I) spectra for which resonance lines are clearly resolved, 

(II) spectra where overlapping of the resonance lines 
blurrs the contrast, 

(III) spectra' which appear to be three-peaked . 

The areas under the absorption peaks, which is a 
measure of the line intensities established the random 
orientations of the frozen spins. Moreover, the total 
resonance absorption in the unsplit spectra equals that of 
the split spectra to within + 10^. This is remarkable 
because it implies that all of the Be atoms that exist in 
the paramagnetic state also exist in the magnetically 
ordered state thereby establishing the freezing to be a 
cooperative phenomenon, at least on a time scale of the order 
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Figure 1 : Resonance Absorption spectra of s*^iiiple AuFe (6,7 &t 

Concentration) at Various Temperatures tafter Eef . C7C| )• 
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of 10 seconds. Last, but not the least, important observa- 
tion was that there is no detectable unsplit paramagnetic 
component superposed on the magnetic hyperfine spectra. 

This implies that below T , the fraction of 'loose' spins 
haying a null hyperfine field is less than 

1.3.2 Ivlagne tization and Susceptibility Measuremerxts 

There aiv three aspects of magnetization Li(H, t,T) 
to be studied - 

i) variation of M with external magnetic field H and 
hysteresis , 

ii) variation of M with time t which reveals irieversible 
and nonequilibrium processes, i.e., magnetization 
relaxation, 

iii) variation of M with temperature. 

We know that magnetization and external magnetic 
field are related by response functions called magnetic 
susceptibility % 

M = H + + (1.1) 

where is the linear susceptibility (or zero field 
susceptibility) and X^, X 2 etc. are nonlinear susceptibi- 
lities. IFow onwards susceptibility will mean linear suscep- 
tibility, if not otherwise stated. But since zero field 
susceptibility is only an ideal concept not realizable in 
practice, one measu3?es very low field susceptibility instead. 
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Gamella and Mydosh [ 73] first mad;.-= successful 
measurements at low field (5Gr) and low frequency (155 Hz) 
and got th^ following interesting results: 

(a) X (T) exhibits a sharp cusp apparc.ntly at a sharp 
temperature f as shown in figure 2, 

O 

(b) The cusp is very sensitively field d«pendent (Piguire 3) 
and g‘,ts flattened even at low fields, eg.50G. They 
concluded that the failure of all previous measurements, 
who observed only broad maximum instead of sharp cusp, 
was due to the use of higher fields which was unavoid- 
able for limitations of accuracy, 

(c) The maximum of a.c. susceptibility always occured at 
a temperature lower than the temperature at which 
Mossbauer splitting occurs. 

later experiments [ 74,75] showed that T. shifts 

s 

with the frequency of the a.c. measurement (Figure. 4). At this 
point one should note that the frequency of the a.c, measure- 
ment can be varied to probe the dynamics of the freezing 
process over varying time scales. 

The significance of using low field is clear. But 
what is the use of the a.c. technique instead of d.c. measure- 
ment? The answer to this question is very closely related 
to the variation of magnetization with time, G-uy [ 76] 
performed accurate measuiement of d.c, magnetization in a 
field of less than twenty gauss. He observed that 



X no*"* emu 






'ipire 3 i Susceptibility Bata for ? /I-.-;'; with o ™ 5 and 8 at. 

showing the Curves for Zero ileld, and for Various 
Applied Fields (after Hef.[73j). 




A.C, Susceptibility of 
of T (after lief . [75] ) 
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(a) for systems cooled in zero external field magneti- 
zation attains a value Mq(T) witlcin a very short time 
after the external field is switched on. However, 
even if the field and temperature are kept constant 
the magnetization increases slowly with time so that 

M(t,T) = 1vIq(T) + Am(t,T) (1.2) 

On the other hapd, if the external field is turned off after 
a time t, ah^EI>j is left whose initial value is apparently 
A m(t,T) and which decays slowly in nearly the same way as 
Am(t,T) develops with the field applied (see figure. 5). 

Guy observed that 

Am(t,T) = s(l) In t (1.2a) 

He distinguished two quantities; 

(i) ideal nonequilibrium magnetization 

M(H,T,t=0) = Mq(H,T) 

(ii) equilibrium magnetization 

limM(H,T,t) = Ivi^(H,T) + lim Am(t,T) 
t -■ 00 t 00 

(b) equilibrium magnetization can be alternatively 
achieved by cooling in an external magnetic field, 

(c) in the a, c .measurement one measures only the rever- 
sible part of the magnetization, and hence the 
difference between a,c. and d.c. susceptibilily. 



ir- 




figure 5 s Tte TTirh-'!“d9poiident Beha^-lour of Magaeti nation 

at 4»'c k'. 2he Lower Curve shows the Levelopnseat 
of Magnetiaatioa with Tiine la an Applied field 
of 200- * The Upper Curve shows the deoay* of 
the IRSl obtained after switching off the field 
of 20C at t = 8 Kin. (after Ref. M )♦ 
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The l::^ste27esis curve differs from that of a ferromagnet 
because the remanent magnetization for a SG is time-depen- 
dent and decays to the equilibrium value (M(t = oo ) = 0 
for H=0)very slowly for all T < T^. The other interesting 
feature of the hysteresis curve is that it is laterally 
displaced from the origin as shown in Figure 6. 

1.3.5 Specific Heat Measurements ; 

First of all, we should recognise that it is extremely 
difficult to separate out the electronic contribution to 
specific heat from that of magnetic part because both vary 
as proportional to temperature T, although separation of the 

3 

phonon part is not very difficult as it varies as T , It is 
generally assumed that magnetic contribution to specific heat 
= (specific heat of the system containing magnetic impurity)- 
( specific heat of the system, under exactly same conditions, 
where magnetic impurities have been substituted by nonmagnetic 
impurities of atomic number very close to that of the magnetic 
impurity). It was observed [ 77] 

(i) that magnetic contribution to specific heat 
exhibits only a broad hump whose maximum occurs at 
a temperature higher than that of susceptibility, 

(Figure 7), 

(ii) that a plot of 0^/T also reveals no anomaly for most 


of the systems. 
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figure 7 s Specific heat of AUq in the temperature 

region 3-50K, The solid curve is the calculated 
iwnmametic contribution to the specif io heat of 
the alloy betvieen 0 and 30K. The insert showa 
toe eusce|tibility results of the sample 
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(iii) that only 30% of the total entropy is recovered upto 
T , the remaining 70% develops above T , 

O g 

(iv) that the specific heat is linear in I at very low 
temperatures. 

! 

i 

; Yexj recently Thomson and Thompson [ 78,79] have 
presented evidence for behaviour of low temperature speci- 

fic heat of SG, exactly as in case of resistivity, to be 
discussed later,. This type of T^'^^ behaviour was already 
suspected by Goles [ 2 ] in 1972. 

1.3*4 Resistivity and Magnetoresistance ; 

Ford and BAydosh [ 80] measured the magnetic or spin 
resistivity P (T) at temperature T, 

P (T) = ^ alloy "* ^pure host 

as shown in Figure 8.How8ver, such a method becomes inaccurate 
if deviation from Mattl^i^aen’ s rule becomes appreciably large. 
Their observations can be summarized as follows i 

(a) residual resistivity cc c for lower concentration 

alloys, 

(b) for c > 5 at % variation of P^ with e becomes slower 

than linear and at about 12 at % reaches a maximum. 

On further increase of c, gradually levels off at 

c > 25 at % Fe . 
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Hgure 8 t Over-all temperature variation of f cm) 

for CuMn alloys with concehtrationa between 0#7““ 
and 9.7- at ^ Mn. Note the changes in soale 
as the concentration increases (after ref. C801)* 
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(c) at low t^^mperature ^(T) risas steeply, initially 
faster than T, and then praotically linear with T in 
the region of T . Yfith further increase of T, P 

o 

reaches a very broad maximum at a temperature T 

mcAji. ' g 

and then gradually falls. 

(d) upto the lowest tomperatuie of 0.5K best fit with 
experimental data was obtained with 

P([r) = Pq ^ 

(e) the range of the fit increases with increasing 

concentration 

(f) the coefficient A varies with c as shown in Figure 9 

(g) T the tempei'ature corresponding to the maximum in 

P , increases less rapidly than c, roughly between 
and o^*^. The qualitative behaviour of T versus 
c is similar to that of T versus c. 

O 

Senou-ssi [ 8l] observed that for zero field-cooled. 

AuFe SG- the transverse and longitudinal magne tore si stance are 
the same. He showed that Transverse Magnetoresistance (TMR)/p 

^ ^ (1.4) 

where a is a constant. Senoussi used this expression and 
Higam and Majumdaxfs [82] experimental values of TIE to 
calculate X . These calculated values of X agreed quite 
well with those measured directly. 

The most detailed study of MR i-n canonical 5G has 
been made by Higam and Majumdar [ 83j84] . They observed 




Mgur. 9 •. Oo.melent of the T^/* depondanoo A (ofieivt’/®) 
against log^^ c for ^Or, and AuFe all 

(after ref.CSOj). 




24 


that — p — a.'^ where n >1. At low fields n £f'2 for all 

the alloys studied. The field H upto which the quadratic 

Vi 

nature of AP /P persists is different for different alloys. 

AP 

Besides at low fields, — ^ — is weakly dependent on temperature 

for all T < T . 

g 


1.3*5 Neutron Scattering ; 

The q-dependent susceptibility X(q) was determined 
by Murani aod coworkers [ 37”'39j85,86] from neutron scattering 
experinents using the relation 


d g 
d9 


a- k^T 


X(q) 


(1.5) 


Some remarkable features emerged - 

(a) the temperature T^ where the peak in X (q) occurs was 
always much higher than the temperature corresponding 

to the cusp in a.c. susceptibility, as shown in Figure. 10-. 

(b) the freezing temperature Tj, is q-dependent- higher the 
q- value lower is T^, as shown in Figure 11. 

(c) there is a distribution of relaxation times. 


1.4 Theories of Spin G-lasses ; 

Mydosh [87] commented that it is easier to discover 
a new S& material than to understand the behaviour of those 
already known. The nature of the SG- transition in even the 
simplest systems e.g. AuFe and GuBfoi is contro-versial. In 
this field of research, as in many other branches of condensed 
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figure 10 


TKe auaceptibtlity %(q) versus f for Cu-8 
at *5^ Mo. alloy. The dashed curve Is the 
>neasured a,c. magnetic susceptihilily of 
tWe saioe large sample as used for the 
neatrou scattering experiment (after ref.[37J) 
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i’xi^ure 11 j Snail aaglo -neutron scattering intensity aa a 
I'i^^iotion of thu temperature for aseries of 
q-values for Au^ alloys containing (a) 10 at 
% F& and (b) TJ at % Fe, I’be arrows mark tbs' 
v-;.'!pora tares of the discontinuity in the scatteilfig 
■*; oh are abov.n function of q in the insets 
(/.-.i'tcr ;ef*i37j, 
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matter pt^sics, theory is lagging far behind experimental 
developments. So far almost all conventional methods have 
had limited success. The solution may require the introdu- 
ction of some totally new concepts. Many new ideas have 
already been introduced (e.g. , frustration) and many more 
have been borrowed from other fields, e.g. the ideas of replica 
and long range order in time borrowed from polymer physics, 
concepts of gauge invarianoe and spontaneous symmetry break- 
ing from quantum field theory. Yet the whole picture, to 
this date, remains unclear. Moreover, a theoretician must 
also decide whether all kinds of SG can be explained by some 
general universal theory or different classes of materials 
showing the SG transition need different theoretical treat- 
ments. At present we can only say that the situation is not 
ripe for reaching any final conclusion. 

In any microscopic theory in physics one has to know 
the nature of the interaction. Before developing a theory of 
SG one has to construct a reliable model for the form of the 
impurity-impurity interaction. Here one once again comes 
across an age old problem in the theory of magnetism. The 
Heisenberg Hamiltonian 




1 s 

2 ij 




( 1 . 6 ) 


is a good description of the system provided there exist 
’good 'local moments at the iiipurity sites. Since we have 
already assumed the existence of such moments in article 1.1 



28 


while defining a SG, one can woik with either the Heisenberg 
model or the simpler models e.g. Xf model or I sing model. 

But the definition of SG in article 1.1 is only a clewer 
trick to bypass the question - "can there exist situations 
where the SG has to be described in teins of itinerant 
model [88]?" Ooles [ 89 ] has shown that there exist such 
situations e.g, AuCo, RhHe, RhGo, etc. These are called 
Stoner glass to distinguish them from the ones defined in 
article 1.1 . The defining criterion for a Stoner glass is 
that the SG transition temperature T^ is much less than the 
Kondo temperature T^^, i.e, T^ « Tj^. Hertz [ 90 ] has develop- 
ed the theory of Stoner glasses. But in this thesis we shall 
not discuss them further, 

Next, one has to know the nature of the exchange 
interactions ■ between the ith and jth spins. The nature 
of this exchange interaction, i.e, whether direct or indirect, 
EKKY or superexchange, has to be decided by experiment. Of 
course, the essential feature of this exchange is the compe- 
tition between ferromagnetic and antiferromagnetic inter- 
actions, as defined in article 1,1. The dominating inter- 
action in cannonical SG with which we are primarily interes- 
ted is RKKY [ 91 - 94 ] in nature. But its explicit form has 
been questioned from time to time. It was De Gennes [95] 
who first showed that in a disordered system the expression 
for the effective indirect exchange is 
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, Jsd (Vo/2) r 2k^R,,Cos(2k^,,)-Si a(2V-,j )^ 3 j/4 

(1.7) 


whe ro 5 


h k..T 


is the mean free path of tho conduction 


electrons and arises due to the effect of disorder j 


volume of unit cell, and 

eV| 


later similar explicit analytical expressions were derived 
by Kaneyoshi [96] and Poon [97] independently, which were 
shown to be equivalent mathematically by Kaneyoshi [ 98] . 

Very o^scently de Ghatel [99] has re-examined the validity of 
De Gonnes' expression and concluded that the latter is a good 
approximation when the disorder is weak thereby confirming 
Kaneyoshi ‘s conjecture based on model calculations that as 
mean free path becomes comparable with the average coherence 
length of disorder, Permi energy approaches the mobility edge 
implying more and moie deviation from le Gennes" approxima- 
tion (called "factorization approximation" by de Ghatel). 

Since the exchange interaction in SG is random 
due to the quenched disorder of the system one has to assume 
or derive the distribution of the relevant random variable 
(either the random exchange inte3?aotion or the random local 
field) ?^nd has to invoke the idea of spatial ergodicity 
[100], common feature of all theories of disordered systems. 
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So far as the configuration averaging is concerned, 
random magnetic systems are more complicated than the random 
non-magnetic systems, because of the added degrees of freedom- 
spin. This causes a difficulty that can he avoided by 
introducing some mathematical tricks which, in turn, introduce 
many more problems whose solutions are still far from satis- 
factory, In random magnetic systems one has to carry out two 
averages - the thermal and configuration averages. So, the 
next obvious question is which should be carried first? In 
other words, what physical quantity should be configuration- 
averaged-partition function, free energy, magnetization or 
specific heat etc? It was Brout [101] who, in his work on 
random ferromagnetism, answered this question clearly for 
the first time . In annealed systems the "impurity" degrees 
of freedom { 1} are in thermal equilibrium with other degrees 
of freedom, e.g. spin i S}, the variables iS] and ^1} can be 
treated on the same footing and therefore the partition 
function Z has to be averaged over all configurations. On 
the other hand, in a quenched system at a low temperature T 
the spin degrees of freedom {S} are in thermal equilibrium 
whereas the impurity configuration {1} is the equilibrium 
configuration corresponding to a very high temperature. 
Therefore {1} cannot be treated on equal footing Y/ith { S} . 

In quenched systems one first calculates the free energy 
corresponding to a given configuration and then averages it 
over all impurity configurations. 



1.4.1 Theories of First Kind 


These are the oldest theories of SG- alloys and are 
based on the idea of a distribution of random local magrietic 
fields. The first among these that explained the low temper- 
ature linear specific heat was that proposed by Marshall [ 102] . 
He assumed (a) that the spins interact via EKKY interaction 
and (b) that the spins are Ising spins (S = l/2). He estimated 
the probability P(H,T), the probability of local field H at 
temperature T,in an intuitiTe manner as discussed below: the 
local field distribution for feiromagne tie or antiferromagnetic 
systems are S -function in ?/eiss molecular field approximation 
but broaden out in the improved Bethe-Peierls-Weiss (BPVV) 
approximation. In case of dilute alloys one would expect it 
to broaden and flatten out more so that the width of the peaks 
become comparable to mean magnitudes of the local field itself, 
Hollowing Kittel's suggestion that this problem is closely 
analogous to the computation of "the line shape for Hlffi on 
dilute sample of nuclear spins distributed randomly in a 
crystal and coupled only by their magnetic dipole interactions", 
Marshall argued that the local field distribution, in analogy 
with HlviB' line shape, should also be a cutoff Lorentzian. Of 
course, Marshall stressed that the exact shape of the distri- 
bution is not important for explaining the temperature and 
concentration dependence of low temperature specific heat 
provided P(H,T) has a width proportional to the concentration 
c for dilute alloys. 
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The magnetic energy is given by 

h = - -I He r dH P(H,T) -iH tanh( £- ' ^ -- - -) (1.8) 


where H = number of sites, and the factor l/2 ensure no 
double counting. Hence, the specific heat is given by 


00 


0^(T) 


2 ..2 


He / dH P(H,T)(i^Y-S^) Sech^ (^) 

^ , T jB 

-00 B 


00 


- j Nc(-^) h(T) + -I He 


J 


-00 


dH {pH - pH tanh } 

B 


3P(h.t) 

3 T 


( 1 . 9 ) 


where 


h(T) 


00 

r dH IH 1 P(H,T) 
-”oo 


( 1 . 10 ) 


is the mean magnitude of H, 

At very low temperatures (T'<^^^-^) most of the spins 
are rigidly pinned by their local fields and with further 
decrease of temperature there can be no further appi^ciable 
increase in the number of such frozen spins. Since we also 
know that only such rigidly alignsd spins make appreciable 
contribution to h, we conclude that the temperature vairiation 
of h is negligibly small. Therefore, the second term in 
equation (1.9) can be dropped. Moreover, in the first and 
third terms only a small range of H(from — p; — to -p—) makes 
dominant contribution, because it is only in this range that 
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thermal agitation dominates over magnetic order* Therefore 
I-'(H) can he replaced by P(0). In addition, for any given 
reasonable variation of P(0) with T, the third term is 
negligible compared to the first. Hence, we finally get, 
after a change of variable 

00 

k T 

G^(T) = He P(0,0) J dx x^ Sfech^x (1.11) 

-00 

and hence the linear behaviour in T. 

later Marshall' s intuitive theory was put on a 
rigorous basis by Klein and Brout [103] (KB). Their main 
contribution was the derivation of the local field distri- 
bution from very basic considerations. They confirmed 
Marshall's intuitive suggestion that the local field distri-' 
bution should be lorentzian. Unfortunately, Lorentzian is 
not a veoy well behaved distribution in the sense that all 
its higher moments diverge, thereby inplying, in case of SG, 
infinite energy per spin. KB realised this drawback. They 
showed that this difficulty can be removed by paying attention 
to the real physical situation where there exists a lower 
cutoff in the relative distance between two impurities 
because the lattice structure is discrete. Y/hen the latter 
fact was properly taken into account the local field distri- 
bution turned out to be Gaussian which still has its width 
proportional to the concentration c, which was the only 
essential feature that Marshall assumed. 
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The maximuia observed in the susceptibility was 

explained by KB as follows : first they showed that there 

exists a characteristic correlation length beyond which 

< S S . > = 0 such that the "impurity spins are fully or 

^ 3 

partially correlated to the spins at the origin if they are 
located within radius E and are approximately randomly 

f? 

oriented if they are located outside of R^. Thus the solid 
is made up of small clusters of impurities. As the temper- 
ature is raised starting from very low temperature two 
competing phenomena take place - the clusters get randomized 
resulting in a decrease in the susceptibility; and clusters 
breaking up thereby increasing the susceptibility. At low 
temperatures the latter dominates whereas at high temperature 
the former dominates. As a result there occurs a maximum in 
the susceptibility. We shall not discuss the quantitative 
details of the above qualitative idea, firstly, because 
zero field susceptibility exhibits a sharp cusp rather than a 
maximum, and secondly, because the physical picture used by 
KB is contradictory to reality as will be shown very soon. 

Klein [ 104] discussed quantitatively the nature of 
the local field distribution and its variation with concen- 
tration and temperature. Since the inc3rease of concentration 
results in a decrease of E^, the effect on P(H) is "to 
remove some probability from the central portion of P(H) 
and to place it into the wings". Unfortunately, the 
variation of P(H) with temperature turned out to be unphysical— 
the most probable field shifts to higher values as the 
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temperature is raised. This is so because of the following 
reason: a spin situated in between two others can screen off 
the interaction. Higher temperature implies weaker screening 
and hence higher effectiye field. This is contrary to 
Marshall's qualitative conclusion: "the effect of correlation 
(at lower temperature) must be to reduce the energy and this 
can only be done by arranging for spins to sit in fields of 
larger magnitude". 

The possible source of error in Klein's theory was 
pointed out by Liu [105] . He stressed that the screening 
effect, if present at all, will be complete only at O^K. 

At any finite temperature screening may be ne^igibly small 
and hence KB approximation breaks down. Klein extended the 
formulation to include the effect of external magnetic field 
[ 10^ and to explain the high temperature properties [ 107] • 

One interesting feature of the internal field distri- 
bution P(H, T) is that it obeys the scaling law [ 108] 

P(H,T) = f: (|| T/c) (1.12) 

for c 0. 

Some other scaling laws follow from this form, e.g. 

= K(T/c) 

= r(Vo) (1.13) 

Such scaling laws are consequences of the power law behaviour 
of the form l/i^ of the interaction. One may wonder what is 


X (T) 
AG 
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the use of these scaling laws. The beauty of these laws 
lies in the fact that K(x) and n(x) being concentration 
independent functions, properties for all concentrations can 
be derived from those of a single concentration for any 
given system, e.g. AuFe . BKKY interaction has an additional 
unique feature that the exponent n=3=space dimensionality 
of the system. This featU3:« leads to another simplification 
and universal behaviour - the properties of all the canonical 
SG systems for all concentrations can be derived from those 
of a single system for a single concentration. An excellent 
discussion on these scaling laws can be found in Rammal and 
Soule tie [15]. 

It appears that the success of the theories of 
first kind was the explanation of low- temperature linear 
specific heat. But this explanation is also not veiy satis- 
factory [109] • One does not get linear specific heat for 
Heisenberg spins in the same theoretical framework, because 

P(H) oA 

However, Eivier [110] counterargued that "the relevant 
quantity, hov/ever, is the finite P(H^) rather than P(H), 
because the local magnetic energy S. H involves a scalar 
product, hence a preferred direction (that of S)", i.e., 

E = j* dH P(H) S . H 


(1.14) 
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where is the Biillouin f auction for the spin S in a 

field Anderson [ 111] rejected this argument saying 

"confusion on this point is widespread. .. .it is quite clear 
that tL.H = pH". The latter situation arises because, in 
equilibrium, is parallel to H. We shall see an alternative 
explanation of linear specific heat in article 1,7.2.1yiore 
about the local field distribution will be discussed later. 

The drawbacks of this kind of theories are summarized 

below 

(a) they predict a rounded maximum in zero field suscepti- 
bility X instead of a sharp cusp, 

(b) the magnitude of X(T), contrary to experinent, is 
independent of concentration, to first order (i.e., 
for low concentration) , 

(c) that the peaks in X(T) are independent of fields (for 
low fields), contrary to experiment, 

(d) that the temperature of the maximum should be directly 
proportional to concentration whereas experiment shows 
weaker- than- linear dependence. 

These points had been cited by Cannella [112] as failures of 
these earliest theories. 
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1.4.2 Theories of Second KiM 

Though the theories included in this categoiy may 
differ from each other in their quantitative details, they 
share a common feature, viz., they show SG transition as a 
genuine magnetic phase transition. There are five steps in 
the study of a phase transition: 

(i) identification of the relevant physical character- 
istic (s) that distinguishes the two phases of the 
^stem, 

(ii) quantification of the distinguishing character- 
istic (s) found in step (i) hy defining order 
parameter(s), 

(iii) developing a mean field theory (MFT) in terms of 

the free energy as a function of the order parameter, 

(iv) -verification that in the long-range limit (i.e., 
infinite range model) MPA gives exact answer, 

(v) improvement of MPT results by incorporating critical 
fluctuations in the theory e.g. by the method of 
renormalization group. 

Unfortunately, in the case of SG none of these 
steps ha-ve yet been completed very satisfactorily. However, 
we shall follow this scheme in our discussion of the SG 
'phase transition' . 
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Steps I and II ; 

Identification of the relevant order parameter is 
very crucial, because a system can exhibit several phase 
transitions characterized by different order parameters. We 
know that spontaneous magnetization is a measure of long-range 
ferromagnetic order whereas staggered magnetization is a 
measure of long range antiferromagnetic order. Hence, spon- 
taneous magnetization and staggered magnetization distinguish 
ferromagnetic and antiferromagnetic states from paramagnetic 
one. In other words, spontaneous magnetization and 
staggered magnetization play the role of order parameters for 
the corresponding phase transitions. After Cannella and 
Mydosh's classic observation of a sharp cusp in the low-field 
a.c. susceptibility physicists tried to explain it as due to 
the onset of antiferromagnetic order. But it could not be 
fitted within the existing ideas of long-range antiferromag- 
netic order because 

(a) X and ^ below the cusp temperatuia were different 
for aatiferro magnetism but identical for SG-, 

(b) the cusp was very sensitively field dependent, unlike 
antiferromagne tism, 

(c) neutron scattering [ 113 ] showed conclusively absence of 
any long-range order. 

Therefore it was obvious that SG is distinct from 
ferromagnetic or antiferromagnetic phases and requires new 
order parameter. Edwards and Anderson [114] introduced, for 
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SG, an order parameter Q, now popularly known as Edwards- 
Anderson (BA) order parameter, defined as 

Q= [<Sj_> (1.15) 

where < > denotes thermal- average and L ] denotes confi- 
guration average. In order to describe the dynamics of the 
system they introduced Ml 5] a dynamic order parameter 

Q(t) = ' S. (0)>] (1.16a) 

av 

Notice that Q.(t) is always positive - it relaxes to its 
static value for t oo only via values all greater than 
zero. The physical significance of EA order parameter is very 
important. In Anderson's own words [6] "we decided that since 
we believed, there was Xio IRO in space, we had to introduce 
a LRO in time". So, the dynamic order parameter €l(t) describes 
on-site autocorrelation in time, i.e., the probability that 
the spin at a particular site will be 'oriented' in a parti- 
cular direction, at time tjif it was 'oriented' in the same 
direction at time t=0. This is, of course, a classical 
interpretation. In quantum mechanical language, '’(t) 
describes the probability that the spin at a particular site 
will be an eigenstate if it was in the same eigenstate at 
t=0. 

Unfortunately, the choice of the order parameter, 
for SG phase, is not universally agreed upon. Recently, 
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Khurana and Hertz [116], from a study of order parameter 
susceptibility concluded the insufficiency of the single 
order parameter Q to describe the SG phase. They also 
indicated the possibility tf introducing generalized order 
parameters involving higher order cumulants. Sompolinsky 
r lltl introduced two order parameters Q and A , to be 
discussed later. Plefka [lis] utilized the generalized 
order parameters 

= I S (1.17) 

i 

Pari si [ 119,120] , in yet another context, remarked that "it 
is quite likely that an infinite number of order parameters 
is needed in the correct treatment". Parisi [ 121] introduced 
a generalized order parameter Q(x) where the variable x is 
defined over the interval [0, 1] . 

Suzuki and Miyashita [122] have discussed in detail 
different kinds of definitions of order parameters for SG. 
They, however, defined an order parameter and many of the 
earlier order parameters turned out to be special cases of 
this order parameter. 

Binder and Stauffer [12^] proposed to analyze the 
nature of the ordering in SG and relaxation from one 
degenerate ground state to another in terms of some new 
order parameter, let the set of phase factors 

, defined by the direction of spins 
C S^ , . . . .Sjj-} in the -1-th ground state, define an H-dimen- 
aional snace with unit vectors 
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5 


( 1 ) 





( 1 ) 




(1.18) 


They thought, "it is interesting to study the "similarity" 
of two ground states" produced by repeated aging and 
quenching in a Monte-Qarlo computer experiment. They defined 
the order parameter 




(1.19) 


which is just the projection of the state 


X = 


1 


Sfj} 




onto the 1-th ground state defined by the unit vector 

^ ^,^(1),..., 41)! 


The question whether the 1-th ground state has a non-zero 
order parameter with respect to m-th ground state is 
described by whether or not 

(1) (m) 

Cos = S .6 

( 1 ) 

vanishes. By starting from some direction 5 one can 
"rotate" in a higher dimensional space which physically 
implies relaxation of the system from one ground state to 
another. Fotice that looks at a particular config- 

uration of all the spins whereas the order parameters defined 
earlier look at a single spin averaged over different 
configurations. 



43 


Kirkpatrick aad Young [ 124,125] discussed the 

drawbacks of the usual definition of EA order parameter for 

their computation on a finite lattice. An external field, 

however small, was required to break tho symmetry in order 

to generate a non-zero <S^> in that particular direction 

and, then, after taking the thermodynamic limit, one can put 

the external field =0, They estimated that < S^> is 

* T * / 

nonzero for h >> h = - — and h-nce h 0 ±n finite 
lattice. One should, however, remember that the field applied 
hj^ = e^h with = + 1 with equal probability so that the 
process of symmetry breaking differs from that in a ferro- 
magnet where a uniform field is used. The utility of the 
symmetry breaking lies in the fact that only states in the 
vicinity of one of the ground states are included in the 
ensemble average. In case of large but finite systems 
Morgenstern and Binder's [126] order parameter 

lim (1.20) 

X av 

with short range interactions is very convenient one since 

S.S. does not change sign when all spins are reversed. Young 
i 

and Kirkpatrick generalized this order parameter to infinite- 
range model 

(i^j) (1.21) 

av 

The superscript (2) suggests that is something like ■tjp 

but not exactly equal to this". 
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Adkins and Ei-vier [127] assumed that the local 
magnetization, a short range order parameter, describes the 
SG- phase. This differs from BA order parameter hj the fact 
that the latter does not assume any SRO in space but only a 
LEO in time, later MedTedev [128] generalized EA order 
parameter to include short range magnetic order by writing 


where even for i i , and 0. . -* 0 as IR. - R. 1 -♦ oo . 

The last condition is necessary to avoid LEO. Medvedev's 
order parameter reduces to EA order parameter if 



Sompolinslsy [ 129] suggested the introduction of staggered 
magnetization associated with the eigenvalue spectrum of the 
random exchange matrix as the SG order parameter. The stagg- 
ered magnetization is defined by 


= i: < i 1 A > (1.23) 

<il?v> being the orthonormal eigenvector of the random matrix 
Jfj corresponding to the eigenvalue 


Step III ; 

Edwards and Anderson [ II 4 ] hoped that their theory 
will play the same role as Weiss molecular field theory did 
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for LRO phases. In order to average the configuration 
dependent free energy over all possible configurations BA 
utilized a technique called "replica trick". The trick 
converts the quenched problem into an annealed problem by 
using the identity 


In Z = lim ^ (z^ - 1 ) 

n 0 


(1.24) 


So, effectively what one averages is Z^. Por integral n, 

=lt 

U=1 

can be interpreted as the partition function of a composite 
system consisting of n non-interacting replica of the system 
and hence the name. More details of the replica method will 
be discussed in the next section. BA had to modify the 
definitions of their order parameters in the replica approach. 
BA assumed 

\ =[<S?> ] (1.25) 

av 

and 

Qctp = [<SJ>< S| >] (a^p) (1.26) 

av 


as the only relevant order parameters, and neglected all other 
combinations including the higher order ones. At this point 
one should notice that Khurana and Hertz's [116] conclusion 
about the insufficiency of the above order parameters, as 



stated in the previous section, can be anticipated here. One 
has to have infinite number of order parameters if the 
successive higher order terms including three or more spin 
operators are not negligible. Moreover, there does not exist 
any mathenatically rigorous proof that the order parameters 
defined by equations (1.25) 25ud(1.26) reduce to M and Q in the 
limit n - 0 . 

EA assumed adhoc that the spins are classical spins 
and the exchange interactions obey Gaussian distribution. 

They were able to show the existence of a critical tempeD>- 
ature T above which q vanishes. T is called the SG transi- 

o S 

tion temperature. Finally they obtained the following express- 
ion for zero field susceptibility 


X(T) = Xq (1 - Q) 


(1.27a) 


whe re ^ 

ES (S+1)(gMB) 


(1.28) 


is the susceptibility of free spins. This expression (1,27a) 
exhibit a cuspat T^, as observed experimentally [ 7!^ . But 
the overall shape of the X(T) curve obtained by EA did not 
agree with experimental results. The slope of X(T) turned 
out to be zero just below T contrary to experiment. 

O 

Fischer [ 130] attributed this discrepancy to the quantum 
natuare of spins and developed the theory for quantum spins 
following EA very closely. Assuming validity of the standard 
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result of linear response theory (LRl), Yiz.j 


X(T) = S{[<S^ S^>] i;<Sj_> < S^>]} 

ij av ar 

and definition of q, he obtained 


xyRS(S4-1 ) 


(1.29) 


X(I) = x^[ 1 - (1.30) 

in the absence of SEO. He also showed that = 0.69 for 

1 ® 

S = ^ which agreed well with experimental result (0,6) where- 
as the classical theory of EA predicted a corresponding 
value of 0.92. Moreover in the classical limit (i.e., 
for classical spin). Fischer's results agreed with those of 
EA. Sherrington and Southern [131] developed quantumr- 
mechanical version of EA model, independently and obtained 
results identical with those of Fischer [ 130]. Fischer [132] 
also extended his calculations in the presence of external 
magnetic field. 

But, unfortunately, all was not well with BA model. 

It also predicted [114] a cusp- in the specific heat at the 
same temperature I , which has never been observed experi- 

o 

mentally. This is the first evidence that the SO transition, 
if at all a phase transition, is an unconventional one. 
Recently Barnes [13^ showed instability of the EA ground 
state with respect to magnons. 

Edwards [134] examined the condition for existence 
of SG phase by quantizing a SG consisting of randomly 
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arranged dipoles of moments of inertia I and experiencing 

random forces. SG phase was shown to exist prowided the 

J I 

dimensionless parameter is greater than a critical value 

h^ 

where the classical ground state energy is N being the 

number of dipoles. However, physically the problem turned 
out to be similar to Anderson localization [13^ • The latter 
aspect will be discussed in detail later. 

Physicists felt uneasy at the use of "replica trick" 
because analytical continuation upto n = 0 involves fractional 
values of n also and the physical significance of fractional 
replica v/as not clear. Sherrington [136] and Southern [ 137] 
derived expressions similar to those of Pischer using sinrple 
mean field theory in much more transparent manner without 
using replica trick. 

The EA model has been extended to include uniform 
uniaxial anisotropy [138] and also random uniaxial anisotropy 
[ 139] * Phe latter exhibits a broadened maximum in the 
specific heat instead of a sharp cusp in the former but the 
broadened maximum occurs at almost the same temperature as 
the cusp. 

Step IV ; 

Sherrington and Kirkpatrick (SK) [ 140,141] applied 
the EA theory to an Ising model where each of the spins are 
coupled to all others by infinite ranged infinitesimally 
weak exchange interactions. In other words, in this model 



4y 


the number of nearest neighbours Z is infinity, SK used 
replica trick and Gaussian distribution for exchange inter- 
actions 






r 

exp L - — ^ 


2 J' 


7-^1 


where 

Jo = Jo/*' 


and 


J J/ 1 


1/2 


CO CNJ M 

with Jq and J both intensive. The factors M and Lr^ woi^ 
introduced to make the thermodynamic limits physically 
sensible. Y/hile calculating the free energy, SK inter- 
changed the two limits K ** ooand n 0. The equilibrium 
values of the order parameters M(T,H) and Q(T,H) are given 
by two coupled integral equations, called SK equations. 

However, SK model suffers from two serious drawbacks 


(i) it gives negative entropy at sufficiently lo?/. temper- 
atures. The authors themselves speculated that "this 
unphysical behaviour has its origin in the interchange 
of limits K •> ooand n -* 0, but the consequences are 
confined to low temperatures". 

(ii) the solution obtained does not correspond to a minimum 
of the free energy though it is an extremum. 


Van Hemmen and Palmer [142], on the other hand, showed that 
the unphysical results appear because of nonuniqueness of 



the limit n *♦0, 

Since SK's original derivation is complicated we 
shall sketch the sin?)ler derivation of Mookerjee [ 143] « 

This method has the advantage that no replica trick is 
required because the configuration averaging is done at the 
end of the calculation. The configuration dependent free 
energy is given by 



(1.34) 
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in thermodynamic limit, where 

/* -ikK(r) tanh(6h) 

E(k) = 1 idh P(h) Cl - e ” } (1.35) 

r 

Equations (n. 34) and (1.35) together constitute self consistent 

equations for P(h). Expanding and retaining only terms up to 
2 

the order of K one gets 


E(k) =jdh P(h) [ik Kq tanh0h)+ ^ k^ tanh^( p h)] (1.3 


where 

Kq = s z(a) 

E 

= S K^(R) 

R 

Defining, 

M = f dh P(h) tanh ( ph) 
Q = \ dh P(h) tanh^( r h) 


(1.37a) 

(1.37b) 

( 1 . 38 ) 

(1.39) 


One finally gets the local field distribution P(h) to be 
Gaussian 


P(h) = 




(h - c K M)^ 
exp [ ^ ] 


2 cQ K 


(1.40) 


1 


Hence, substituting (1.40) into (i.38) and (1.59,) we.get,^, ^ 


t ] I' , 


83813 
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00 

J 

-00 

T M t/Q X 

e^/2 tanh(-|-+-f- 

+ f )dx 

(1.41) 

and 






00 

r' 

\ 

-00 

-ic2/o o 

e /2 ^ c ^ 

|)dx 

(1.42) 

where , 






a Pgff iJ-gH/kg, (1.43) 

^eff effective Bohr magneton number. 

The two coupled integral equations (1.41) and (1.42) 


are solved selfconsistently to get the equilibrium values 
of M and Q, This analysis led to the following conclusions: 

(a) three phases-paramagnetic (M=0, Q=0), spin glass 

(M=0, Q;^0) and ferromagnetic Q^O) appear in the 

phase diagram (Bigure 12), 

(b) all the transition lines in the absence of external 
field are second order, 

(c) specific heat is zero at T = 0 as expected from 
the rrao dynamic s , 

(d) specific heat displays a cusp for transition to SG- 
but exhibits a discontinuity for transition to ferro- 
magnetic phase. 

The temperature and field dependence of M and Q in 
the SG- phase have been studied by solving equations (U41) 
and (1.42) numeric ally for AuBe (6.6 a.t.%) and CuMn (0.7 at.^) 
by Mookerjee and Ghowdhury [144] • Klein [145], Kaneyoshi 
[146] and Plefka [147 ] gave alternative derivations of 
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equations (<.41 ) and (1.42). 

We want to stress that although the final equations 
(t.41)a'ad (1.42) are the same in both Mookerjee's and SK appro- 
aches there are some fundamental differences between the two 
approaches* iPirst of all, Itookerjee's model is based on EKKY 
interactions whereas SK used infinitely weak gaussian distri- 
buted infinite range interaction, though both are infinite 
ranged. Secondly, Mookeroee's approach is not based on any 
replica trick and so the interchange of the two limits 
K —r. 00 and n 0 is not required. The configuration averag- 
ing is done at the end of the calculation. The free energy 
is a true minimum with respect to single spinflips as 
computed before configuration averaging. But the stability 
analysis against arbitrary numbers of simultaneous spinflips 
has not been investigated. The entropy before configuration- 
averaging is positive definite and hence expected to remain 


> 1/2 


so even after averaging. 

All the above MBT predict T^ or scaling 

laws support T a' ° whereas experimentally observed exponent 

s 

varies between wide ranges. Explanations for deviation from 
the MF predictions were given by Larsen [148] . He took 
account of the facts that the conduction electrons responsible 
for a-d interaction may have finite mean free path due to 


(i) self damping, i.e., scattering by the localized spins 
( X oc "d) 

scattering by other defects ( X independent of concen- 
tration) . 
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Besides, the fluctuation in the nearest neighbour distance 
also adds to the deviation from &IBT prediction. 

It is vexy interesting to note that Japanese groups 
worked out, independently, and almost simultaneously with 
EA, the B/iET of glass- like phase (GLP) [ 149*-'i52] and random 
ordered phase (ROB) [ 155-156] for short as well as long range 
interactions. Both G-IP and ROP later tumed out to be 
eq.uiTalent to MPT of Sd. Takano [1571 and Sakata et.al. [152 ] 
compared SbP with ROP. SLP is characterized by [<Sj_>] ^^=0 
[<S^> ^ 7 ^ 0 whereas ROP is chara.cterized by 

[<Si> 3^^ = 0 aid [ I<S^>1 /^O and hence are qualita- 

tively equivalc^nt. The differences in the quantitative 
results is due to the different matheuoatical approximations 

adopted. 

Be Almedia and Thouless [158] first showed that the 
SK solution does not correspond to the minimum of the free 
energy and they concluded that the permutation symmetry 
between replicas must be broken in the ’’SO phase". The 
Almedia-Thouless instability turned out to be a special case 

of Pytte-Rudnick instability [159]. 

So, one feels uneasy at this point because the 
authors of references [ 14G],[ 142] and [158] blame different 
aspects of the SK model for the disease. Moreover,Kosterlitz 
et.al. [160 ] showed that for the spherical model identical 
results are obtained with and without use of replica tnck. 



Khurana and Hertz in an alternative (mathe- 

matically equivalent) diagrammatic approach, without using 
replica trick, showed that this instability arises due to 
the insufficiency of the order parameter Q to describe the 
SG- below Ig. In this approach, the instability is signalled 
by a negative value of the order parameter susceptibility 
D below Tg, where D is defined by 

Q = D h^, 

2 

h being the mean square value of random field h.. For 
stability D must be positive definite. Unlike the piascri- 
ption of breaking the replica symmetry, Sompolinsky [ 117] , 
Hertz et.al. [161] and Hertz [162] suggested dynamics as a 
cure- for this instability. ¥e shall first discuss replica 
ssrmmetiy breaking schemes and then the dynamical approach.. 

Blandin et.al. [163], Paris! [119-121 ] and Bray 
and Moore [ 164-166] suggested different schemes of replica- 
symmetry breaking. Here we shall discuss Parisi^s scheme 
which is the most popular among all, 

Paris! [119-121, I 67 -I 71 ] suggested the following 
parame t rization 

Q =0 

Q = if I ( a /mj_) ^ I( P/m^) and 

I ( a = I OAi+P 

where I(x) is an integer valued function, i. e, I(x) is 
the smallest integer equal to or greater than x. Q^(i=0,K) 



are real numbers and m. (i= 1 ,K) are integral numbers such 
m. 

that is an integer (i>1). Then putting in the basic 

i 

requirements (i) that the configuration-averaged free energy 
be finite, (ii) that it be a minimum with respect to Q and 
(iii) that its minimum at high temperature be located at 
Q = O 5 he obtained 

1 

X(T) = 5 ^ (l-STb) 

o 

after taking the limit K ^ oo so that the order parameter 
Q(x) is now a function of x defined over the interval [ 0,^0 • 
Prom numerical calculations for K = 1j2 etc*, and then 
extrapolating to K '»■ 00 Parisi [168] obtained 

X 

X(T) = (1 - X) (1 - 5ax] 

° (1.27c) 

whe re , 

Q(x) is linear in x for x < x, and 

QM = «Bax = %IA ^ ^ 

This abstract scheme has been explained pictorially by 

Rammal and Soule tie [15] . 

There is, of course, lot of arbitrariness in this 
scheme so far as the choice of the values for the various 
parameters are concerned. However, Parisi [168] showed 
that entropy becomes less and less negative as K is increased 
from zero (no symmetry breaking) to higher and higher values 
(severe symmetry breaking). Therefore, one expects that the 
best possible choice for the parameters should correspond to 



the most severe S 5 aimietry breaking among the replicas. Thouless 
et.al. [172] and De Dominicis and Kondor [ 173] showed that 
Pari si's solution is marginally stable. 

Though the meaning of x was not clear from the original 
formulation, Gabay and Garel [174] interpreted it in the light 
of two fluid model of superfluidity as folloY/s; 'a fraction 
X (T) of the spins has condensed and a fraction 1-x(T) remains 
normal. This condensed state spans all x (0 < x <x), so that 
X appears as the probability that a given spin has ai BA order 
parameter lying in the interval (0,Q(x))". The condensed part 
consists of overblocked plaquettes, forming 'walls' , whereas 
the normal part are frustrated forming 'domains' * In this 
light "SK solution can be viewed as a solution without walls". 
The hint for this interpretation came from the fact that 
X = -^ at T = 0 predicted by Parisi Toulouse hypothesis, to be 
discussed later. Parisi' s order parameter Q(x) , as discussed 
above, is a function of one variable x(0 < x £1). Recently 
Jonsson [175] ' presented a scheme whereby the order parameter 
be.comes a function on a unit square, i.e., q(x,y), (0<x<1, 

0<y <1). However, physical interpretation of x and y is 
not yet clear, later on Sompolinsly [117] related x to the 
spectrum of relaxation times. We shall discuss Sompolinsky' s 
approach later in this article. 

Very recently Parisi [176] and Houghton et.al. [ 177] 
have suggested interpretation of x and Q(x) in terms of the 
probability distribution of the overlap of the magnetization 
in different pure states of the system. 
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Now let us discuss the altexnative (dynamic) approach 
suggested by Hertz et.al. [ 161] in order to remove the instabi- 
li1y. Using a diagrammatic expansion, they showed that in 
the limit w 0 the self consistent equation for the dynamical 
response function X(w) has two solutions - the solution 
consistent with fluctuation-dissipation theorem is the same 
as the usual equal time spin- spin correlation function of the 
static formalism and leads to instability^ the second solution, 
on the other hand, though violates fluctuation-dissipation 
theorem (HDT) stabilizes the ground state against order para- 
meter fluctuations and yields a finite positive value of the 
oirder parameter susceptibility. 

Next we shall discuss Sompolinslsy’ s approach and then 
correlate the works of various authors as proposed by Hertz [16i§ 

Sompolinsky [ 117 ] investigated the static (statics 
as a limit of dynamics) properties of SG using a dynamic 
approach. He assumed that there exist a distribution of 
relaxation times, instead of a single one, all of which become 
infinite in the thermodynamic limit. In order to correlate 
with Parisi’s results, he parametrized the relaxation, times 
t in a decreasing order by an arbitrary paranBiter 0 <x <1 
such that 

Q(x) = [<S^(0) av (1.16b) 

measures the correlation retained at time scale t^. The 
maximum value 


dd) = %,A 
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The most interesting result was that equilibrium Talue ,x(0) 
of susceptibility was given by 

X(0) = p[ ^ ^ ^ 

whereas the quasi equilibrium value X (1) after a time large 
compared to microscopic times but small with respect to 
macroscopic time scales was given by Jischer’s relation X(1) = 
p (1 - a(x) is an additional order parameter: 

A(x) > 0 as X >• 1 

Also Q(x) 0 as X 0. 

But a peculiar feature was that "since the scale of x is 
arbitrary the theory can only establish a relation between 
Q(x) and A(x) but cannot determine both of the functions, 
except at the uniquely defined end points 0 and 1" [117]. 

Sompolinsky free energy functional was rederived by 
de Dominicis et.al. [ 178] by the replica approach in a static 
treatment. This was another indication of a possible close 
relation between Paris! ’s replica symmetry breaking scheme 
and Sompolinsky ’ s dynamic approach. Purther support for such 
close relation came from the observation [179] that the free 
energy could be written as the solution of a differential 
equation-one of the special oases, viz. 

A (x) = - xQ ’ 

yields Paris! 's differential equation [168] . This relation 
was also obtained by Sompolinsky, because as already stated, 
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his theory could not determine Q and A independently but 
established only the relation between them. 

Hertz [162] showed later that the SK spin glass is 
only in dynamic equilibrium. He developed the dynamical 
theory for a finite system. Since no finite system can 
exhibit a phase transition, he assigned a finite width s 
to the 6 -function part in the spin correlation function 
C(w). The spin correlations die out for t>> e ( s-»0 as 
the size of the system -*• co ) and susceptibility becomes 
identical with that of free paramagnetic spins. Thus Almedia- 
Thouless instability is suppressed only for tiuB scales 
t « e"""' . Of course, for infinite systems, which are relevant 
for comparison with experiments, this theory remains dynamica- 
lly stable for all times because e -» 0. Nevertheless, no 
stable equilibrium description is provided. Hertz interpreted 
this in the light of the large number of metastable states in 
SG and will be discussed later. 

In addition, assuming that equilibration proceeds in a 
stepwise manner — successively larger and larger spin clusters 
flipping at successively longer and longer time scales — 

Hertz was able to interpret Parisi's Q(x) as "a measure of 
the local persistent spin autocorrelation for a metastable 
state where equilibrium has been reached with respect to all 
flips of spin clusters upto a size (1-x) N, but not for flips 
of larger clusters” [ 18#] « 
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Although, the mathematical structure of this theory 
looks identical with that of Sompolinsky, there is a crucial 
difference — no explicit adhoc violation of fluctuation- 
dissipation theorem is necessary though such a violation was 
concluded by Sompolinsky [117] j and by Hertz et.al. [ 161] , 

Recently Bray [ 181] examined the effect of including 
a site-dependent magnetic field proportional to the eigen- 
vector with largest eigenvalue of the random matrix and 
concluded that it only suppresses strongly, but cannot remove 
Almedia-Thouless instability. 

Parisi- Toulouse (PaT) Hypothesis ; 

PaT hypothesis [182-184] (also called projection 
hypothesis) states that in the SG phase the entropy does not 
vary with the magnetic field, i.e., 

S(T,H) = S(T) 

implying that equilibrium magnetization 

M(T,H) = M(H) 

which together with the hypothesis 

Q(T,H) = Q(T) 
leads to P(T,H) = P.j(T) + 

The predictions of this hypothesis agree well with results of 
monte-carlo simulations [ 184] . If extrapolation to infinite 
system is' allowed, Ma and Payne’s [185] numerical computation 
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of entropy on a finite system with, long ranged interactions 
justifies PaT hypothesis. However, PaT hypothesis has been 
shown to be approximate [ 186} « 

nonlinear Susceptibilities i 

If paramagnet-to-SG is a genuine equilibrium phase 
transition then X should diverge at T with a universal 

^ o 

critical exponent. However, observed variation of this 
exponent with temperature led Pahnle and Egami [187]to the 
conclusion that the divergence of Xg "has nothing to do 
with a real critical phenomenon". We shall not discuss the 
nonlinear susceptibilities any more in this thesis. 

Improved Mean Field Theories ; 

The Be the approximation for short range site random 
as well as bond random I sing model has been carried out by 
Katsura et.al. [188,189] , Ghatak and Moorjani [ 190] , Klein 
et.al. [ 191] jMorita [192] and Horiguchi [ 19^ '. An attempt 
to cure the diseases of the SK model was made by Thouless 
et.al. (TAP) [194] by applying a diagrammatic method and 
without using the replica trick. They observed that "it is 
simpler, and perhaps more physical, to observe that the set 
of diagrams contributing in order H are just those which 
would remain on a "Be the lattice" or Cayley tree, all 
diagrams with no loops. Any diagram containing a closed ring 
is necessarily of order H/Z or less, since the internal 
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connection reduces the number of site summations by one". 
Moreover the I sing model on Be the lattice is in principle 
exactly soluble by the so called *Bethe' or ^ cluster’ or 
Bethe-Pierls-feiss' (BIW) method. The final equations for 
local magnetization, known as TAP equation, is 

= tanh ( p h^ (1.44) 

where 

h. = E J. . M - p. (1.45a) 

0 

The physical interpretation of equation ( 1 .4^) is veiy clear - 
it subtracts out the Onsagar reaction field from the local 
magnetic field acting on a spin* This can be seen as follows; 
the field created by ith site at the jth is and hence 

induces a magnetization where is the local 

magnetic susceptibility at the jth site. This induced moment 
produces a reaction field J?. M. . at the site i. Summing 

X J -L J J 

2 

over all 3 ' -s, the total reaction field at site i is S 

X • . M. and hence the reaction field = p E J? . ( 1 IL t 

33 1 ^ j 

because, from IRT 


Xjj = p (1-1^) 


In the case of mr- component vector spins [195-196] equation 
( ) geae rail ze s to 






(1.45b) 
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ITakanishi [197] developed two-spin and three- spin cluster 
theories of SG. The equation for local loagnetization for 
two-spin cluster was exactly identical with TAP equation 
because the mathematical approximations made in both aare 
equivalent. However, the result of three-spin cluster theory 
was an expected generalization of equation (1.45): 


h- 


SJ 






- SJ?. 
13 


X. 
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% - 


J 


X. X, 


io "3k ki "j "k 


T''". 


(1.45c) 


This subtraction of the reaction field is very crucial 
for SG as shown by Gyrot [198] because the magnitudes of the 
three fields - Weiss field, Onsagar field, and the cavity 
field are of comparable magnitude in SG. In fact the reaction 
field can be larger than cavity field in these systems. The 
physical importance of Onsagar field was, of course, realized 
much earlier when Brout and Thomas [ 199] showed that the 
apparent contradiction and violation of fluctuation theorem 
in naive MPT can be removed by Onsagar' s method. 

Next we come to the crucial evaluation of TAP model - 
what have been gained by this model? Do they really cure 
the diseases of the SK model? TAP solution is a saddle point 
of P and not a minimum [11]. On the other hand very recently 
Plefka [ 200] has shown that the entropy in the TAP model is 
indeed non-negative . He also obtained an upper bound for 
the entropy. So far as the crucial question of the occur3rence 
of the phase transition is concerned TAP asserts sharp phase 
transition like SK model. Plefka [ 201 ] obtained the following 
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lower bound for EA order parameter Q 

p - I t^ for 0 < t < 2/3 

0 > < 

1 - t for 2/3 < t (1.46) 

where , 

t — CpJ) ~ fp 

c 

/V 

Since q = 0 is not allowed for t < 1 , SG seems to be a 
genuine thermodynamic phase distinct from paramagnetic and 
ferromagnetic phases. Moreover, Pleffca’ s result agrees with 
TAP's original result: 

Q = 1 - t for t < 1 

= 1 -■ a t^ for t 0 (1.47) 

The slight disagreement in the value of a ( = 1.665 in TAP 
and 1.5 in Plefka) affects the TAP results by little amounts 
which fall within the error bars of computer simulation. Also 
notice that all these apparent successes (even if partial) 
have been achieved without using replica tricM. 

TAP condition for convergence of their expansion 

J ( 1 - M^) < T (1 *48) 

was valid only near T . later Plefka [118] presented a 

o 

general convergence condition valid for all T: 


(p J)'"^> max {(1 - 2Q2 + Q4)? 2(Q2 “ Q4)} 


(1.49) 
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Howe-vfer recently Owen [202] has corrected some errors in 
Plefka's calculation [118] and showed that the condition 

mmm O 

( p J) > (1 •* 2 Q 2 + Q^) is both necessary and sufficient 
for all temperatures. 

But TAP model has its own draw backs. It predicts 
a cusp in the specific heat exactly similar to the SK model. 
Besides, in the finite Z analogue of the TAP model [ 203] specific 
heat becomes negative near T = 0 because of negative slope of 
entropy, for sufficiently large values of Z, although entropy 
at T = 0 is zero, fohlfarth [204] argues that since 


H(-|^)h = (4f)H = 


slopes of M-T and X-T curves must approach zero as T 0 
if third law of thermodynamics is to be obe 3 ?ed, because the 
latter claims that "as the temperature tends to zero, the 
magnitude of the entropy change in any reversible process 
tends to zero" [205]. TAP clearly violates this law. 
Wohlfarth accuses the relation ( 1 ,47) as the possible cause 
of this anomaly. So far as the stability of TAP solutions 


are concerned Bray and Moore [206 ] computed the eigenvalues 
of the Hessian matrix for each solution of the TAP 

equation. Por stability all the eigenvalues of the Hessian 


matrix must be positive definite. Though, all the eigenvalues, 


indeed, turned out to be positive? the eigenvalue spectrum 
extended down zero for all T < T^ implying that the solution 
is always on the limit of stability. 
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Let us now compare and constrast the SK solution and 
the TAP solution of the infinite ranged Ising SG model. 

Sommers [207], using a diagrammatic expansion of the free energy, 
summing the relevant diagrams and averaging over the ring 
diagrams and tree diagrams at successive stages, showed that 

(i) the free energy resulting from averaging over only the 
ring diagrams is identical with the TAP free energy, 

(ii) averaging over both ring as well as tree diagrams 
indicates the existence of two solutions one of which 
is the SK solution. The other solution does not suffer 
from the unphysical properties of the SK solution. 

Sommers' theory predicts that the entropy vanishes 
as 

exp [-2 (2A J/T] as T - 0 (1.51) 

2 

in marked contrast to S '=»= T as T -* 0 as shown by TAP and 
confirmed by numerical experiments. Bray and Moore [ 16^ > 
in addition to the above point, raised a serious objection 
against Sommer's solution; "an exponential decay suggests 
a unique ground state separated by a gap from the first ex 
excited state, whereas one expects on physical grounds a 
quasi- continuum of nearly degenerate ground states". 

De Dominicis [208] from a diagrammatic expansion of Martin- 
Siggle-Eose type, showed how to recover SK and Sommers' 
solutions in two special cases. Garel [ 209] pointed out the 
analogy between the properties of the infinite range Ising 
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SG and the spinodal behaviour of a ferromagnet. In case of 
a ferromagnetic I sing model, the system can ;iump from the 
unstable branch (SS’) to the stable branch (See Pigure 13) 
whereas in the case of an I sing SG "the system is always on 
the verge of 'jumping’ but that there is no thermodynamically 
stable state onto which the system can jump because of frus- 
tration. Ihe system therefore chooses the lowest available 
energy state, neither stable, nor unstable, that is a typical 
spinodal situation". Plotting q, against its conjugate field 
H' (see Figure 14) we see that the curve is similar to the M 
vs H curve in Figure 13 except for the fact that the system 
does not get stable branch to jump on in the former. Q < 
are not allowed as was originally shown by TAP [194]. oK 
solution falls in this forbidden region. 

From the simple consideration of the temperature 
dependence of the order parameters 

Sap ^ f s M? 

1 ~ 

«-ta = S’ Slin Sj_(t t) 

t - 00 

and Q =[<S->^] , Young showed [ 210] 

av 

(a) the SK model is non ergodic, 

(b) a time average gives the same result as a single solution 
of the TAP equation, 




figure 13 t The ferroiaegjietlo Islng model in a magnetic field 
(T<,X ), The branob SS* is thermodynamioally 
unstable, (after ref, [209J ) • 
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(c) '‘ensemble averages correspond to averages over many 
TAP solutions", and hence, "interfeience between 
different solutions causes the ensemble average to 
decrease more rapidly with T than the results of one 
solution" . 


In short, the non-ergodic behaviour arises because of 
the splitting of phase space into many separate configurations. 
This behaviour is not a characteristic of only spin glasses 
but of any glass transition e.g. liquid-glass transition. 

Jackie [ 211] has compared this ergodic-to-nonergodic transition 
for glasses with that of SG-. 

Another similar and related question is: Are the 
famous Maxwell's thermodynamic relations violated in SG? 

This is a very controversial question. Pogle et.al.'s [212] 
calorimetric measurements apparently indicated such a failure. 
They calculated susceptibility from their measured values of 
G(H,T) using Maxwell's relation 


{ 


3 





_ 1 

“ T 8H^T 


( 1 . 52 ) 


The calculated values of susceptibility did not agree with 
the experimentally measured values. But a deeper analysis 
[ 213-21^ revealed that their data were consistent with 
Maxwell's relations. On the other hand, from an extensive 
numerical calculation Soukoulis et.al.[ 217, 21^ concluded 
that Maxwell's relations were violated in all cases except 
for the field cooled samples. Of Qourse, the question 
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remains unsettled. 

Kaneyoshi [219] introduced a new method of eyaluating 
the SG transition temperature T beyond its mean field value 
by using the method of differential operator and carrying 
out a new type of expansion for the order parameter q. near T . 

O 

A full fledged cluster model of SG was developed by 
levin, Souhoulis and Grest (1SG)[ 220-22^. The Hamiltonian 
was divided into inter- cluster and intracluster parts as 
follows 

=- l.<x . Sx + ^ S.^ . 

V i 

where represents the spin the ith site in the a;th 

cluster.,!?^ and are the intra cluster and inter cluster 

exchange interactions respectively, 

A very remarkable fact was that the moment per cluster 

turned out to be temperature dependent which in turn implied 

a temperature dependent size of the clusters. An equivalent 

in 

model will be treated later/this thesis where the volume of 
each cluster will be temperature dependent, ISG claimed that 
this model can explain almost all the experimental results - 
susceptibility, sp. heat, resistivity and neutron scattering 
experiment. However, the genesis of these clusters, which 
has to be universal in order to explain universal properties 
of all SG, is not convincin^y clear. Besides Murani [223] 
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criticized Soukoulis et.al.’s works on neutron scattering. 
Soukoulis et.al. subtracted a suitably chosen static compo- 
nent from the total scattering cross section; the static 
component decreasing mono tonic ally from its maximum value 
at T = 0 to zero at T = Tg* The residual (i.e., dynamic) 
part exhibits peak at unique Tg independent of q. However, 
Murani argued that "the resulting l-inde pendent peaks in 
the dynamic component are 'artifically' created additional 
anomalies in as much as they are assumed to be there to 
begin with". Such peaks can be obtained even if there were 
no peaks originally in the total cross section as shown in 
^gure 15. Moreover, the sharp change of slope at 1 in the 

O 

total scattering cross section calculated by Soukoulis et.al, 
has never been observed experimentally, a reply to this 
criticism Levin et.al. [ 223] attributed the latter discre- 
pancy to their random phase approximation and not the physics 
of their theory. But their answer to other objections of 
]V[uT^a.ni seem to be less convincing to us. 

It is interesting to see that Morandi ard Gorbelli 
[ 224] argued that "no transition with a continuously vanishing 
order parameter is possible for any physical value of the 
concentration without violating some necessary conditions 
on the structure of the internal field distiibution function" . 
They attributed the incorrect conclusions of theories of 
SecovtdL kind to neglecting the spatial fluctuations. HoTvever, 
from a detailed study of the local field distribution, and, 



-I ^ 



15 * SoheHuatlo Tariatioa of th» lateaBl'^ ^th temp®rftttt3C« 

, . for one q value, Jhe chain curve shows a pos^hl* 

Biono tonic fom for the total steasured inl^nsiljj'’ (note 
' no sharp peak originally), Jhe solid ourro represent# 

the static (ot Bragg) contribution with its sharp 
discontinuity at Tgn which is circled* The resultant 
dynaoBiio conq^onent is. given by the dashed curve toother 
with tbB chain curve vatjove * 0 ( 5 ; )• ^e blear that 
the sharp peak formed at the junction of the two curves 

ct results of the assumed sharp 
(also circled) in the static 
assumed to be the same for all q 

values, it is evident that the dynamio component will 
have sharp maxima at the saiM temperature for «Q1 

q (in addition to asqy other sharp fea-hires present 
in the total intensity curves), (after ref* [223}). 


(eircled) is the dire 
discontinuity at Tga 
component. If Tesn lo 
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hence, the upper and lower bounds of the transition temperature 
Gorbelli and Morandi [ 225] concluded that a "weak first order 
transition” takes place at a temperature 



because ^(T) vanishes discontinuously at 

o 

Step V ; 

Several attempts to understand the critical behaviour 
of the SG- transition have been made using the conventional 
techniques only with partial success. Most of these results 
were mutually contradictory and, hence, inconclusive. Since 
in this review we are mainly concerned with only the MPT we 
shall not discuss these attempts. The best references are 
the review articles [8] ,[11], [18], The upper critical 
dimensionality seems to be six whereas the lower critical 
dimensionality four, though these points remain highly 
controversial. One of the interesting results was obtained 
recently [ 226 ] applying the ideas of hierarchical Ising 
model. It was observed that the renormalization group (RG) 
trajectory is chaotic. By chaotic EG trajectory one means that 
"as successively longer length scales are viewed using the 
EG, strong and weak correlations are encountered in a 
chaotic sequence". The other interesting result is that 
the upper critical dimensionality is eight for m-vector SG 
in the limit m oo [227] • Besides Morgan-Pond [ 228] has , 
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recently, sbown that the size of the critical region for long 
range (but not infinite) Ising spin glasses is proportional 
to 1/Z where Z- is the number of interacting neighbours. When 
compared with phase transitions in other magnetic systems 
this result is unusual, 

1.5 Frustration in Spin Glasses 

G-ene rally, frustration arises in SS due to the contra- 
diction of the interactions. In canonical SG, for example, 
due to the oscillating natU 2 ?e of the EKKY interaction and 
random separation between different spins it may so happen 
that the information reaching at an arbitrary sLte-i- from 
other sites are in conflict with each other leading to frus- 
tration of the ith spin. Since spins connected by frustrated 
bond do not satisfy the lowest energy (ground state) config- 
uration, the energy of the system is always higher than the 
corresponding unfrustrated system. It will be shown in this 
article that frustration distinguishes nontriwial disorder 
from trivival disorder. Let us consider spins on a square 
lattice with nearest neighbour exchange interaction with 
Ij I = 1 , as shown in Figure 16, where + stands for ferro- 
magnetic and - for antiferromagnetic bond, notice that the 
plaque tte in Figure 16a is unfrustrated whereas that in 
Figure 16b is frustrated because all the bonds cannot be 
satisfied simultaneously in the latter. The number of anti- 
ferromagnetic bonds in Figure t6a is even but odd in Figure 16b. 
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In a classic paper [229], Toulouse utilized this fact to 
define frustration function 0 over any closed contour (C) 
of connected bonds; 

5 = I i sign J . . (1.55) 

(C) 

A plaquette is frustrated or unfrustrated depending on whether 
S' = - 1 or + 1. Bethe Lattice does not contain any closed 
ring and hence totally unfrustrated. 

An interesting model of disordered spin system that 
apparently looked like a SG was introduced by Mattis [ 230] . 
This model is defined by 


H = 


S 

id 


Jij Si Sj 


where disorder is introduced via 


= ±1 randomly on sites. However changing the variables 

S. - s! = e- S. we recover an ordered problem. In other 
1 1 11 

words, Mattis model does not contain frustration and hence a 
model of trivial disorder. However, the disorder in Mattis 
model for quantum spins is nontrivial because the transfor- 
mation is not necessarily canonical [ 231-2331 • 

A very important consequence of frustration is the 
degeneracy of the ground state. This can be understood 
geometrically as follows; 



"a frustrated bond belongs to two plaquettes (in two dimen- 
sional models) and such bonds form a set of strings linking 
the frustrated plaquettes. The strings are constrained to go 
through the plaquette centres, and the total length of these 
strings must be minimal to obtain the lowest energy configu- 
ration. There are in general many such minimal pairings, 
their number is the ground state degeneracy" [ 234] . 

Perhaps, the most important impact of the frustration 
model is that it showed a close relation between quantum field 
theory and condensed matter physics (and possibly with many 
other branches of human knowledge) [ 235-242] prowing once 
again the unity of physics. Psychologists hawe studied the 
stable configuration of a group of human beings in the presence 
of attractions and repulsions between them [ 234 ] . The geometry 
in a SG has similarities v/ith the "parallel transport" of a 
tangent vector in a curved surface. The misfit between two 
different lines of transport is expressed by the ‘curvature* 
and 'frustration' in the two cases. In this sense, a frus- 
trated plaquette is 'curved' whereas unfrustrated plaquette 
is 'flat'. The Mattis model is flat. 

The advantage of applying techniques of gauge theory 
to SG is that frustration can be expressed as a monopole [243]. 
Rammal and Soule tie [15] calls frustration "the third concept 
in the physics of disorder" - the other two being the concepts 
of 'percolation' and 'localization'. Percolation concept is 
associated with the 'dilution disorder* whereas ‘localization’ 
describes ‘spectral disorder’. In the same spirit we can say 
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that frustration is associated with "comnuni cation disorder” 
among the various spins. 


1 .6 Dynamical Theories of first Kind ; 


(a) 

(’d) 


There are two aspects of spin dynamics in a SG-; 

Very Low Tenperature (T«T ) or short time behaviour 

O 

and 


Low Temperature (T<T )..Qr long time behaviour. 

& 

Very Low Temperature (T<< T ) or short time behaviour# 


At low temperatures the system gets "trapped" in 
one of its low lying metastable states, at least for a short 
time. Consequently, the corresponding excitations, commonly 
known as spin waves, reflect characteristics of that particular 
metastable state. In other words, these spin wave excitations 
arise as small amplitude collective oscillations of the system 
about the given metastable state when the system responds to 
an external perturbation. Such excitations have been studied 
analytically [244-253] as well as numerically [254-261], 
Halperin and Saslow [244](HS) developed a hydro- 
dynamic theory of spin waves in SG. The fundamental assumption 
of the hydrodynamic theory is that the state of the system with 
small amplitude, long wavelength fluctuation from equilibrium 
is essentially determined by the long-wavelength variations 
of the various relevant quantities, e.g. , magnetization density, 
energy density, etc. In that sense HS theory of spin waves 
is not a microscopic description. The important results of HS 
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treatment are; linear dispersion law, i.e. to oc k and that 
damping is proportional to k^. One must remember that these 
results hold only for k -* 0 as the latter is the underlying 
assumption of the hydrodynamic approximation. Andreev's [2453 
treatment, in the limit of negligibly small external field 
and relativistic interactions reduce to the HS treatment. 

These linear HS modes were also obtained by IDzyaloshinskii 
and Volovik [240]. Saslow [246] extended HS hydrodynamic 
theory to include remanence and anisotropy effects. Saslow 
assumed that the external field H is applied along the same 
direction as the cooling field KL (i»e. the direction of 
remanent magnetization M^). later Henley et.al, [247] removed 
this restriction and took into account a finite angle between 
H and H^. In the HS theory two macroscopic variables descri- 
bing the dynamics were the magnetization m and infinitesimal 
angle © by which the spins are rotated from equilibrium. 
Saslow [248] introduced an orthonormal triad of unit vectors 
(a , Pj i)> called the spin triad in terms of which individual 
spin orientations could be specified. The equilibrium orien- 
tation for (n, p, q) is called the anisotropy triad (H, P, Q) 
[248]. Recently Saslow [249] has shown that in the presence 
of strong uniaxial anisotropy there exist two ground state 

A 

solutions for n (i.e., the direction corresponding to the 
cooling field and hence remanent magnetization) - one in 
which n lies in the H - H. plane, called the planar solution, 
and in the other (n, p, q) is rotated from (H, P, Q) by 180° 
about the axis midway between H and H^, called the 
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^ -rotation state. The state considered by Henley et.al. 

[ 247] was identified with the planar solution. 

So far as the macroscopic description of SG are 
concerned, the latter has similarities with the description 
of amorphous solids by the theory of elasticity. The spatial 
displacement vector, velocity and deformation in the elasticity 
theory are the analogues of the rotation of all spins through 
the same angle, the temporal derivative and the spatial deriva- 
tive of the rotation angle, respectively, for SG, However, one 
remarkable difference is that the different displacements 
commute whereas different rotations, in general, do not commute, 
Takayama [250] gave a microscopic analysis of spin 
waves in SG using the standard technique of Hols tein-Primakof f 
transformation on the same line as Sherrington’ s [251] treatment 
of spin waves in Mattis random magnets. The approximations 
made in this calculation were valid for k<<l/Ei^ where is 
the m»an distance between the spins. 

The linear dispersion law was in agreement with HS, 
but the damping rate turned out to be proportional to k and 
hence much smaller than that predicted by HS. Takayama [ 250] 
also predicted a possible energy gap in the excitation spectrum - 
the gap arises due to the inclusion of magnetic dipole intei^ 
action. However, Takayama *s approach was criticized by Barnes 
[ 252] . One should expect three distinct spin wave branches 
due to the spontaneous breaking of the 0(3) symmetiy whereas 
Takayama obtained only one. Barnes [252] obtained all the 



three branches in the liirit k -* 0 for a three-dimensional 
isotropic SG- with an axial anisotropy energy. The two trans- 
verse branches have identical gaps while the longitudinal does 
not. The effect of an external magnetic field is to split 
further the two transverse branches. In general, the disper- 
sion is linear for the longitudinal branch. The quadratic 
dispersion for the two transverse branches reduce to linear 
dispersion law for vanishingly small anisotropy. The gap 
between the macroscopic theories of Andreev [245] and HS 
[ 244 ] and the microscopic theory was bridged by Ginzburg [253] 
whose microscopic equations reduced to the corresponding 
hydrodynamic ones in the appropriate limit, 

Ohing, Huber and coworkers [254-257] and Krey [258] 
carried out extensive numerical calculations on the existence 
of spin wave excitations in SG and dynamic structure factor. 
Their main conclusions weie that, for low enough k, the modes 
were delocalized for all SG, as predicted by the analytical 
theories with a remarkable contrast between the metallic SG 
(e.g. CuMn) and the insulating ones (e.g. Eu^ Sr^^^jin the 
case of the latter the modes are delocalized over nearly the 
entire band whereas for the former only low-energy modes are 
delocalized. Similar results were also obtained by Walker and 
Walstedt [259-26(] whose main quantity of interest in this 
connection were the density of these low-lying excited states 
for SG with EKKY interaction and the effect, on these spectra, 
of the damping of BKKI interaction, dipolar anisotropy, cubic 
■crystal field ^isotropy and uniform external field. The 
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value of the stiffness constant for RKKY-type SG was found to 
4/5 

vary as c and about 25 times smaller than the corresponding 
value for a comparable ferromagnet [261], consistent with the 
predictions of HS, One of the assumptions of Watstedt [ 26 1] , 
namely, that the exchange stiffness energy varies linearly 
with area was later confirmed by Banavar and Oieplak [262 ] 
thereby lending further support in favour of the existence of 
spin waves in SG. 


low Temperature (T< T ) or long Time Behaviour 

- - - - — ^ - - - -■ - ■ - -- 


Thermally activated or quantum mechanical tunnelling 
commonly known as spin relaxation becomes predominant at higher 
temperatures or long time scale. These phenomena arise due 
to the flipping of single spins or clusters of spins as a 
result of which the system 'hops’ from one metastable state to 
another or to the global ground state, Bor the time being, 
we shall be concerned only with the relaxation of the order 
parameters in SG. 

The simplest approach is the use of langevin equation 
together with a phenomenological Ginzburg- landau free energy 
functional —the so-called time dependent Ginz burg-landau (TDGl) 
approach. In case of SG there is no universally accepted form 
of the free energy functional. However, the T^ free energy 
is believed to be good enough for this purpose. Kumar and 
Barma [263] first adopted this procedure where the free energy 
functional near fg is given by 
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5 (Q) = J [ I + I e + I ] (1.56) 

where 



ipo 

and J is the r.m.s. width of an assumed gaussian distribution 

of random exchange integrals. The important point to note 

is that both odd and even powers of q, appear whereas for 

ferromagnets, from symmetry considerations, only even powers 

of the order parameter M appear. This is because dimensionally 

2 

Q is equivalent to pi , Por temperature less than, but very 
close to, T Kumar and Barma observed power law decay; 

s 

Q t“''/^, (T "Tg) 

whereas for T > T the decay is exponential: 

o 

Q ^ exp (- e^tA ) ) (T > "Jg) 

Kumar and Barma used TAP free energy where M = 0. Therefore, 
they could not derive any expression for magnetization relax- 
ation. Later Shastry and Shenoy [264] generalized the free 
energy expression to include terms involving even powers of K 
as well as terms that couple M with Q: 


where SF- is Kumar and Burma's [ 263] free energy (1.56) 


^ V = 
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and 

*liQ = - b m2 Q - 0 

In another TDGL approach Ma and Rudnick [265] rederiTCd the 
-' 1/2 

t behaviour of the spin correlation function. Some 

fundanBntal errors of the latter treatment were later pointed 

out by Sherrington [266], Besides, Ma and Rudnick’ s model 

[ 265] does not contain any frustration in explicit manner. 

Hence it was not clear how such a model can represent a SG 

phase. Hertz and Klemm [ 267] removed this deficiency of 

such a model considering short ranged random exchange inter- 
— 1 /2 

actions. The t tail, was recovered also in this model at 

and very near T^ for T < T . 

8 8 

Hext better approach is the Master equation approach 
which was first applied to one dimensional Ising model by 
Glauber [268] and hence popularly known as Glauber model 
(also called Kinetic Ising Model), Application of Glauber 
model to SG was first carried out by Kinzel and Rischer [269 ] 
in very simple fashion, which yields exponential decay of 
the self correlation function for T>T , and power law 
(t ) decay at T = T . This method was generalized for 

O 

any spin S by Riess et.al. [ 270 ] who observed exponential 
decay both above and below T and hence much faster than 
the observed power-law or logarithmic decay. They attributed 
this discrepancy to the assumption of a single common 
relaxation time for all spins. In reality, different spins 
can have different relaxation times. 
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Kiazel and K-scher's [269] theory was extended by Mody 
and Rangwala [271] to include cluster effects. Bray et.al. 

[ 272 ] used the same Glauber model to study the average auto- 
correlation function 


a(t) = ^ 

and EA order parameter 

Q(t) = ^ 


s < Si (0) s^(t) > 

i 


■’ t 

O 


which are related by the expression 


Q(t) = f dt' (t-f) a(t») 

o 

They derived an equation of the type 


+ 1)<Sj_ (0) Sj, (t)> = <Sj_(0) Sj_(t) tanh(ph^(0)) Sj_(0)> 


(1.57) 


and then irade the crucial approximation that the same equation 
holds even vd.thout ensemble average. At low temperatures, 
the autocorrelation function decays logarithmically whereas 
the EA order parameter Q(t) tends to zero for long tine t in 
agieement with Monte-Carlo results. Bray et.al. [272] approach 
is probably the only one that yields logarithmic decay even 
in the homogeneous MEA. It is not clear whether it is 
consequence of the unconventional approximation used. Very 
recently a novel idea of Sompolinsky and Zippelius [273] led 
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to a power law decay of the form t”^ for the autocorrelations 

-1 

where the exponent 2^ 2 depends on tenperature as 


V 


(T) = 4 - 


JL 

It 


(1 - 4-) + o((i - 


1. 

T. 


)2) 


(1.58) 


A closed expression for the generating functional for 
the dynamic correlation function was derired by Schuster [274] 
using De Domincis [17 ] method based on an extension of Martin- 
Siggia-Eose approach [ 275] . A close relation between SG- and 
many other systems exhibiting l/f noise (e.g, semiconductors, 
dielectrics etc.) [276,277] has been indicated by Kogan [278] ; 
the kinetic I sing model was argued to be fairly good descri- 
ption of all these systems. 

There is a general trend in the literature to explain 
the slow decay as a result of spin dynamics of the kind 
discussed in article 1,6b. Howe-ver Bray and Moore [ 279] claimed 
that such slow decay processes can be dealt with the spin 
dynamics of the kind described in article i .6a. Besides, Saslow 
[ 280] included both diffusive modes (discussed in article t 6b.) 
and relaxational modes (discussed in article 1.Q>)in the sane 
frame work. For k -» 0, rather than having two propagating 
spin waves, he obtained one diffusive mode and one relaxation 
mode 5 the latter having a relaxation time T.j. On the other 
hand, for ck» only propagating spin waves with linear 
dispersion were obtained. It is now well established fact 
that relaxation of both m and q are non-exponential (slower 
than exponential; either power law or logarithmic) [76,24,141 ,281] 
for all T < All the simple mean field type approaches 
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result of thermal fluctuations exactly as in case of Brownian 
particles suspended in a liquid. 

In a system of uniaxial superparamagne tic particles 
iremanent mgne tization relaxes as 

-t/z 

= Mg e . (1.59) 


where is the full magnetization when the field is switched 
off and the relaxation time z is given by 


1 





-KV/kgl 


(1.60) 


Zq 10~^ sec.) arises due to the random thermal fluctuations, 
K is anisotropy constant and V is the volume of the particle. 

The appearance of V in the exponent is the most crucial factor 
in this relaxation, Bor a particular measurement, characterized 
by a typical measurement time particles having volume 

greater than a critical size appear frozen because their 
relaxation tine will be greater than the measurement time -i^. 

As an example, a spherical iron particle of radius 115 
possessing only first order crystalline anisotropy will have 
a relaxation time of sec. at room temperature. On the other 
hand, a particle of radius 150 A° under the same circumstances 

Q 

will have a relaxation time 10 sec., and hence, extremely 
stable. Usually z = 10^ sec. is taken as a rough criterion 
for stability and this occurs at a barrier height approximately 
equal to 25 k^T, The temperature at which such stabilization 
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(popularly called Hocking) occurs is called the "blocking 
temperature " . 

Wohlfarth [ 28^ , in a paper entitled "SG exhibit rock 
magnetism", drew analogy between the blocking of superpara- 
magnetic clusters and SG "freezing", replacing single domain 
particles of the former theory by non-interacting clusters 
of spins in the latter. The evidence for the formation of 
such clusters at temperatures T»Tg comes from the observed 
deviation of the ^(T) curve from the Ourie-Weiss behaviour 
[73]. Qualitatively speaking, blocked clusters do not respond 
to external magnetic field and hence do not contribute to 
susceptibility. This explains a sharp drop in X(T) at 
T_. Though Wohlfarth could not show explicitly the physical 
origin of the anisotropy energy he gave a list of possible 
sources [286]: 

(i) surface anisotropy 

(ii) magneto- crystalline anisotropy 

(iii) shape anisotropy 

(iv) magnetostrictive strain anisotropy 

(v) other magnetostatic dipole-dipole effects 

(vi) ani so trojic exchange effects. 

Wohlfarth [287] realized that a distribution of blocking 
temperature is closer to physical reality. Therefore, if 
f (T) is the distribution of blocking temperatures, the suscep- 
tibilil^y is given by 

T 

X(T) =1 i f(T') dT’ 


(1.27e) 
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where Q is the Curie constant. Physically, this expression 
(1.27c) stands for the fact that only those clusters will 
contribute to the susceptibility whose blocking temperature 
Tg < T i.e., those clusters which have not yet 'frozen'. No 
first principle derivation of f(T) was known at that tine, 
though it could be determined from the experimental values 
of X(T) using the relation 

f(T) = |^*^^^^/lim ( XT) (1.61) 

T •* 00 

Van Puyneveldt and Mulder [ 288] used this relation to 
calculate f(T) for CuMn (0.7 at %) and Co 0 .Al 205 sSi 02 (14.3 
at fo of Go). 

In order to explain the existing neutron scattering 
experiments Murani [ 289] introduced the concept of distri- 
bution of relaxation times g(T; ) in a similar manner. The 
expression for X( t ) turned out to be similar 

x(t ) = ^ g( T') d -c' (1. 211) 

■I • 

where G ' is a different effective Curie constant. Equation 
(1,27f) expresses the fact that only those clusters whose life 
time is less than the measurement time contribute to suscep- 
tibility. 
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1*7 . 1 Role of Time and Length Scales of Measurements in 
Theories of Third Kind ; 

E-very experiment probes a certain length of a system 
for a certain interval of time, A table of characteristic 
tins and length scales of some common experiments are given 
in Table 1 . 


TAELS 1 


Experiment 

Characteristic time 

C harac te ris ti c 
len; 2 :th 

a, c. suscepti- 
bility 

Mb*ssbauer 

-5 -1 

10 sec, -10 sec, 

-q -8 

10 sec. -10 sec. 

1 

10 "'^cm. 

p SR 

- 11 -5 

10 sec. -10 sec. 

10' ’'^cm. 

Neutron 

scattering 

- 1 P — 1 1 

10 “^sec. -10 sec. 

1A. 

ESR 

-9 -10 

10 ^ -10 ' sec. 



?/e know that blocking occurs when % > Therefore, T^ 

depends on - greater the frequency lower is Tg for given 
V. This explains the frequency dependence of the cusp- 
temperature in the a.c, susceptibility. The frequency 
dependence of T for not only the canonical SG but also of 
insulating SG [290] and transition metal- transition metal 
SG alloys [291] can be explained by this model. Moreover, 
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T (neutron scatt.) > (Mosabauer) > T ( a. c susceptibility ) 

S So 

because (neutron scatt.) < (Mossbauer) < (a.c. 

susceptibility) . On the other hand, in neutron scattering 
experiioants, higher the value of q sinaHer is the size of 
the clusters explored and hence lower is the blocldng temper- 
ature. BesidesjZibold [ 292] and Hardiman [295] suggested the 
idea that ferromagnetic clusters aia responsible for the 
different magnitudes of the frequency shifts in type I and 
type II SO. The consequence of the different time scales of 
Mossbauer spectroscopy and SE have been discussed by 
Hartmann-Boutron [294] • 

Here one must notice that the frequency dependence of 
X (T) and difference between the transition temperatures 

3, * C • 

f rom ^T) and neutron scattering are not only universal 
features of SG but also for dipole glasses. In latter kinds 
of systems, e.g. (KBr)Q^^Q (KcN)q^^q, ^ne measuies the frequency 
dependent dielectric constant instead of susceptibility. Such 
measurements by Eowe et.al. and by Bhattacharya et.al. [ 295 ] 
showed characteristics similar to SG. Hence, these can also 
be explained by theories of thiard kind. 

One additional feature of the model is that it does 
not exhibit any sharp anomaly in the specific heat because 
blocking is a gradual process where the spin degrees of 
freedom start freezing from temperatures much above T^ and 
continue up to very low temperatures. 
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This model also explains the Yo gel- Pule her law [ 296] . 
But one remarkable feature of this derivation was that the 
clusters had to be temperature dependent. A detailed compari- 
son of the dynamical properties of rock magnetism and spin 
glass has been carried out recently by Souletie and by Bormann 
et.al, [ 297 ]. Malozemoff and Imry (MI) [ 298 ] measured d.c. 
susceptibility by using SQUID susceptometer with slowest 
cooling rate possible . The susceptibility peak turned out to 
be independent of the cooling rate apparently supporting 
theories of secend kind. But such claims have been criticized 
by Mydosh [33] on the following grounds - 

(i) Mi's experiment used static field cooling thereby 
'preparing' the sample in an equilibrium state which was 
retained so long as the cooling field was kept constant, 

(ii) The time scales of a.c» susceptibility measurements is 

•“1 3 

10 sec. as compared to 10 sec. in SQUID magaetization 
measurement and the frequency dependence of M and X may be 
significantly different in these two regions. Eecent experi- 
ment of V/ebb et.al. [299] that shows the ferromagnetic resonance 
field depends on the rate of cooling as is expected from 
theories of third kind. 

The main drawbacks of such models are the following: 

(i) the susceptibility peak is rounded, though quite 

sharp, unlit© the sharp cusp observed experimentally. 



97 


(ii) while the superparamagnetic particles in rocks are 
indeed very dilute, they may not be so in SG and 
hence intercluster interaction may not be negligible , 

(iii) if chemical clustering is assumed to be the dominant 
source of cluster formation it will be inconsistent 
with the observed scaling laws, 

(iv) the origin of the anisotropy energy is not clear, 

(v) there is no prescription for calculating fCfg) from 
first principles. 

One must notice that Wohlfarth and his group were mainly 
concerned with susceptibility and specific heat (though of 
course , they derived Vogel-Fulcher law) . An equivalent pheno- 
menological model was suggested by Tholence and fourmier [ 300] 
much earlier. The latter model was applied for explaining 
irreversibility and remanence effects in SG. From a study 
of the remanent magnetization Tholence and Toumier concluded 
that a gaussian form is a good representation of the distri- 
bution of uncompensated moment of "clouds" of spins, where a 
"cloud" containing n spins has a moment proportional to n*. 

A rigorous derivation of the latter statement will be presented 
in this thesis. 

So far as the genesis of the clusters is concerned, 
there have been two distinct views - one group led by P.l- 
Beck views the cluster formation a consequence of compositional 
short range order (mictomagnetism) and hence have temperature- 
independent size? the other group views the cluster formation 
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as a dynamic process and hence clusters possess temperature 
dependent size. The latter category will be discussed in 

the Chapter 6 of this thesis. 

Although Tholence [301] and later other authors have 
shown that T obeys a Vo gel- Fulcher law [8 ] rather than 
Arrhenins law, the controversy is not yet closed. 

1.7*2 Two level Systems (FIiSl and Anderson- Halperin-Varma 

(AHV) ModeT ; 

We have been discussing relaxation of magnetization 
in SG-. But that is not the only^of relaxation studied in such 
systems. The other important relaxation phenomenon is energy 
relaxation, Berton et.al, [ 302] emphasized the difference 
between these two relaxation processes. In their experiment 
the SG system ’'acq.uired some energy by following a "thermo- 
dynamic path" 0(H,t), T(t)) (characteristic of the history 
which the sample has gone to reach the experimental conditions)" 
and the variation of this energy with time at a given fixed 
temperature T in zero field was called the energy relaxation. 

The interpretation of these results and the comparison of 
these two types of relaxation needs the introduction of the 
phenomenological theory of SG as TBS [ 303] first introduced 
by AHV p04] as an attempt to explain the low- temperature 
linear specific heat in all glasses including SG. In this 
TLS picture it becomes obvious that a splitting of the two 
levels is essential for energy relaxation but not so for 
remanent magnetization relaxation. 
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AHV’s idea consisted mainly of the follovd.ng parts- 

(a) in a disordered solid (disordered structurally or 
magnetically) certain atoms or groups of atoms can 
have two available positions (or, 'directions’ of 
spins in SG-) which correspond to two minima of the 
free energy separated by a finite barrier 

(b) due to disorder the splitting of the two minima A is 
a random function v/hose probability distribution can 
be assumed to be a constant. 

The energy can be written as 

00 

n ( A B) f ( A E) d( A E) 

00 

= I n( A E ) 
o 

and hence, the specific heat is given by 

2 5 

C kgT n(0) 

i.e. linear in T. Though this adhoc phenomenological model 
worked well, the physical origin of the two-levels was not 
very clear. The physical origin was later explained by 
Philips [305] for chalcogenide glasses, by Krishnan [306] 
for met glasses, and by Cyrot [198] for SO. Oyrot's idea 
involved two basic ingradients: 

(a) that of Onsagar's reaction field, and 
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(b) that in SG the magnitude of the reaction field 

be greater than that of the cavity field. In any 
SG there will be certain finite number of sites 
where condition (b) will be satisfied. At those 
sites reaction field will be along the Weiss field 
be cause 


Weiss field = Gavity field + Reaction field 


Hence , depending on the direction of the cavity field in any 
Ising SG there can be two energy levels available - one with 
energy - <S^> “ ^i (spin parallel to the cavity field) 

and the other with energy + < Sj^> -X ^ (spin anti- 

parallel to the cavity field but parallel to the Weiss field), 
the first being the stable and the second is metastable state, 


^i = 


13 


X. 

3 


The same argument is not so obvious in case of Heisenberg 
model. However, Gyrot [198] considered two possibilities that 
can introduce uniaxial (or, better, unidirectional) anisotropy - 
inclusion of dipolar coupling and anisotropy of the reaction 
field itself. In such cases, all but one of the available 
states w-1 1 1 be metastable. 

Similarity of low temperature behaviour of specific 
heat and similarity of the explanation in terms of TLS ’^ave 
often been cited as an evidence for SG as a nonequilibrium 
state and hence the terminology glass. The elongation of a 
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rubber as it responds to an exteriially applied stress is 
similar in many respects to the magnetization in a SG. 
Compliance in the former is the analogue of the susceptibility 
in the latter, and both exhibit Curie-like behaTlour much 
above glass transition temperatures. Besides, the elongation 
that remains after withdrawing the stress below in a stress- 
cooled rubber can be called thermoremanent elongation (T.R.E.) 
by analogy with Thermoremanent Magnetization (T.R.M.) Even the 
time dependences of T.R.M, and T.R.E. seem to be identical, at 
least qualitatively. 

Experimental evidence for TLS nature of SG comes not 
only from the low temperature linear specific heat but also 
from the simultaneous study of magnetization and energy 
relaxations in these systems [3071. In order to study the 
energy relaxation experimentally the sample is connected with 
a source of heat and the flow of heat from the latter to the 
former can be regulated and measured. The rate of loss of 
energy (i.e., energy relaxation rate) is equal to the rate of 
gain of energy (which can be measured directly) provided the 
temperature remains unaltered. Eor magnetization relaxation 
no splitting of the two levels is necessary whereas for energy 
relaxation such splitting is essential, McAlister and 
Freeman [ 308] measured the total field cooled magnetization 
M^ along the direction of cooling field. Then rotating the 
sample by 90 ° very quickly again measured the ma^etization. 

The difference between the two values could be explained by 
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yi'ot s work [198] has been extended [309] to understand 
th o igia of the ELS. Apparently, two level picture indicates 
Ising character of the spins. Why does, then, a Heisenberg 
SO exhibit apparent I sing-like behaviour? Ihe potential energy 
of a claosioal vector spin can be written as 

U = - S“ - 1 2°"*^ ^ qP 

IG i 2 ^i ^ ^i 


where is the cavity field at the ith site, 





X being the local susceptibility tensor with a frozen spin 
S^. The local anisotropy arises due to the anisotropy of 
local susceptibility. Therefore, there exists a local easy 
axis (chosen as the polar axis) of magnetization at every site 
"against which the cavity field will flight". If the 
anisotropy fisi(i >> transverse component then we 
shall get the TLS picture-the two-levels beir:^ degenerate 
provided = 0- An alternative proof for TLS for the XT 
model was given by Villian [310], 


1*8 Metastable .qtatAc . 

Once we have understood the physical nature and origin 
of the metastable states, we shall next like to count them. 

One obvious method is to use alternate ’warming* and ‘cooling’ 
algorithms for a computer generated finite lattice with a 


■npr-fii o,nl ar* tir-no f 




T-vTJnrrvr \ 
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and count all the metastable states in which the system 
settles after each 'cooling' process. Such computations 
have been carried by Das Gupta et.al.[ 311] , and both the 
postulates (a) and (b) of the phenomenologioal AHV model were 
confirmed within the limits of computational error. 

Analytical methods are also available for counting 
the number of metastable states. Since the metastability 
condition is defined by 


2l. J. . 

13 


= X . 
3 1 




the average number of local minima (metastable states) is 
given by 



where implies taking trace over all spin configurations. 

-1 

The prefactor ensures exclusion of the trivial degeneracy 
that arises due to the simultaneous flipping of all spins. 
Notice that the factor [ ] av average number of local 

minima for a given configuration and hence the trace and the 
integrations are to be carried out as specified. Tanaka 
and Edwards [ 312] using this method showed < 8 q > "fco of 
the form e°^^ ( a <1 fo^ spin dimensionalities) for the 
Ising and XY models with gaussian bond distribution, N being 
the number of spins in the system. On the other hand, 

EA [ 115 ] showed < gQ> = 2^ in two dimensions. 
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Here again, as in the case of free energy calculation, 
a controversy arises as to the configuration averaging - 
whether 5 ~ ) should be first configu- 

ration-averaged or the configuration averaging should be 
carried out at the end. The reason will become clear very 
soon. Since the average number reveals only a gross feature, 
it is better to study the distribution of the me tas table 
states Hg (E). One can calculate both In [HgCSjg^^nd [In 1^(E)] 
which are the analogues of free energies for annealed and 
quenched disorder. Bray and Moore carried out these computa- 
tions for long- ranged SK model [313], short-ranged (nearest 
neighbour) Ising model [ 314] , vector spins in the SK model 
[ 315 ] , m-component TAP model [ 3 I 6 ], and observed, in all these 


av 


cases, that there exists a critical energy E above which both 

Q 


[In H (E) ] and In [H (E)] 

av - av 

become identical. This result was interpreted as expressing 
absence of correlations between the metastable states for 
E > Ec* "The onset of correlations between metastable states 
which develops at E may be regarded as a phase transition 
with energy taking the role usually played by temperature" 
[314]. Roberts [317] solved the problem in the presence of 
external magnetic field for SK model at T=0. Hertz's dynami- 
cal theory [162] can now be interpreted as followsi- the 
barrier between any two local minima diverges in the infinite 
system and hence any finite-time experiment probes only 
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fluctuations around a particular miniinum. In any finite 
system g (the theory remains stable for t << e”^) depends 
on the barrier height. Even the infinite system is never in 
equilibrium although it takes infinitity long time to hop from 
one metastable s'oate to another. However, Kaplan, on the 
basis of lyons-Kaplaa theorem [318], showed that "in the 
classical Heisenberg linear chain with arbitia,ry nearest 
neighbour interactions and nonzero but otherwise arbitrary 
next nearest neighbour interactions there are no metastable 
states" [319]. Therefore, he commented "that large competition 
plus randomness does not automatically imply, the existence of 
a large number of metastable states". Kaplan [ 319] also 
compared critically the definition of metastable state adopted 
by earlier workers with his definition and pointed out the 
errors committed in these papers. Most of the earliest works 
defined hb testability with respect to sin^e spin flip whereas 
the correct definition would be 'stability with respect to 
deviations (or, flips) of arbitraiy number of spins'. 

Finally Kaplan concluded that these earlier papers "leave 
unanswered the question of the existence of large number of 
low-lying local minima in two and three-dimensional vector- 
spin model SG-". Therefore the validity of Bray and Moore's 
objections against the exponential decay of entropy in the 
Sommers model (see article I. 4 . 2 ) remains debatable. 
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1.9 SO Transition as Looalizatioa Delocalization Problem 

Anderson [320] first made an attempt to apply the 
concepts of metal-insulator transition and mobility edge to 
the SO transition. We know that a magnetic ordering takes 
place in ordered systems when X (T) (q) = 1y 

being the largest eigenvalue of the matrix 1 . If J(q) is 
maximum for q = 0 the ordered phase is ferromagnetic and for 
it is antiferromagnetic. One must notice that in such 
cases the corresponding eigen function is extended throughout 
the macroscopic system. This situation changes drastically 
if the matrix J is random as in case of SO. A random matrix 
has both localized as well as extended eigenstates. If the 
eigenfunction corresponding to the largest eigenvalue of J is 
localized, as is usually the case, one gets a hypothetical 
transition temperature- corresponding to a localized eigen- 
function. It is these eigenstates which so convincin^y 
mimic the "clusrers" of the Ne'el interpretation of SG 
phenomenology above [ S ] • since no phase transition 

can take place in finite systems the hypothetical transition 
mentioned above is meaningless. The SG transition was associ- 
ated with the mobility edge and the random frozen pattern 
corresponds to the first extended eigenstate of the exchange 
matrix [ 521] . This result, howe'ver, has been contradicted 
by Sompolinsky [129] who showed that ’’SG 03rder is not 
associated with macroscopic condensation of modes, but is 
characterized by a distribution of magnetizations throughout 



tho whole spectrum, each of the F eigenstates a.cquires a 
magnetization of 0(1)”. However, very recently Ueno [322] 
showed that there exist an infinite number of condensed 
modes in SK model when analysed in the light of TAP approach. 
Besides, such a treatment led naturally to the Pa-T hypothe- 
sis. Bray and Moore [ 323 ] recast the same problem in terms 
of the eigenstates of the susceptibility matrix ^ rather 
than those of the exchange matrix J. In this case the 
critical temperature is that at which an eigenvalue of 
first vanishes. Ihe density of states vanished at zero eigen- 
value, in contradiction with the result of Hertz et.al. [ 321] . 
The most rigorous quantitative formulation of these ideas 
was proposed by G-inzburg [324]. But here we are not concerned 
with localization- delocalization problem. What is more 
interesting for us is the first indication that probably the 
correlated behaviour of the spins and the cluster picture 
are two different manifestations of the same phenomenon. 
Anderson [ 8] suggested that the situation in SG- is more like 
a percolation theory than a localization theory. Besides 
Hertz et.al, [321] commented "the near divergence of the 
susceptibility associated with the localized states finds a 
natural physical interpretation in terms of spin clusters”. 

1.10 Computer Simulation ; 

Some results of computer simulation have already been 
reviewed in some related contexts earlier in this thesis. 
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The results of computer simulation by Binder and coworkers 

have been extensively reviewed in the references [ 3? 21-24 

and those by TiTalker-Walstedt and by Ghing-Huber and coworkers 

have been reviewed in the References [25-27] • So far as 

question of the nature of the SG transition is concerned 

computer simulations give inconclusive answer. However, 

there is a strong indication that a phase transition may take 

place, if at all, at T=0. Similar results were obtained 

by Reed [ 325 ] and by Banavar and Gieplak [326,32'^ . Reed 

argued that since a SG can be considered as a unit cell and 

since the critical temperature was observed to fall with the 

increase of unit- cell size, T should vanish for a SG. 

e 

Banavar and Gieplak, on the other hand, applied same type of 
arguments applied by Abrahams et.al. [328] in the study of 
Anderson localization. This investigation is based on the 
dependence of the free energy on the boundary conditions. 

This result can be interpreted in the light of percolation 
model because at T=0 radius of the sphere of correlation 
diverges and hence all the spins belong to the infinite 
cluster. Recent results of Kinzel [329] supports the view 
that "SG properties can qualitatively be understood in terms 
of the relaxation of a distribution of independent clusters", 
like those in the theories of INblrd kind. Palmer and 
Pond [330] performed zero temperature computer simulation of 
the SK model in the absence of a field and Bantilan and 
Palmer [331] extended it to include nonzero external field. 
The results of these investigations agreed with Paris! 



replica symmetry- broken theory. 

Ivlost of the computer simulation v/ork done so far 
assumed nearest neighbour gaussian distributed random exchange 
interactions. Only few attempts [332-334] have been made to 
work with realistic RKKf interactions. Besides, -various type 
of algorithms have also been tried [335,336]. Pahnle [337] 
has carried out computer simulation of the distribution of 
magnetic dipole fields - the distribution is approximately 
Gaussian. Computer simulations suffer from the drawback that 
they are performed (i) on finite systems, (ii) for a finite 
time and (iii) usually for classical spins. 

1.11 Some Other Properties 

In this section we shall discuss briefly some lesser 
studied properties of SG. Yery little is understood about 
the transport and galvanomagnetic properties of SG, The low 
temperature behaviour of the resistivity in SG was 

explained by Rivier and Adkins [ 338 ] in terms of the scattering 
of conduction electrons by the diffusive elementary excitations 
having veiy long life time. However, lischer [339], taking 
into account the Dzyalo ;?hinski-Volovik modes [240], showed 
that P (T) = A with positive A and B in dis- 

agreement with experiments as well as Rivier and Adkins [ 338 ], 
The possible reason for disagreement with Ri-vier and Adkins 

3/2 

was pointed out and the observed T ' behaviour was attributed 
to the existence of predominantly ferromagnetic clusters. 
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In the papers cited above resistivity means only two 
contributions to electron scattering - 

(i) scattering due to potential fluctuations, and 

(ii) s-d scattering. On the other hand, Lisyanskii 
and Trush [340] calculated the following contri- 
butions to resistivity: 

(a) electron-electron scattering (relaxation time 

(b) electron nonmagnetic atom scattering (relaxation 
time T;j_) 

(c) electron-magnetic atom scattering (relaxation 
time Tg) 

(d) electron- phonon scattering (relaxation time '^p^) 

in the limiting cases of hydrodynamic relaxation { q / ^ 

and kinetic relaxation ( !)• Campbell et.al. 

[ 341 ] have developed a semi phenomenological model for resis- 
tivity of SG. This model utilized the elementary excitations 
studied by Walker and Walstedt [259]- 

The Kondo-effect in SG have been studied by Larsen 
[342] and lischer [ 343] using different techniques. The high 
temperature In T tail of resistivity and giant thermopower 
[ 344 ] were interpreted by Pischer to be "reminiscent of the 
Kondo effect". 

Theory of magnetoresistance was developed by 
Mooker 3 ee[ 345 ] neglecting the effect of spin dynamics. 



PART B 


HOMOSEHEOUS PEAH BIELD THEORY 
OB 

COOPERATIVE PHASE TRANSITION MODEL 



CHAFER 2* 

ISIHG SPIH gLASS WIFH EKKY 
IHEEAGFIOH 


In this chapter we shall investigate some particular 
featuies of canonical SG, e.g. AuFe and Gulin, which will be 
modelled as a quenched random Ising spin system in three 
dimensions. The spins, occupying the lattice sites randomly, 
interact via indirect Ruderman-Kittel-Kasuya-Yoshida (HKKY) 
interaction, as discussed in Chapter 1, Undamped EKKY 
interaction is effectively infinite- ranged. The interaction 
J(r) is non- random and randomness lies in the occupation of 
the lattice sites. Such infinite range model [143] is more 
realistic than the infinitely weak infinitely long-ranged 
Sheriington-Kirkpatrick (SK) model [ 140] . 

2.1 Local Field Distribution ; lorentzian or Gaussian ?: 

As discussed in section 1.4.1 of Chapter 1, local 
field distiibution (IIL) in Ising SG was derived by Marsimll 
[ 102] , Liu [105], Klein and Brout [103 ] and Adkins and Eivier 
[ 127] to be lorentzian in mean field approximation (MPA), i.e.. 


A part of this chapter has been published in 
J. Phys. PI3 (1983) 365 and the other part has 
been submitted for publication. 
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P(H) = 


+ (H - 


( 2 . 1 ) 


with a mean the width A is proportional to c, the 
concentration of the magnetic constituent. On the other 
hand, Mooker^ee [ 143] and others [145-147] showed that P(H) 
is Gaussian given by equation (1.40) which could be recast 
as 

1 


P(H) 


V It A 


exp[- (H-H^)W] 


(1 ,40a) 


where , 

A = (2c Q 

Our aim in this section is to decide whether the LED is 
truly lovertzian or Gaussian and see under what circumstances 
they agree with each other. 

LPD (2.1) reduces to (1.40a) and vice-versa provided 


(K- H^)' 
a2 


« 1 


( 2 . 2 ) 


In other T/ords, IPD (2.1) and (1.40a) agree with each other 
when condition (2.2) is satisfied. Our aim in this section, 
in other words, is to discuss the physical implications of 
condition (2.2) in some detail. It will be valid in the 
following situations: 

(i) large Q, (ii) large c, (iii) large and (iv) small 
(H-Hq)^. All the conditions (i) - (iii) follow from the 
same condition — that A should be large. 
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Q increases monotonously with decreasing temperature (to be 
shown in the next section) thereby implying indirectly the 
•validity of condition (2.2) at very low temperatures. 

Situation (ii) should be interpreted with caution. The 
concentration c must be much less than the percolation 
threshold for long range order so that even short range order 
can be neglected. Situation (iii) will arise when the 
magnitude of exchange energy is large . This is physically 
obvious as large J would imply large free energy per spin 
as usual for Lorentzian dist2?ibution. Condition (iv) implies 
that the centre of the distribution can be approximated by 
a lorsntzian. Although conditions (i) - (iv) look very 
simple, they have never been stated clearly and systematically 
in the literature. 

Inadequacies of oust two order parameters M and Q 
to describe the SG fphase* have been pointed out by Khurana 
and Hertz [ 116] , by Parisi [II9-I23 and others [122,129] and 
hence the approximate nature of (I.40), 

Finally, we conclude that IFD in Ising SG is neither 
Lcrentzian nor Gaussian but a much more complicated function 
given by the infinite series (1.34). However, in the 
following chapters in this thesis we shall use (I.40) rather 
than (2.1) because the former is much more physical than 
the latter, and poses no problems of diverging moments. 
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2.2 Variation of Order Parameters M and Q with Temperature 
and External Magnetic Field ; 

Our aim in this section is to find out, for gi'van 
canonical SG alloys (characterized by given 1^, Pgff)? 
the correct physical self-consistent solutions of the two 
coupled integral equations (1.41) and (1.42) numerically 
by an iterative process. 

The characteristics of the given alloys were taken 
from various experiments (Table III of Reference [85] ) 
and are shown in Table 2. 


TABLE 2 


Alloy 

t 

\ Concentration of magnetic 
{ constituent (in at.%) 

t 

1 . . 

■'= 

(in »K) 

r 

I T 

! g ^ 

5 (in °K) 

j_ 

1 Poff 

f 

1 

T 

AuEe 

6.6 

2.0 

25.5 

6.0 

CuMn 

0.7 

5.0 

8.0 

5.0 


The best way to calculate T^ is from the intercept on the T 
axis of the high temperature susceptibility X(T). Peff 
was also estimated from susceptibility data. 

The choice of the initial input values for M and Q is 
very crucial for ensuring convergence to the correct physical 
roots, because the transcendental equations (1.41) and (1.42) 
have infinite number of solutions. However, the correct 
self-consistent solutions for T = 0 were obtained as follows: 



116 


Let, 


and, 


then, 


T M + <• = A 


Tg '^Q = cr , 


M = 




71 


00 

I 

-00 


'xV2 


<7'3C + A 

tanh ( — ^ ) dx 


Q = 




2 71 -CD 


00 ^2 

f 


Cfx+ A 


e"^ tanh^ ( — ^ ) dx 


Again., putting + A = y, 


f tanh(y/T) dy 


M = 


V^o 


-00 


Q = 


1 I tanh 2 (y/T) dy 


VTlc a .5^ 


Since lim tanh (+ x) = + 1, 

X -» 00 


00 9 P 

M(T^ 0 ,H) = — - — [ J /^cr 


2 ’ 7 r 0 o 


00 


dy 


- J e-(y+^)^/2c?k 

0 


= I " ) 


T.M + « 


^2 CJ- 


(2.3) 


where Br(x) is the error function, and 


00 


2/^2 


Q(L=0,H) 


J_ J M7-^r/2a^ 




-CD 


dy = 1 


(2.4) 
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Substituting (2.4) into (2.3) we get 

T M + a 

H('r=0,H) =Br (-S_ ) (2.5) 

which can be easily solved self- consistently. In case of 
Auge (6.6 at.%), for example, M(T=0,H=1 KG) 0.01 together 
with Q(T=0,H) = 1. These values for I>;i(T=0,H) and Q(T=0,H) 
were used as Initial inputs for M(T=0.02, H) and Q(T=0,02,H). 
In other words, we move up the T scale, for fixed H, solving 
for M(T) and Q(T) self-consistently using M(T-6 T) and 
Q(T-6 T) as initial inputs. This is based on the assumption 
that M(T) and Q(T) vary smoothly with temperature, i.e 
M(T) and Q(T) differ infinitesimally from M(T-6 T) and 
Q(T-6 T) respectively for infinitesimally small 6 T. 
Therefore, in this procedure jumping from the physical branch 
to unphysical branch of solutions, during iteration, is very 
unlikely. Besides, from analytical solutions of (1.41) and 
(1.42) we know that in the SG phase, for H = 0, M(T) = 0 for 
all T and Q(T) = 0 for I > In-- These facts were used as a 
check of our computer program. This procedure was repeated 
for various values of H(from H=1 KG to H=16 KG), each time 
beginning with the corresponding M(T=0,H) and Q(T=0,H) =1. 

The variation of M and Q with temperature T for two 
fixed values of H (H=0 and H = 16 KG) are shown in Figure 17, 

Since Q is dimensionally ‘similar' to M (both'^'S ) 

2 

we plotted Q against M , But no general trend emerged. 



The temperature variation of the order p 
nod v(T,H) without ext»--n:ia] fif-’.d 
symbol 3 ) ‘-lid vdth 16K'1 applied external 
(full symbols) for the two alloys j^?e ( 
and CiiMa (0.7 at 1j; = 0) is zero 


firm:- 

(op 

i 

4 ~* > . . 
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Surprisingly, AQ = Q(H) - Q(0), turned out to be approxiimtely 

2 

proportional to M over wide; ranges of M and the trend was quite 
general, as shown in iigure 18 for G\Mia (0.7 at.?&). This fact 
will be utilized in Chapter 5» 




GHAPgER 3 

MAGMETOEESISMGE WITHOUT SPIN DIHMIGS : J 


sd 


As shown in equation (1.7)j related to 

"^sd* therefore, in order to calculate JgjQQ- from micros- 
copic considerations, must be known. But in the 

Heisenberg type Hamiltonian for the s-d interaction J , is 

sd 

treated as a coarse-grained parameter. It is much more 
convenient to estimate the latter by comparing experimentally 
measured quantities (e.g, magnetization, susceptibility, 
resistivity, magnetoresistance etc,) with the corresponding 
theoretically calculated values, rather than calculating it 
exactly from detailed microscopic considerations. In this 
chapter we shall estimate in AuPe and CuMh SG by comparing 
experimental data on magnetoresistance [ 85 ] with the corres- 
ponding values calculated using Mookerjee's theory [545] . 

An approximate theory for magnetoresistance of spin 
glasses was first introduced by Bfookerjee [ 543 • This theory 
was formulated using Boltzmann transport equation formalism 
in spatially random scattering potential. This approach 
closely followed that of Beal-Monod and Weiner [ 546] for 
dilute alloys. 


Contents in this Chapter have been published in 
J.Phys. P15 (1985), 565. 
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Let us first clearly systematize the assumptions 

involved in the theory, as follows [144]: 

(i) there is no short range order (SEO) due to chemical 
clustering in these alloys, 

(ii) the temperature- concentration range under consider- 
ation is far from the Kondo regime , so that the 
magnetic atoms have well defined local magnetic 
moments which scatter the conduction electrons, 

(iii) the dominant interaction seen by the conduction 
electrons are ion-core potential fluctuations 

V = - Vg because of alloying and s-d coupling with 

well defined local magnetic moments at the magnetic 
sites, 

(iv) the mean fi^e path of conduction electrons 6» a, 
the lattice constant, so that v/e are far from the 
mobility edge and aH conduction electrons are 
delocalized, 

lasistivity can be calculated using the Boltzmann 
transport equation which is valid in the limit 
k6» > 

(vi) Bom approximation can be used. The latter approxi- 
mation requires ka » 7t, which is not necessarily 
guaranteed by assumptions (iv) and (v) . 
the potential fluctuations and s— d coupling are 
isotropic and local, 
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(viii) the concentration of the magnetic component is 

small, so that there is no mictomagnetic clustering. 
The dominant order parameter of interest is only 
the on-site freezing, so that 


[<S><S>] = Q(T) 6 

(ix) the effect of spin-orbit interaction on magneto- 
resistance may be neglected, 

(x) the spin dynamics is ignored so that 

/♦- i CO t 

H(co) =j[<S.(t) S.(o)>] e dt 

av 

may be replaced by Q(T) 6 (co), i.e., the spin system reaches 
equilibrium instantaneously. This last assumption mil be 
relaxed in Chapter 5 so as to incorporate the effect of 
spin dynamics. 

Mookerjee [345] showed that 

o' SurH 

AP = P(H) -P(0) = - c Rq Jg^[M(H) tanh( -g-y ) 

B 

+ 2iQ(H) - Q(0)}n - — S(S+1) } 

Y 

2 

^M(H) Q(H) tanh ( J (3.1) 


where , 
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n being the number of conduction electrons per unit volume. 
Por Au!B^e and GuMn the ratios of the s-d coupling to the 
potential fluctuations, i.e, 1 J/V 1 are about 0.19 and 0.16 
respectively [ 347] , so that we may neglect second and higher 
order terms in in (3.1)* Thus 

o - SRtjH 

AP = - c Sq tanh(^^) +2{Q(H) - Q(0)}] (3.2) 

M(H) and Q(H) were calculated as reported in Chapter 2. 

was calc’ulated by ITigam and Majumdar [83] assuming that the 

Permi energy E^, and the number of electrons per unit volume 

n are the same as those for the corresponding host metals. 

Hence, AP /c R„J , was calculated using (3.2) and shown in 
o sd 

Pigures 19a and 20a for AuPe (6.6 at ^) and Culfa (0.7 at %) 
alloys respectively. The corresponding experimental values 
for ^p /c Rq) for the same two alloys are shown in Pigures 191^ 
and 20b for comparison. 

let us nov/ discuss the general features of the H and 
T dependence of magne tore si stance in iuPe and G\Mi 3G. (Theory 
and experiment show qualitative agreement with a quadratic 
H dependence for small values of H, going over to H (1 <cc< 2) 
as we go to higher fields. Quantitatively, the agreement 
is also fairly good. (Theoretical curves show the following 
features for the magnetoresistance isotherms: highr- temperature 
magne to resistance is small and the isotherms lie near the 
base line; as the temperature decreases negative magneto- 
resistance increases and the isotherms lie below the high 



ic field (kG) 



(b) fixpeiiiiantal Ap/qH for AuP e aa reported by Migaa and ^ajusaar 
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temperature ones. At still lower temperatures magneto- 
resistance again increases and the isotherms cross those 
for high- temperatures and lie above them. In other words, 
AP as a function of I has a maximum. Sine© at these 
temperatures P (T) itself increases with T but with no 
sharp features, we also expect AP /p to have a maxiiojim. 

This feature was reported by Nigam and Majumdar [83] • 
Theoretically, this maximum is around 20K for Auge and 
between 8 and 6K for Oulvin. Experimentally, this maxinam 
teirperature is much lower: below 9K for AuEe and below 4.5K 
for GuMn. More careful theory and analysis of experimental 
data are required to obtain quantitative agreement. 

Finally, was estimated from 


J 


sd 



AP - 


) 


expt. 





2 ^theo. 
sd 


1/2 

] 


where taken from Nigam and Majumdar^s [83] 

experimv-sntal data. We present the calculated values of 
in Table 3* The variation of with H and T arises due to 
both the approximations involved in the theoiy as well as 
random scatter of the experimental data. Besides, in 
estimating various parameters that enter the theoretical 
calculations, experimental data had to be taken from various 
sources. As a result, one could not estimate and compare 
the relative accuracy or reliability of different data. 
Ideally, either there should exist a comprehensive microscopic 
theoiy which would enable us to calculate all the parameters 
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abinitio, or all experiments should have been performed 
on the same samples, preferably by the same group. In 
Ta.ble 4 the values of predicted previously by various 
authors have been shown together with ours. 
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GHAPTJSR 4* 


SPIN ]3Y¥A1\!IGS BBSISTIYITY ; BTiJ3£QUi.GY OF 
SINGLE RliLAZATION TIMS 


It is Jiow experimentally well established that there 
exists a distribution of relaxation times (LliT) in SG 
[ 85 j 288]. On the other hand, all mean field theories 

eg. [270, 271] of cooperative phase transition from paramagnet- 
to-SG phase, except recent dynamical theories of Sompolinsky 
[, 117 ] and Hertz [162], yield a sin^e relaxation time. Even 
in the Sompolinsly-Hertz dynamical approach to infinitely 
weak infinitely long ranged SK model DRT appears only for 
finite sized sys terns | in the theimodynamic limit all the 
relaxation times diverge. 

In order to establish inadequacies of single relaxation 
time we present some further evidences in this chapter. 

We shall show how single relaxation time is not only 
inadequate to account for the observed relaxation of order 
parameters but also insufficient to explain the variation 
of resistivity with temperature. 

Our brief discussion of relaxational spin dynamics 
will have some new features. 


^Submitted for publication. 
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4.1 Spin Dynamics 

The system under consideration is modelled by a 
collection of n spins randomly distributed over IT lattice 
sites. The spins interact with one another with predomi- 
nantly indirect HKKY exchange interaction. Tirst let us 
consider classical Ising spins with a hamiltonian 

^ - S J( 1 r^ - ^ 1 ) S^ (4.1) 

Defining M(t) =I (t)l and Q(t) =[c^(t)]_ , as before, 

-L ci V 1 Q/'V' 

whe re 

\(t) = 2 S P(...,S.,...?t) 

{S} ^ ^ 

and writing the usual Master equation for the spin config- 
uration P( . . . , . j t) we get 

fl = + [<tanh(Ph3^)>] (4.2) 

and 

( ^ ~ T ~) ft = +[c^j_(t)< tanh( ph^)>]^^ (4.3) 

B 

where , 

h. = ? J. . S . 

1 3 13 3 

Equations (4.2) and (4.3) are of central importance in the 
discussion of the spin dynamics. We carry out the config- 
uration averages in (4.2) and (4.3) iu the Bethe-Pierls-Weiss 
approximation to get 



133 


rJl-> M 

dT 

Jd 


- (1 - M(1;) - j (^)5 l^(t) 

rp 5 v2 

“ -V ^ M(t) Q(t) + 


(4.4) 


and 


2 2 

i = - (1 - ^ (^-5) 


We shall solve equations (4.4) and (4.5) for low concen- 
tration (SG- regime ), alloys at various temperature ranges. 


For such alloys < fn-* 


(a) High Temperature (T fg) Solutions ; 

Ife know that the asymptotic values are M = 0 for 
all T and Q = 1 for T > T„. Hence, if we start with small 

g 

M .q-n d Q at t = 0, they vd.ll continue to be small at all 
times, So, in this regins 



(-^) 

dM 

dt 

= - (1 - 

• ^) M(t) 

(4.6) 

and 

f ) 

dt 

= - (1 - 

3j2 ^ 

- -f) Q(t) + (m^) M^(t) 

(4.7) 


Solving (4.6) and (4.7) with the boundary conditions 


M = M(0) at t = 

0 and Q = Q(0) at t = 0, we get 


M(t) 

= M(0) exp (-tA^) 

(4.8) 

Q(t) 

= [Q(0) - DM^(O)] exp(-2t/Tg;^ 

(4.8) 


+ 11^(0) exp(-2t/'^) 

(4.9) 
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where, 



II 

k^T T 

-^(1 - -^) 

Y] ^ T^ 

(4.10) 


-1 

v _ 
g> 

k_T T^ 

—(1 - -f) 

ri ^2^ 

(4.11) 

and 

-1 

D 

gi2 

1 - -S_ 

2! T„ 

(4.12) 


c 


Such a behaviour of M and Q is characteristic of paramagnetic 
phase. The novel feature of our solutions is the coupled 
relaxation of M(t) and Q(t) as can be seen from eq.uation (4.9). 

As T Tg, Q exhibits critical showing down. However, 

M does not shov/ any critical slowing down near T^. Physically 

o 

this is expected, because, as discussed in Chapter 1, only the 
order parameter susceptibility diverges at T whereas linear 

o 

susceptibility associated with M exhibits only a cusp at 




At temperatures very close to (but not exactly at) Tg? 
both the first and third terms in (4.5) cannot be neglected* 
For convenience, we choose our initial condition as M(0) = 0. 
Then (4.5) reduces to 


( 







-f) Q(t) “ -f Q(t) 
T T 


Completing the squares on the right hand side we get, 


^kgT^ dt 
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where , 


and 


Q = Q + b 


b = .i (^5 - 1) 

T_ 


Q}ho2?efore , 


Q(t) 


- 2k T X 

2b Q(0) exp -[- 1 ^ (1 - ^)t] 


2kg T T 


Q(0)[1 - exp (1 p)t}] + 2b 


which is still exponential for -very large as well as -very 
sxoall t. At 1 = 1 

s 

dt M (t) - (^) Q^(t) + .... 

-1 

For large values of t, Q(t)'^ t , since the first term is 
exponential and dominates only at short times. 


(b) Low Temperature (T <T_) Solutions: 

In this regime Q(oo ) 0 although M(oo ) is still zero. 

Let, Q(t) = Q(co ) + q(t) and M(t) = M(oo ) + m(t), where m(t) 
and (l(t) are small amplitude fluctuations around the equili- 
brium values M(oo) and Q(oo )5 at time t. 



136 


Then from (4.4) we get 

2 

^ft^ - T^^ V *2(o° ) (4.13) 

Jj rj]^ 

On the other hand, from (4.5) we get 

T^ 

Q(oo ) =0, or Q(oo ) = ( 1 ^) (4.14) 

g 

It is easy no see that the latter solution is the only stable 

solution for T < T^ while the former is the stable solution 

g 

for T>Tg. Hence substituting (4. 14) into (4.13) and solving 
it 

: k^T 2T T T^ 

M(t) = M(0) exp[- (-^) (1 - + -J^)t] (4.15) 

Similarly, from (4.5) we get 

= - (1 - ) 

? p 

T ijt 2T^ 

- (1 - -§) ^(■fc) + m^(t) - Q(co ) iCt) 

(4.16) 

Substituting (4*14) into (4.16) and solving 

m^{0) 

Q(t) = Q(oo ) +[Q(0) - Q(oo ) - 2 ] 

(2 + -^ - | 2 ) 

TT T 

exp (-2t/T^ 
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m2(o) 


+ 


(2 + 


T 


B 


T'l. 


_S_ 


_ 2T> 


2k^’S 

exp[- (_^_)(T 


2T T 

— # + — ^)t] 

(4.17) 


where , 


— (-1 - 1 ) 
T ] rpS 


It is easy to verify that (1 - 


2T 

c 

T 


4* 


T T T 

-^) = Cl -1^)2 


_ fj- 

+ ^ (~2 1) > 0 for I < I i.e., for both T <T <I and 

T-.''TcvIg, and hence M(t) in equation (4.15) and the last tera 
in (4.17) always decay to zero. Thus, even below T both M 
and U decay exponentially to their equilibrium values in the 
homogeneous mean field theory. There is a unique single 
relaxation time 
earlier. 

Veiy close to, but below, T , we have 


T for T < T v;hich is not true as stated 

O ‘ ' & 


( 


-Jl—N ^ 
2k2T^ dt 


T 


(-f - 1) l(t) - 


T 


T 


' <i3t) 


and hence 


Q(t) = Q(co ) + 


2kT.T tJ 

2b’{Q(0)-Q(oo )}exp[ ^ 

2k T T^ 

Q(0) -Q(oo)}{1-exp[- 2b» 


'n rp‘ 


where , 


b’ 


:i - 

it 


Therefore, for very large as well as very small t, relaxation 
of Q is again exponential. 
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Both real experiments [76 ] and computer simulation C 2811 
revealed that the relaxation is governed by either power law 
of the form 


M(t) = M( 0 )t"^ (4.18) 

or by a logarithmic law 

M(t) = M(0) - S(T)ln t (4.19) 


However, it is difficult to distinguish pov/er law (4.18) vath 
a« 1 from logarithmic law (4.19). We shall derive (4.19) 
using a DRT in Chapter 8 . 

The study of the ^If- correlation function R(t) = 
[<Sj^(t) Sj^(0)>]g_^ involves setting up a Master equation for 
the conditional probability P( 5 ... ?t j . . .,S^, . , , 0 ) , 
Proceeding exactly as before we get 


^_IL) Mill 

^k-nT^ dt 


T^ T 

- (1 - -f ) R(t) + ( 7 ^) M(t) M( 0 ) 

T ^ 

T^ 

- (-f) Q(t) R(t) 


( 4 . 20 ) 


The structure of equation (4.20) is exactly similar to (4.5) 
except for the factor (p/ Zk^T) being replaced by (p/kgT). 
It is therefore easy to see that the equilibrium values for 
R and Q are the same. This is the justification for using 
Q(oo) as the static spin ^ass order parameter. Hence, 
assuming M(0) =0, the solutions for R(t) are 
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R(t) = E(0) exp (-tA ) for T > T 

g 

and 

R(t) = S(co ) + [ R(0) - E.(oo )]exp(-t/T J for T < f . 

^ S 

faking" "bne Fourier ■feransform and remembering "thaf R(oo ) = Q 
we ge 1 5 

E(o) = Q6(o) +[E(0) - Q]-jA— (4.21) 

“ tg 

vYhere , 

r _ r for T > T 

= r for T < T . 

o'- g 

notice that we have not used eigenstate space technique 
used by earlier authors. Besides, our treatment yields 
coupled relaxation of M and Q. 

4.2 Resistivity ; 

In this section we shall apply equations (4.21) to 
calculate electrical resistivity in canonical spin glasses. 

We shall consider two separate sources which contribute to 
conduction electron scattering and hence to resistivity, 
s-electons scattering against the potential fluctuations in 
the alloy and by localised impurity moments on the magnetic 
atoms. The hamiltonian of such a system may be written 


as 
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H = + (1/N) 2 2 S E n Y + 

° k k' cr i V' ® ®kV 


-(1/N) n^ ^ 

•i(k-k’).rj_ 




®k4 \f ^ ® 


( 4 . 22 ) 


Hq is the Hamiltonian of the s-electrons in the host, 
the potential difference between the host and the impurity, 
'^kk’ coupling between the conduction s-electrons 

and the moment carrying d-electrons of the magnetic impurities 
and n^ is the occupation probability of the site r^^ by a 
magnetic impurity. We shall consider both and Jjj. ^ t to 

be isotropic V and J. 

The transition probability X"''(ko' ••k* o * ) is given in 
the Bom approximation [339 ] 


r(ko- 


k’ o') = i(kcr) [1 -f (k’ 




(- 

2 


J - 


00 

J au 6{“-e^„ + e 

""OO 




I / 


00 

r dt e 


imt 


'(J 


-00 


[J^< s. (t).s.(0)> + +0 JY <sS] (4.23) 

where f is the Bermi function and o is the unit vector along 
^ , We shall calculate the spin flip and non spin flip 

scatterings separately. 
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(a) Non-spirx ELip Scattering 

In this case ^ = o'' • Therefore, neglecting any short 

ranged order, the probability of transition is gi'ven by; 


r(kcr-k'cr) = A + B + G 


(4.24) 


whe re , 


K 


i 


. » - f {^,)i 




a 

°° ^ nr rr 

) J 4t<S^(t) SjO)>e 

uu “"(X) 


271 


icit 


J 


r, ® -i 1 

H 1 “i i -f ("],.)! 

-00 ktr — 




whe re 


00 


^zz “ 2n -00 


X dt <sf(t) sf(0)> 


Using the relation f(£jj-) 3 -f(sjj--“)} = 


X'- ' 1 

f(^) - f(^ 


- «) 


1 - e 


pco 


We get, 


A — 


£ 


B n . 
1 1 


00 

/ 


-00 


da 


f(®k) ■ ^ 


_ Ci>) 


1 - e 


Lx 

O 


E 


(4.25) 
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2 OD 


00 




■.<y 

^cr -CO 


ll "■’ 

1 k ' 0- 


£(^) i " £v _) 


k' 1 i ^"-k tr 

Using the relation f(£^) i 1 - f(^T_)I 




k 


B 


b“' 5 (^ ] 1_ 2 a 1 S we 

P ^d£^- H i i F k^-"^ k>0 £0- 


k..-^ 


ll 

1 


r - Q-^ /' df 1 S n. S U(e ) -] 
L P 1 1 1 rr ^ k 1 J 


(4.26) 


Similarly, 


G - (^^''^)L- P ^ (4-27) 


In equations (4- 25) - (4-27) N(\^) is the density of. states 

of electrons with energy %. and spin cr , Since M =[< sS-l ,, = 0 
in the spin glass phase G = 0. Bow carrying out the config- 
uration aye raging, the total non- spin flip scattering contri- 
bution to the transition probability becomes 

,6 L S r{kcr-k' a)]av 
fe'o - ■ 




143 


f{C . u) 


Rzz (“) 


H (^, 
+ ^ k 


-w) {■ 


f(V 




■}] 


( 4 . 28 ) 


where c = 


1 


^ n. is the concentration of mametic 


1 i i 


impurities , 

N^=H( s ]^f) + 


(b) Spin Plip Scattering ; 

In this case a / 0 ’ , and 




p(k0-» k' 0 ') - { 1- f 2 


s 

i “-i 


00 4 T 

1 tH^(“) +Eyy(o)} 


The refold, 

p [ S , rtea>4‘a')]„ 


/df N r o P 


J B 


C 


f ( £ -m) -f ( £ ) 

J " V'“>> ^ 


00 

J 

-00 


(4.29) 


In the spin glass phase because 

of average isotropy of the spin orientations. Now, using 
the relation that 

- (df/d T J = P [ » P o o ) lav 

'.< 4 > (y tmmi 0 " 

and resistivity p =(“^) we have 

ne 



calculated, using single relaxation tine, shovrs cusp at T 
contrary to the smooth heha-viour observed experimentally, 
Ihis is shovfn to be a consequence of the inadequacy of 
cooperative phase transition model with single relaxation 


tin© • 

Chapter 5 is an extension of Chapter 4 to vector spin 
glasses in the presence of external magnetic field. In this 
Chapter we improve Mooker;iee's theory of magnetoresistance 
incorporating effects of spin dynamics. 

Part G of this thesis begins with Chapter 6 where v/e 
develop the percolation model of paramagne t-to-spin glass 
transition and calculate the transition temperature as a 
function of concentration c taking the effects of finite mean 
free path and frustration properly into account. In this 
Chapter ?/e derive Vegel-Pulcher law and explain the variation 
of y/ith time and length scales of probes. Besides, we 

to 

calculate the low- field low-frequency a.c. susceptibility and 
'sho?/ tho.t linear response theory holds. 

In Chapter 7 the mathematical relation betv/een the 
percolation model and Wohlfa,rth's model is established and 
interpreted physically. Comparing our theoretical results 
vj'ith experimental data we calculate the distribution of 


blocking temperatures. 

I The logarithmic relaxation of magnetization is derived 
analytically in Chapter 8 so as to reveal a deeper connection 


with l/f noise in various other physical ' systems. 
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Part B highlights some further inadequacies of the 
coope^ra bive phase transition model of iidY^ards™ Anderson type 
vvith single relaxation time. Part C proposes an alter^ 
natiTe theoretical framework whose predictions are 5 at least 
qualitatively^ closer to the experimental results. The 
latter paves the way for an unified approach to the theories 
of spin glasses by incorporating cooperative effects on the 
one hand, and features of a gradual freezing process on the 
other. 


Part B begins with Chapter 2 where we investigate 
the physical conditions under which both lorentzian and 
Gaussian local field distributions are good approximations 
for I sing spin glasses with HKKY exchange interaction. In 
this chapter we also solve Sherrington-Kirkpatrick equations 
for two canonical spin glasses numerically and self-consis- 
tently so as to find the variation of the order parameters 
v/ith ' temperature and external magnetic field. 

r.iooker jee ' s theory of magnetoresistance has been 
applied to tv/o canonical spin glasses in Chapter 3* 

Comparing the magnetoresistance thus calculated with 
the corresponding experixental data we calculate the strength 
of the s-d exchange interaction, Jg(j* 

The results of our investigatijen of relaxation spin 
dynamics in spin glasses using Master 


been applied in Chapter 4 to calculate 


equation approach has 
the resistivity by 


Boltzmann transport equation method. The resistivity thus 
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Because of the exchange interaction among the moment 
bearing magnetic atoms (or ions) in random magnetic alloys, 
various kinds of magnetic ordering take place depending on 
the tempox'atui-e and the concentration, c, of the magnetic 
constituent ( s) . In this thesis we have investigated the 
nature of the transition from paramagnet-to-spin glass state 
and the behaviour of various physical properties across the 
transition. 

Recently there have been several attempts to regard 
a spin glass as the prototype for many more complicated and 
less understood physical systems. Unfortunately, even the 
theoi-Gtical understanding of the physics of spin glasses 
is far from complete. There have been two distinct and 
apparently irreconcilable approaches - the cooperative 
phase transition model and the model of gradual blocking 
of superparamagnetic clusters. Both of these have had 
limited success. All the successes and limitations of ^ 
these models have been critically reviewed in part A of 

this thesis. 
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-1 


e 


[V^ 






3 P 


00 


-00 


d^R(o ) i^( -a)C 


f(e^ ~") - f(£T^) 




0) 


-1 

(4.30) 


1 ] 


Yla have considered because due to the derimtive of the 

Permi function in expression (4.30), it is clear that only the 
electrons near the Permi surface contribute dominantly to 
resistivity. Assuming, further to to be small 


f(Sy -^0 “ f(ep) ^^pw/4 

which is true for reasonable values of (note that 11=1 
throughout) , 


-1 

e 


= A c^p) 

P 00 

+ 3pJ s 

-00 


^ R(c^) 
(eP“ -1) 


do] 


(4.31) 


ito equation exactly similar to (4.31) was derived by 
Campbell [348] from intuitive-qualitative arguments. Since, 

Ii(oo ) = Q(T), and R(0) = (4.21) 


can be written as 


R(m) = Q(T) 6 (w) + [ - Q(T) ] -g 

6> +Y 

where 

■i —I 

_ “ ‘ or T according as 1 > 1 „ 

Y ^ g> < S 

or T < 1 . low onwards we shall use the symbol % for 

e o 


Then 
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00 

/ 

-00 


W R(co) 
"m ^ ' ' 
e - 1 


a. CO = +l'Lv-S+i 


i2 




IX 


s 


- <50 


05 dW 


(e -1 ) (o) +v ) 


Lot us f irs 0 chock the asymptotic behaviour of eq.uation 
(4»51)* AtT = 0, Q = g /3 and T; = Oi»s. y oo , Besides 
as L -*■ 0, p -* 00 . Hence in this limit 


II 

7T 



G 



dJJ 




dw 


(1 


“f- 


pw 

2\ 


+ 4. 


= 0 


3 '. 


+ , , ) 


and, we have 


p(0) = (-^) c [V^ + J^S^] (4.22) 

ne 

which is the correct loshida limit [ 349] • On the other hand, 
as 1 -* oo , = 0 and V ->-0, In this case, pw is small and 

we have 

PY ^ to d C' c-JL ^ <ito ^ ^ 

n -00 (epu _i)(|+y2) ■■ « -CO 

Hence , 

p(oo) - (-^) c [V^ + S(S+1)] (4.23) 

ne 

which is, again, the correct To shida limit [349]. 

It is interesting to notice that the integral above is 
exactly the same as encountered in computing dynamical suscep- 
tibility in the microscopic theory of spin fluctuations 
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(see integral ( 3 . 8 ) of Ha^n.hnan [ 35c] ) . ine mtegiai 
my be split into two parts as 


00 

_L. r 

7 [ 

-CO 


d 




2 ^ 


% i 


0 

r 


6 Cl) d f ,) 


~oo (eP''1il)(w2+Y^) 


00 

+ ; 

0 


J 6 ) do 




How, 


o 


P to d to 


CX) (d^“-1 ) (0)2+,^) 


00 

/ 


joto d to 


(1 


p- pw., 2 2 . 

® ) (d> + Y ) 


00 
= / ' 


pto eP^ dto 


0 - 1 )(co^+y^) 


ft: 


Y 


CO 

/ 

0 


(gtoeP"^ -Yo I dM 

^ r~r 

(e -1)(o+,/ ) 


00 


+ / 


Pto 


dCO 


(eP“ - 1 )(o% 2 ^ 


CO r> 

^ ___ptodco 


^ 2 2 , 

( to + Y ) 


00 


+ f 

o 


pto d CO 


(eP“ - 1 )(co 24 ^. 2 ) 


1.3 . 


X. 

% 


00 

/ 

-'CO 


do) 




-X 


CD 


Po) 

(to^ +Y 


iL. 


CO 


0 




pto dto 


•i)(to^ +r ) 


4 " 
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The upper limit of Integration in the first integral is teJoen 
bo be ' 0 (kgig) , following Eamakrishnan [ 350J , we absoi'b 

the effect of this term in • As a result H gets renormalized 
to . Hence j finally, resistiyity is gLyen by 


P(T) 



cH^( 




3Pj2[3m +(311+11 


•Q }' 


0(1? 


^ J 


where , 


(4.24) 


0(T,y) = In y - (2y)"’' - >7 (y) 

with y = ,-3/-7rT and 'S (j) is the digamma function. 
Eamakrishnan [ 350 ] has giyen an approximate form for 0{y) 
yalid oyer a large range of yi 


0(y) =-- (2y + 12 j^)~^ 


(4.25) 


let us now inyestigate the behaviour of P(l) at some 

particular temperatures, yiz. 1 = 0, o and 1 = 1. Eor 

S 

breyity let us call 


= i(T) 

ftT Tik^l y+3 ^ ‘ 

(i) Very Low lemperature (1 0) ; 

ils 1 -*■ 0, ly ■* 0 and hence X(l) "* 0. Besides, 

Q(1 = 0) = sV?. 

Therefore , 

p (0) = (^) oS {Ej,)[t2 + S^3 

ne 

which agrees with the Yosbida limit (4^22) . This agreement 

f-i?4tPhl i cshA cs 4!^*^ ^^r£>l n A n /^-f TYiQ 4o -Prv'y'm fA OC;^ 



148 


for T << 


(ii) '^'^ry Hig:fa Temperature r (''P 
In this regime k^Tx = - 


IS T -*■ 00 . 


Therefore X(T) 


and hence k-nl t ^ 
T‘^ 


as T -» 


■Thus . 


P(oo ) C-^) cl (e ) [ ^ J S(>yi . l - 

ne ^ ^ +(3 ai ) - 


which does not agree with the Yoshida limt (4.23), and hence 
the inappld.cability of the approximate form (4.25) for T» T 

g' 

(iii) As T -♦ Tg from both abowe and below, I(T) i go that 


(Tg) = (^) cN^(ej,) [ V2 + j2g(s+i)] 


= P (co ). 

lot us now examine the slope of P (T) near T . 


imow that 


2k-n X 


■B ^ >^ 


T> for T > T 


(■^ + 


2kp T< 


) T< for T < T . 
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Since 


1 

T 

d TT- 


« 






for f near T , 

xy ^ 


2 

'B '*'> 


2kT, %. 


dT 


for T > 1 


and 


d r 


dT 


B 

B 


for T < T . 


Thus the derivative of is negative above T^. and positive 

o 

below Tj^. This is physically consistent because increases 
from zero to infinity as the temperature is decreased from 
large -values to T , while x decreases from infinity to zero 

O 

as the temperatures decreases from' Tg to zero. Notice that 
is the lifetime of de-viations from the frozen equilibrium 
configuration of spins below Tg. 

Now, 

^ [ X - Q} + Q ] 

dT 0 

, C M|±li - a } § + (1 - X) ® 


5 


dT 


=[ 


a ^B ^ 


2 2 2 


+ (1 


dT 




+ 

(■jck T X + 3) (mk^Ti; + 3) 


fdjE) 
• ^dT ^ 


2 , 2 „2 

a ^B ^ (^)IC 

^ MT 


('n:k2TT+ 3)' 


S(S+1) 

3 


- Q} 
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71 k t; 

r -- 5 -- 


(itkgli: + 3 ) 


2 2 2 
TC T 


+ 


Tlk-nT T 

^ (I_ 

( TTk^T t: + 3) ^ 

2 v2 m2 2 

Ti; kjg T T 


(Txkgl^i; + 3) (Tik^TT 4- 3) 


(f- 


Tl 


2kgT; 

__ 


■)] 


C^^-Q 1 + ( 1 -X) §, 


= [ 


2 7i:kT2l 

( 71 kgT X + 3 ) 


- 


tg t. 


T1 


■) - 


o 2 , ™2 2 

2 Ti kgi X 

+ 3) 


■ (- 
2^T 


V - 2 

T) 


f sl|±li_ Q5 . (1 -X)f 2 


2Ti;k2TT: 


(TUCglT^ + 3) 


(4 


n 


)[1 - 


7T;kgT X 

(nk^T r+ 3) 


{Si|+il _ QJ + (-1 


X) 


dT * 


^ (• 


I 

T 


^JL) [ s(s+i) 

71 3 


- Q] 


6 TikgT X 

+ 3) 

+ (1 - X) II . 

2S(S+1 ) m ^ m from abo^ve as X (T ) = 1 

Tl TgTi: g S 

d rye S(S+.11 n } + o 1 = 

7j L xt - y 1+ r]Tr^ 


Similarly, ^ 


for T •* from below. 

s 
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Therefore ^ exhibits a discontinuity at T = T^. in other 
words (4.24) exhibits a cusp at T^. Experimentally no such 
cusp has been observed so far in any of the canonical SG 
[80,341] • s-'t'tiibute this discrepancy to the oversimplified 
picture of hanogenoous MET of cooperative phase transition model 
is ^ single relaxation time* We shall shovir in 
Ohaptex* 9 how a distribution of relaxation times giYes rise 
to better agreement with experimental results. 



GHAPIER 5 * 


SPCT. DYUAMIOg II {IHE PHBSEjICE OF EX!ISBIir.AL 
MAGMEPIG PIEOQP AED MGIilETOIiESISPAlGB 


In this Chapter we shall study the spin dynamics in 
the presence of external magnetic field along the same line 
as adopted in Chapter 4. The results of this investigation 
will be utilized to calculate the magnetoresistance of 
canonical SG. 


5.1 Spin Dynamics in the Presence of Pield 

Generalizing equations (4.4) and (4.5) to -vector SG 
in the presence of external uniform magnetic field H, we get 


dE (t) T„ T 

— ^Tt = ” (1 “ “f) (0) 

B T 


(-f) ""5)^ s^(0y>] 


ar 


( 5 . 1 ) 


where (t) = [<^'^(t)^(G) ] 


av 



* Yarious parts of this Chapter have been 

i) accepted for publication in PhysJiCtt.A (in press) 

ii) submitted for publication. 
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because the precessional motion of a spin S in an external 
magnetic field H is given ^7 ff =y(S x H) , where y is 
the gyromagnetic ratio. 

Next we solve equations (5.1) and (5.2) in the low 

concentration regirae (T < 1 ) at high tempeiatures (T>'I >1„) 

^ a go 

and low temperatures (!r„>T>T^ and T >T >1). In the 

O ^ go 

high temperature regime we can drop the third term on the 
right hand side of equation (5.1) and socond and thiid terms 
on the right hand side of equation (5.2) Let us assume that 
the external field is applied along the z-direction. Then 
the solutions for magnetization are given by 

t(0)i 

■) V = r 

i ( 

( -iq(o) ) 

and, 

M^(t) =1^(00) +lpl^{0) -M^(oo)] exp-(-tAo) (5.4) 

wher-e M 2 (oo ) is the equilibrium value of Notice that 

M^(oo ) and B/y(oo ) both vanish. The physical interpretations 
of equations (5.3) and (5.4) are very clear - the precessional 
motion of the magnetization vector is described by the sine 
and cosine factors in these expressions whereas the thermal 
decay is expressed by the exponential factor involving only 
one characteristic temperature Similarly, solving (5.2) 


Cos(7'Ht) exp(-t/'r ) 

G 


Sin(rHt) exp(-tA^) (5.3) 
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we get. 


• > 




(' Cos(y Ht) exp(-t/T,) 




+ 


( V'°>( 


Sin(YHt) exp(~t/ t ) 


-M^(t)/ 


+ mUo)-< >[i 
/M(t)r 

^ j j 


■3xp(_ ] (5.5) 


where , 


-1 


= (1 - 


T T, 




+ 


E. 


V'°> 

-E^y(O) 


Gos(y Ht) exp-(tA ) 


Sin(YHt) exp-(tA ) 


+ DM (0) 
«/ 



k T 

r 1 / B Gv , T 

[ 1 - exp (- 7 ^)t ] 




(5.6) 


and, 


exp(-tAg) 




+ DM^(O) M^(t) [ 1 - exp (- ] 


(5.7) 
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In the low teii^jerature regime we write S^^j^t) = 
H^j,(co) + rp, 3 ,(t) and M^(t) = M^(co ) + m^(t) and gat, as in 
Chapter 4, 



dt 



I 

1) + (^)mj^(t)m2/0) 




'\x Jar 

(5.8) 


and, 


dm,,(t) 


/ h . — n 
dt 




+ (^) ( m (t) X 2)^ 


(5.9) 


Equations (5*8) and (5*9) can be solved exactly as equations 
(5.1) and (5.2) So the relaxation of E turns out to be exponen- 
tial at all but one temperature. At I = T , Ep, 2 Xt)/^ t ^for 

o 

all \}j V Sit not too short time scales, exactly as in Chapter 4* 
Unfortunately, the experimentally observed decay of 
M and 2 £>- 3 ?e much slower than exponential (probably lo^mthmic) 
and do not agree with the solutions obtained above. The 
correct behaviour will be derived in Chapter 8. 

Let us now calculate 2(u). Proceeding exactly as in 
Chapter (4) and using the forms of 2('^) derived in this 
Chapter we get, 
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{ S(S+ll _ _ 

2 TTT 1- 


+ 


1 


(-^f- + (c^+yEf 


+ 


t^-Q{T)3^ 1 

!_ 


2 n'V 


(:^)^ yH)' 


( 5 . 10 ) 




^ S(S+1 ) 
Q(T)6(w) + — ^ 


- Q(T)} 


%% 


[ 


1 1 

, ^ .2 2 J 

( ) 


(5.11) 


E. 


^(“) = - V = ^-iWc 


(;r) "+ (^ “tHT (^) H) ^ 


(5.12) 

Notice that (5.10) -- (5.12) reduce to the correct limits for 
H = 0 deri'ved in Chapter 4. 


5 .2 Magne to re si stance 

We proceed exactly as in Chapter 4. Eor non spin 

f 

flip scattering { a = a ) 


s r'(^ rr-^k’cr) = A + B + C 

f »— •>«« 


(5.13) 


where , 


A = (J A) 


X 






f(ek) - 

^71 


IN 




oy) 


(5.14) 
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B = (vV») P"'' (|^) I “i 

0 = (2 J T/H){ - <S?>ffS(e,^ )} 

ic — . CT 

Carrying out the configuration averaging 

[ S rCkcj -k'a )] c ) 

k' ay ^ k 


( 5 . 15 ) 

( 5 . 16 ) 


+ 2 5j V M (e ) - J% /R^^(o) 

k -00 


) f 


-M) - f 


60) 

- 1 


} d.o> 


( 5 . 17 ) 


Next we consider spin-flip- scattering. We define, 


and, 


whe re , 



ISTotice that R ^_ is associated with down up spin flip 
scattering, whereas is with up -► down. Hence, taking 
the conservation of energy into account, we get 
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r( kcr - k’ cT ' ) = f (kcr) H - f (k'cr ' ) 3 (J | n. 


1 1 


/ d6^[R^_ (o)) 6 (o)- + e^,^, -y H) 


+ (“) S (w- + Sj^-i f + Y^)] 


Hence , 


[ E r(kcy - k'<y')f 
k’ ' t 


2 H ^krf^ 


CO 'f(S. - +YH)-f(€; ) 


/ dw{- 


f(e, -w-^/H)-f(£, ) 


''jr=wTZT ^ “TT^^WTT 


(5.18) 


How onwards we shall consider ejj- ^ because only the I’ermi 


electrons contribute most dominantly in resistivity, 


H_^(u) = E_(u) + E,„(u) + ifIW(“) - H„(“)! 


E-+(“) = Sxx^“) + Syy'“^ ■ 


Therefore , 


R_^_(63) 


2Q(T)6(«) + 


3 nr 


{ -T 2 2 ^ 

(t) + (co-yH) 


C 1 

nr ,1,2 . . 


(r) +(“-YH)" ii) + (co+yH)' 


(5.19) 


and , 
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s_^(“) = 2Q(T)6(“) + i ' f ' 


(t*-) + (m+tH)' 


_ _SL ( i 

nr . 1 » 2 


+ 


1 


nr ' . 1 s 2 , TT s 2 / 1 v 2 , 

(^) + (co+yH) (-) + (co-yH) 


j-}(5,20) 


Nows 


00 ■ f(S ■■“+YH) - f(%) 

I, = _/^doE^(c,)S H) ■ I-;- i 


f {e„ + YH) - £ (6 ) 

= 2Q(T)f pr-g-g ; 

(e ■] ) 


+ 


2S 


(S+1) ^ 1 

loo ^ ^ 


Q 

TIT 




/ du[ 


-CX) 


] ^ 

/ » N ^ // ■ 


.p - 1 ' (o>-yh) 


+ 


1 


(4)^ + (c-'+rH)^ 


•} 


(5.21) 

S'irst three integrals in (5.21) can he carried out analytically, 
as was done in Chapter 4, assuming 'x; + yH to he small, The 
fourth integral can he written as 


00 

/ dw 

-CO 


r f(^ -“+TH) -f(=y)T, 

I- .l^o-rTIJ : ^'7T 


1 


p(co-rH) __ ^ =■" ^1^2 ^ (03 -yH)^ 


} 


In order to do this integral analytically, we approximate 
it hy 
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CO 


f(Ej, - UHYH)-f( 

^ ^ P(u-yH) ' 

“00 e'^'- ' ^ “ 1 


.]C-T^ — 

(:^) + (0“YH) 


where , 


[ (-h^^ + “o ^ 


where we have replaced witihin the bracket by a constant 
independent of w . In this approximation, 

si-, - — +1 3 

- (-jlf) * ) ! (5.22) 

where j2^(T,t) has the sane approximate form (4.25)* Similarly, 


00 


I2 = / dw R (co) { 

“00 


f (Sy “«“YH) 'f(S) 


P(.<^+yH) 




Therefore, 

* R = p- ^ ^ 3 

(5.23) 


where , 

f II 

Rotice that in the limit H = 0 ,t = x ==t. 

Using the approximate form (4.25) for 0(11, x) we get 
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' ^ v3 m3. '2 '' ”2 , 2 2 2, '2 "2s 

( . .. » , oo ^ ^ ■^'1^ )+^-rc t T +T ) 


0(T,t: ) + 0(T,'v ) .-. 


JL 


o 2 2 • " / f " s 

n ^ T r + 37rk^T( x -n )+ 9 


^B- 


(5.24) 


In order to compare theory with experiment we make a further 
approximation ; 

4 Y H T^«'l 

so that retaining only terms upto the order of H, 

t;' = (1- SVHWj, t:^) 

and, 

" = (1 + 2yH«„ t:^) 

and hence, 

0(T , t;') + 0(T,Tr") = 2iZ)(T,T) 

Hence under this approximation, resistivity is given by [ 349 ] 


p(H>'!:wy4w + ’ ’ 


Wq - 


where , 


'.To = 0 Y^H (ej,) + 3opHH(ej|,)[ , y) 

- ^ rtl,T) ] 


and, 


Therefore, 


a = 2cJVN(ep) 


q2 2 

P(H) - [ W„ - (f-)M^ ] 


(5.25) 
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As T •* 0, T -* 0 and hence Y ■* 1 , and as T -» T ,t ■* oo , 
hence Y 0, Por very large T, t; ■* 0, and again Y 1, 
Besides A Q -* 0 for both T 0 and 1 -» oo . Thus 

n 

(4 M - 3 AQY) > 0 for all T, and since terms of the order 

—K are very snail, AP is negative at all T, in agree nent 
V 

with experiirent and its value determined dominantly by 

(4 - 3 AQ). 

Using the self-consistent solutions for M(H) and 

2 2 

Q(H) obtained in Chapter 2 we plot (4 M — 3 ^Q) against M 

for fixed temperatures in figure 21, Thus, clearly AP 

2 

is proportional to M at least up to 12 KG. This in excellent 
agreement v\rith experimental results [385, 386]. Besides, 
plotting (4M^ - 3 '^Q) against T for fixed H we obtained 
qualitative agreement with Nigam and Majumdar's [83] experimen- 
tal data which show a maximum near T . 

g 




PART G 


THE PERCOLATION MOREL 
OP 

SPIN GLASS TRANSITION 



GHAPOBE 6 


SPIF GLASS TBMSITIOII AS A EBHOOLATIOU PEOBIBM 

In this chapter we develop the basic framework of 

our percolation model [351]. let us first describe the 

'mechanism' of formation of 'locked' spin clusters, A free 

spin in contact with a heat bath, maintained at a temperature 

T, has a finite relaxation time t^ . For CulSn an estimate 
-13 

of t^ = 10 sec. has been obtained from remanence data 
[ 352 ] . The presence of a second spin interacting with the 
first has a retardation effect changing the relaxation tine. 
So long as J <^ 3 ^? these two spins can be regained as 
'independent', albeit with lifetimes longer than a free spin. 
However, if J>kgT, then in whichever configuration we begin, 
the spins eventually 'align' (either feiro magnetically or 
antiferromagnetically depending on the sign of J) and from 
then onwards move as a single entity locked together. We 
may, therefore, define the locking criterion as J_> Q-nd 

hence the locking temperature 

"'1 = JAb 


Contents of this chapter have been published in 
J. Phys. P13, 431 (1983). 
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Por PiKKY interaction, the locking condition becomes 

AR ^ exp(-^/6 ) > kji (6.1) 

the distance between the spins, and 
i Cos(2k R) 

h k-tj! 



m* ’ 

= volume of unit cell = V 

2 ° 
h kj 

~ * j 

2m 

= Perm! momentum, and 
= effective mass of electrons. 

The exponential damping factor arises due to the scattering of 
the conduction electrons in a disordered system. If the 
dominant scatterer are the spins themselves, rather than 
other defects the damping is called self damping, and 

6 cc c , 

On the other hand, for damping introdubed by defects, 6 is 

independent of concentration [ 14^ , We assume that our 

3 

system remains metallic and 5 >> V . 

This criterion (6,1) itself is insufficient to 
define the locking of two spins [ 4C3 . We have not yet taken 


where R is 


and 


a" 


'P 


k 


P 


m* 
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frustration into account. In this chapter, for the first 
tinB, effect of frustration will be incorporated into the 
percolation model of spin glass transition. 

Because of the oscillating nature of the EKKY inter- 
action and random occupation of the lattice sites by the 
magnetic constituents, the geometry of the locking bonds 
inside a locked cluster may lead to frustrated plaquettes. 

The frustrated bonds weaken 'lo eking ^ because the terminal 
spins, due to frustration, can still flip independently ewen 
if condition (6.1) is satisfied [4-0]. Moreover, a set of 
'loose' spins may also exist even inside the locking sphere, 
if the local field that 'pins' the spin vanishes at some sites. 
We shall assume that such spins constitute a negligibly small 
fraction. This assumption will be justified in Chapter 11. 
Thus, in addition to condition (6,1) our criterion for the 
formation of a locking bond involves the examination of the 
bond topology locked according to condition (6,1), and the 
exclusion of all frustrated bonds from the cluster. 

Bext we geometrize the problem exactly as done by 

Smith [ 353] • Attached to every magnetic site is a sphere 
1 

of radius where R is determined by the equality in 
condition (6,1). Two sites are bonded by a locking brnd 
if their spheres overlap, and the environment of the bond 
is such that it does not form part of a frustrated plaquette. 

M n- cluster consists of a set of n such sites, each site is 
bonded to at least one other site within the set and no 
site is locked to any other site outside the set. In such 
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clusters, unlike ferromagnetic or antiferromagnetic clusters, 
all the spins do not necessailly align parallel or anti- 
parallel, but depend upon the sign of the local magnetic 
field. If we try to flip a spin within a cluster, all the 
other spins in it try to flip with it. Since the sites may 
be distributed a.t random, amorphous as well as crystalline 
substituional alloys may be taken into account. 

As the temperature is reduced, the 'clusters’ grow in 
size^ more and more smaller ’clusters' coalesce together to 
form bigger clusters. At a critical radius, (i.u, at a 
critical temperature T ) an infinite cluster of locked spins 

s 

first appears and percolation is said to have taken place. 

T is identified with the paramagnet to spin glass transition 

O 

temperature. Such percolation process, where the concentration 
of sphere centres is kept fixed while the radius is varied, 
has been studied numerically by Kertesz [ 354 ]. Other kinds 
of percolation of spheres have been reviewed in our paper 
[ 351 ]. 


6.1 versus c; 


The crucial quantity of interest in the study of 
this percolation is the volume fraction 


p(T) = 1 - exp a ] 


( 6 . 2 ) 


where n is the number density of participating (non frus- 
trated) sphere centres. Let n^ be the number density of all 
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atomic sites and c' be the fraction of these atoms that 
are magnetic and non- frustrated. Then n = ° is 

small, frustration is small so that c' = c. As c increases 
but remains below c^ (percolation threshold for long-range 
magnetic order) frustration increases and c' becomes pro- 
gressively less than c. We shall examine o' = f(c) as a 
function of c subsequently. Prom equations (6.1) and (6.2) 

AE'^(Tg) exp[ -E(Tg)/6 3 = k^Tg (6.3a) 

and 

y = E^c' (6.3b) 

where 

y = Tfly - p 

Eliminating R from (6.3a) and (6.3b) we get T as a function 

o 

of c. We shall examine the relation between T and c in 

O 

the following special cases. 

Case (i) : Infinite Mean Pree Path and Very low Goncentration : 
In this case 

AR"2 (Tg) = kjTg 

7 = 

6 -*00 and c' = c. Hence y = Ac/k^T^ implying 

This is in agiaement with the numerical data of 


and 


be cause 


c, 

O 
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Kiazel [354] and the scaling ideas of Souletie and lournier 
[ 35^ . In fact, the scaling laws of Souletie and lournier 
follow from two facts — “4“ law of interaction and c -♦ 0 

[l08] . Both these facts have been utilized in case (i) 
and hence the agreement is not so surprising. 

Case ( ii ) i Infinite Mean Bree Path and o' <c ; 

In this case, exactly as in the earlier one 

c’ . 

O 

Thus frustration lowers T^ for a given c, because, less 
magnetic sites are available for cluster formation; conse- 
quently I^(Ig) hs-s to be larger or lower than the value 
when frustration is absent and all magnetic sites participate 
in cluster formation. 

Case (iii) s Ifo Prustration (c* = c) and Finite Mean Bree Path ; 
Now, 

^ exp [- y''/3 

— 1 

Hence for self damping, (i.s.s^c c ) 

T <rc c exi[ - b 

s 

Thus, both frustration and finite mean free path lead to a 
weaker than linear concentration dependence of T . However, 

O 

there is no evidence of a unique power law of the form 

T cc c^ which holds across the whole concentration range, 
g 
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Moreover, the functional dependence is not universal, as 
c' and 6 both depend on the structural and chemical nature 
of the system under consideration. This is contrary to the 
universal T ec law obtained in the mean field theories 

of Edwards and Inderson [ 114] . Real experiments on canonical 
spin glasses yield p varying between 2/3 to 3/4 [73] . Thus, 
our prediction turns out to be closer to experimental 
results. 

From equations (6.3) we may write 


Tg (c,6) = TgCcjOo) exp[-E (c,Tg)/6] 


where E(c,T ) is the critical radius for a given concen- 

O 

tration. This immediately gives us a direct physical 
interpretation of the ‘characteristic length* 


0(c) = - 


d(ln T^) 
d(l/6 ) 


(6.4) 


of Kinzel [354] and Kinzel and Fischer [ 356] . 

Using p(T„) = 0.28 as obtained by Kertesz [354] 

s 

and equations (6.1), (6.2) and (6.4) we get 
n ^ 

(' s l S) o' = 0.0784-2 X ^ (6.5) 

® [ 0(c)/r^ 5 

For a given c, we put /(c)/r^ for QuMn from Kinzel 's data 
into equation (6.5) and evaluate c'. Figure 22 shows c‘ 
plotted against c. The behaviour confirms our earlier 
analysis . 



- 0.01 



Fii’ur^ 22 : Plot of o' i^c) against c, ODtained from tbe 
interpretation of the numerical data of -r 
Kinzei f354-J . The units of c* are (n^T^ if/S 
and are explained in the text,*^r^ afid 
are as used by Klazel. 


m m 
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6,2 Frequency Dependence of 
Yogel-iPuXoher Daw; 

Almost all types of glasses (including spin passes) 
obey Vogel-Fulcher law [8] 


r = Tq exp ^ C^-^) 

gv j o 

rather than Arrhenius law 

t: = Tq exp [-^] 

where is independent of the frequency of measurement, 
although T depends on it crucially. The latter dependence 

o 

will be studied in this section. Since % diverges at 
Tg(ci>) = T^, is a true transition point revealed by a 
static measurement. Hence, we identify T^ with T^ in 
section 6,1, Thus all cooperative effects will manifest 
themselves via T^ (i.e., T^), 

We now derive (6.6) from the basic model developed 
in section 6.1. let us concentrate on the cluster, say 
ith, that grows to infinity at T , Hence the number of spins 

o 

in the ith cluster should be given by 

n^ - (T - Tg)” ^ 

which diverges at T = T , as required. On the other hand, 

e 

the life time of this cluster can be written as 


'z = X ^ exp ( IFj^l/kjT) 

where is the free energy of this ith cluster, i.e., 
is the free energy needed to overcome the barrer. 
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Since the cluster is in thermal equilibrium with 
the free spins at any temperature the free energy per spin 
within the cluster is In (2S+1), so that 

X = 'V ^ exp [ n^ In (23+1)] (6.7) 


i.e., 


T = T ^ exp [ 


A 


(T - T) ' 


] 


where , 


A = (const) ln(2S+l). 


Pora.c, susceptibility measurement, the cluster will appear 

frozen as soon as t; > where v is the frequency of 

m 

measurement. Thus 


V 


m 


^o 


exp[- 


A 


1 T„(a)) - T^(0)}Y 


•] 


( 6 . 8 ) 


We call (6.8) Vogel Fulcher type law; for Vo gel- Fulcher 
law Y = 1. From experimental measurements Tholence [301] 
found Y = 1 , although he maintained that more high-frequency 
measurements are required to either confirm this or obtain a 
more reliable walue of Y • I'or Be the lattice y= 1 [3531 . 
Ohalupa [357]} on the other hand, found y = 15 . However, 
his ideas were different from ours. He maintained that the 

■z 

number of clusters is q/ s ( correlation length) or 
that clusters have a characteristic volume. This is not 
true in our percolative picture. Cluster sizes in our 
picture grow inhomo geneously and we have a distribution of 
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cluster sizes and their corresponding relaxation times. This 
fact will be utilized in a later chapter. 

Shtrikman and Wohlfarth[ 296] derived Yo gel- Pule her 
law within their model of blocking of superparama^etic 
clusters. In the strong coupling regime, i.e,, when inter- 
cluster interaction is not negligible, Shtrikman and 
Wohlfarth had to assume temperature-dependent effective 
cluster size. This idea seems to be closer to ours. The 
relation between Shtrikman-Wohlfarth approach and ours will 
become clear in the next chapter. 


6,3 


q-Dependence of T 


£ 


Neutron scattering data in spin glasses were analyzed 
by Murani [37 >8 5] to show that x(q) exhibits a maximum at 
Tg(q) which decreases monotonously with increasing q- value. 
Bhargava and Kumar's C358] treatment based on EA model yield 

q-independent T . levin et.al.'s cluster model C 220-2223 also 

s 

yields a q-independent T . The latter work has been severely 

o 

criticized by Murani [223] (see section 1.4.2 of Chapter 1 for a 
discussion). On the other hand, MurardI 37>85] argued that 

i 

at finite q the characteristic length scale probed is^-^ 
and hence only clusters of size _> l/q will appear frozen. 

At Tg, v;here an infinite cluster first forms, there are still a 
finite fraction of spins in clusters of size l/q which have 
lifetimes smaller than the measurement time Only at 

lower temperatures all these clusters join the infinite cluster 
in being frozen. This would lead to a Tg(q) lower than 
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Tg. (cL=0)* this section we shall quantify Murani's ideas 

o 

using the standard scaling laws of percolation theory [359^* 
Let C(q) be the correlation function of our 
percolation theory. This is directly related to the ’pair 
connectedness’ function of Essam [359]. G(£-r',p) is defined 

as the joint probability of 

(a) the point r being occupied by a sphere centre which 
belongs to a finite cluster and 

(b) the point x Q-lso being occupied by a sphere centre 
77hich belongs to the same cluster. 

p(T) is the volume fraction. If P(p) be the percolation 
probability i.e., the probability of a sphere joining an 
infinite cluster, then the probability of the first event 
mentioned above is p„ = p [l - P(p| . Moreover, defining 
'0 ( = 1 if £. and are coniiected in the sense of 

(b) and zero otherwise . IJnen 


G ( £ - £' » P) = Pa<'n 

Uso, we have 

n. (p) = IKp (r - r' )la >= -T— S G(E,p) 

-^a R 

where, 

R = £ - ^ 1 ' . 

Let us define C(l,p) as the Pourier transform of C(R,p)? 


Pa =0 (0,p). 


then 



(6.9) 


The life time of a cluster is 

■ G(1jP) 111 (2S+1) 
T = r exp[ 


] 


as seen by a probe of characteristic length l/q. As q-®-0, 
C(£»P) 

— r -* n. and (6.9) reduces to (6.7) which was valid for 

^a ^ 

q=0 probes. 

Note first that Pg_ is a function which does not 
diverge at any temperature, from the scaling ideas 


C (£, p) = ^ ^ g(q I ), 


where g is a universal function for problems in this class, 
with g(o) equals constant and g(r) = tor large 

values of r . This gives c(0,p ) = g cc. (T-T , 

diverging, as expected, at T : and, C(_2.,p) = g 

O vAy 

for p = p.. If 'C™ is the characteristic time of our probe 

W 111 

then X (£) has a maximum at T ( £) which is a solution of 

g 


In ( v „/ ) 

r '■ m' 0 ^ 1 

In (2S+1) 


G(a,p) 


P 


( 6 . 10 ) 


a 


With the weak temperature dependence of p^^, near T^ the 
right-hand side behaves as g^q whereas it ^es to zero 
as T 00 ; in between it is presumably a smoothly decreasing 
function. The right hand side of equation (6.10) is plotted 
as a function of T for various values of q in Figure 23, 

From the intercept shown in the graph for the solution of 
equation (6,10) it is clear that T (£) decrease with 

O 

increasing q, in agreement with experiment. 



3 



Scbeniatic plot of C (ii»T) as a function of T for 
various of q fjie intercepts iniicate 

solutions for r (q)* The intercept of th© 

horizont^ line ocoSrs at la (T,,/t )/ln(2St1) 

(point A‘)» At O' ' 
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6 .4 Linear Response Theory and Zero- Field L ov^- Frequency 
AG Susceptibility ; 

Since zero- field susceptibility is an idealized concept 
and a finite, although small, external magnetic field is needed 
to measure susceptibility, we shall use the terms low- field 
susceptibility and zero-field susceptibility interchangeably. 

A well known result of linear response theory (LET) 
is 

X(T) = iSs 

- [< S3^(t)Sj(o)>] ) 

t**co (6.11) 


which expresses the response function, namely, low frequency 

a. c. susceptibilily X(T) in terms of fluctuations. Here 
HS(S+1) 

0 = is the Curie constant. Neglecting short 

range order, 


L 


[<Si-Sf ]aT= S 6i,3 = »S{S+1) 

Ij 


ay 


Similarly, 


( 6 . 12 ) 


lim 
t -00 


S 

i3 


[<Si(t) S.(o)>j 


av 


S P(T) Q(T) = FP(T) Q(T) 

(6.13) 


where P(T) = fraction of spins included in the infinite 
cluster in the thermodynamic limit, and 


JL 

NP 


sr<s.(t) Sj_(o:^l 


i^in j. 

clus. 


ay 


Q 


lim 
t - OD 


(6.14) 
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is the EA order parameter. Notice that while carrying out 
the spatial averaging in (6.14) we have averaged over the NT 
sites included in the infinite cluster instead of all the N 
sites. This is so because only the spins in the infinite 
cluster have non-zero autocorrelation in the limit t -* oo-. 
Substituting (6.12) and (6.13) into (6,11) we get 



Similarly, in case of classical spins of unit ’length' , 
i.e., S^ = 1, we get [351] 

X(T) = I [ 'I - 5(2)) Q(T) ] (6.16) 

It is often argued that the disagreement between 
Fischer’s and Pari si ' s expressions for X(T) indicates that 
LET breaks down in spin glasses because Fischer used relation 
(6.11), a standard result of LET. But one should notice 
that in addition to LET Fischer also used homogeneous MFT of 
cooperative phase transition model (GPBu) (i.e. P(T) =1 
for all T <_Tg). In our model LET holds good although GPTM 
was assumed to be inadequate [ 360] . 

In a finite- sized spin glass, the lifetime of the 
biggest cluster, even below T , will have a finite (although 

o 

large) lifetime Hence, in such cases P(T) =0 

for t shall have X(T)= -^i.e,, a Curie-like 

behaviour. Such behaviour of susceptibility in spinglasses 
have been observed in computer simulations on finite 
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systems [24,561 ] . In other words, spin glass state is not 
an equilibrium phase in the static sense but only in dynamic 
equilibrium for t << So we should write (6.16) as 

'^'de^^^= I [1-P(T) Q(I) ] (6.16a) 

On the other hand, all the spins (P=1) appear frozen to any 
probo of rery short characteristic measurement time for 
all I < Hence the ideal non-equilibifum susceptibility 

will be given by 

'Vxe'^^^= f [ 1 - Q(T)'] (6.17) 

which is identical with EA result. 



GHAPOISR 7 * 

HELATIOI? BBOMEEIT PERGOLATIOH MODEL , 
^A 'MOPgL MB WOHEiPABTH MODEL 


In this chapter the mathematical relation between 
our percolation model and EA model and that between perco- 
lation model and Wohlfarth model will be established and 
interpreted physically [362]. 

7.1 Relation Between Percolation Model and EA Iitodel 

Putting P(T) =1 for all T < I in equations (6.16) 

o 

and (6.15) we get EA and Eischer‘ s results ( 1.27a) and 
(1.30) respectively. Phis is expected because P(T) = 1 
for all T T implies cooperative freezing of all spins 

O 

at T which is, indeed, the physical picture of the SA 

o 

cooperative phase transition model 

7 .2 Relation Between Percolation Model and Wohlfarth Model 

In the Wohlfarth model the magnetic susceptibility 
at temperature T is given by [ 287] 

X(T) = ^ f f(T') dT* (7.1) 

o 

Various parts of this chapter have been 

(a) published in J.Phys. P13> D19 (1983) 

(b) submitted for publication. 
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where G is the Curie constant and f(T) is the distribution 
of blocking temperatures, i.e. f(T) AT is the fraction of 
blocking temperatures available in the temperature interval 
between T and T + At. 

The physical meaning of (7.1) is very clear. At 
temperatures T only those clusters whose blocking temperatures 
Tg < T will contribute to the susceptibility, i.e,, those 
clusters which have not yet ^frozen’ . Alternatively, rewriting 
(7.1) we get 

00 

X(T) = § - I / f(T') dT' (7.2) 

since according to the definition, f(T) is normalized, i.e., 

00 

/ f(T) dT = 1 (7.3) 

o 

If all the clusters were unblocked, as is true at high 
tenpe natures only, X(T) would be equal to the first term in 
equation (7.2). However, at lower temperatures the second 
term, which represents the contribution of blocked clusters, 
has to be subtracted. The aim in this section is to find 
the relation between (7.1) and (6.16) . 

let us, for the time being, assume that 

f(I) = - ^ (7.4) 

This is the main result of this chapter and will be justified 
very soon. Substituting (7.4) into (7*1) we get 

X(T) = I [ P(0) Q(0) - P(T) Q(T)3 


(7.5) 



But P(T=0) = 1, according to definition, and Q(T=0) = 1, 
for classical spins of unit length. 

Hence (7.5) reduces to 

X(T) = f [ “I - Q(T)] 

which is identical with (6,16). 

Let us now interpret equation (7.4) physically. If 
there is a large (small) fraction of blocking temperatures 
f(T) AT in the temperature range T to T + ^ T, then as one 
goes from T + A i to T the change AP in the function P 
will be large (small) because a large (small) number of spins 
will coalesce with the infinite cluster. Haturally, the 
question arises %B to why we should take the change in the 
product P'5, i.e. A (PQ) rather than sicgjly A p in equation 

(7.4) . This is not yery susprising. Although nonzero P 
’indicates' spin glass ordering, it is insufficient as a 
'measure' of the order. Any system of 'particles' inte3>- 

i 

acting yla — ^ type attract! ye interaction may undergo a 
r-^ 

percolation transition of the type considered in Chapter 6 . 
However, such a geometrical transition does not necessarily 
inply a spin glass type ordering telow the percolation 
threshold. In other words, EA order param.ter Q is a 
measure of this ordering (some sort of spin stiffness). 

There foie, one should take the product function PQ as a 
measure of the spin glass ordering. The negative sign in 

(7.4) accounts for the fact that PQ increases with decreasing 
temperature , 
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Por quantum spins S, (7.4) gets modified to 

f(T)=- (7.6) 

Proceeding as before and remembering that Q(!r=0) = S(S+1) in 
this case (see Appendix 1) we get from (7.1) 


:(T) 


Or. P( T.)_q (P) -j 


s(s+i 


which is identical with (6.15). 

‘The results Q(T=0) = 1 for classical spins of unit 
'length' and Q(T=0) = S(S+1) for quantum spins were derived 
in appendices .within the framework of EA model. But since 
r(T=0) = 1, the results of BA treatment and percolation model 
should be identical at T=0. 

Another similarity between the Wohlfarth model and 
percolation model is the fact that the former even allows for 
a temperature-dependent cluster size in order to be consistent 
with the observed Vogel- Pulcher law [296] , though it does not 
bring percolation into the picture. 

It is very interesting to note that an equation 
exactly similar to (7.4) has been derived by Murani 
for the distribution of relaxation time ^ as 

g( ( 7 . 7 ) 

where £-(u ) is the time- correlation function in the limit 
momentum transfer q -* 0, and g('c) is the distribution of 
relaxation times. Since there is a one-to-one correspondence 
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between blocking temperature and relaxation times, we think 
that equation (7*7) supports our view. Moreover, is 

related to g ( t) by the equation 

X.(T) = f- 7g(ir')at' (7.8) 

o 

which is formally similar to equation (7.1j and can be similar- 
ly interpreted. Since only clusters with relaxation times <% 
contribute to x(r) the upper limit of the integral on the 
E.H.S. of (7.8) is t * 

It appears at first sight that since even finite 
clusters can freeze in Wohlfarth model it will invalidate the 
derivation of equation (6.15) and (6. lb) from (7.1) But it 
is not so. finite cluster freezing iirplies B(T) 0 even 
for But since B(T_) has not been utilized in the 

o 

above derivation our treatment remains valid. The only 
modification necessary for finite cluster freezing is to 
call P(T) in equations (6.15) and (6.16) the 'fraction of 
spins blocked' instead of 'fraction of spins belonging to 
the infinite cluster' . We summarize this section saying that 
it is the percolating spin-clusters that mimic the 'grains 
fins' o.f the Yifohlfarth model so convincingly. But unlike 
the Wohlfarth model, spins belonging to different clusters 
interact via EKKI interaction as described in Chapter 6. 





7*3 Calculation of f ( 7) [ 3631 

Although Wohlfarth model does not pro-vide any 
prescription for abinitio calculation of f(l), the latter 
can be obtained from experiiaentally measured -values of 
using relation (7*1)j i.e., 

f(T) = (T><) (7.1a) 

T 00 

Van Duyne-veldt and Mulder [288] calculated f(l) from their 

susceptibility data using relation (7.1a). Surprisingly, 

f(T) thus calculated exhibits a sharp peak, near T , the 

a 

temperature corresponding to the susceptibility maximum, 
krom this observation Shtriki^ian and Wohlfarth [ 296 ] concluded 
that such apparently cooperative effect around T arises 
due to inter-cluster interaction. It is difficult to account 
for the observed behaviour of f(I') within the original frame- 
work of non-interacting clusters of the Wohlfarth model. In 
this section vre calculate f(T) by an alternative method 
comparing the experimental data on susceptibility with our 
theory (i.e., using equations (6.16a), (6,17) and (7.4j). 

As a result we gain an insight into the origin of the 
cooperative effects in the gradual freezing of spin clusters. 
This coopei'ati've effect arises from the intercluster inte 2 >- 
action which is taken into account in the percolation model. 
We can rewrite equation (6.16a) as 

^de = f (1 - Q + A ) (T.8) 

(which apparently looks similar to Sompolinsky* s result [117]) 
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and 


A = (1 - P) Q 


(7.9) 


On the other hand 


X 


ne 


(1 


Q) 


( 7 . 10 ) 


Yeshurun and Sompolinsky (YS) [364] found the same express- 
ions (7.8) and (7.10) for and X^^ although the micros- 

copic description (i.e., Sompolinsky ' s model [liyl) underlying 
their derivation is different from ours. However, this 
difference does not affect our reasoning here because neither 
YS nor we calculate h and Q from microscopic considerations. 
Using equations (7,3) and (T.10) and the well known fact [76] 
that a field-cooled spin glass achieves its equilibrium state, 
YS obtained A and Q from their experimental data. Yrom 
these A and Q we calculate PQ (shown in figure 24) 
using relation (7.9), i.e. PQ = Q - Using PQ thus obtained, 

we calculate f(T) from equation (7.4). f(T) plotted in 
figure 25 is in remarkably good qualitative agreement with 
that obtained by VM. [288], from PQ and Q values we calculated 
P as shown in figure 24 . It is typical of any percolation 
probability s the probabilily rises very sharply near the 
percolation threshold and then flattens off. 
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• PQ 



Figui© 24 J F, PQ and Q shown as a functtort of tempara-tore for 
amorphous SG (5'®o,64 *^0.36^75 ^16 ^6 ^3* 



Vi) 



Pigur© 25 * IJlstrlbution of blocking temperatures f{f) for 
amorphous SG (^^0.64 *%.56^75 ^16 % 



OHAP!I!BB 8 


SLOW HBLAXATIOU OB MAGMETIZATIOl 
IH SPIB GLASSES MD 1/f lOI^ 


Our aim in this chapter is to derive the logarithmic 
relaxation (initially exponential for veiy short time) of 
magnetization analytically within the framework of the perco- 
lation model of spin glass transition developed in the precee- 
ding chapters [ 365-368]. It reveals a deeper connection with 
l/f noise in a wide variety of other systems [366], 

<1 

8.1 Relaxation in the Absence of External Magnetic Field 

logarithmic decay was derived earlier in other similar 
contexts by Reel [ 369] , Street and Woolley [370] and Gaunt 
[371]. In case of rock magnetism Reel [369] assumed the 
existence of a distribution of the volume of magnetic 'particles^ 
and coercive force while Street and Wooley [ 370 ] assumed a 
distribution of barrier heights. Both these approaches are, 
of course, physically equivalent. A distribution of relaxation 


Various parts of this chapter have been 

i) accepted for publication in J.Phys.B (in press). 

ii) accepted for publication in Sol. State Gommun. (in press). 

iii) accepted for publication :in Phys.Lett.A '^in press). 
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times appears very naturally also in our percolation model - 
it is associated with the distribution of cluster size. 

The distribution P(n) of cluster size n, i.e., the 
number of finite clusters containing n spins can be written 
as 


iP(n) 


IT 


5xp (- /,n) 


( 8 . 1 ) 


5 3>- 


where is a normalization constant , X depends on tempe; 
ature and \x is another positive constant [ 365 ] . Value of p, 
for various lattices are well known [ 372 , 373 ] , ha(T-T^)^. 
Then, for normalization 


~ a ' 
o 


n. 


■jnaz 


n 


max 


n 


^ 1 


, oay, 

' ( 8 . 2 ) 


where n^ ^^y is the size of the biggest finite cluster and H 
is the total number of spins in the system. Putting An = s, 
(8.2) reduces to 


' o 


(p, -2) ^max _s (i_p) - 

'■ / <3 as = Ns 

K 


and hence 


11 


= /[ ! X - ^ V'V2 


*( 1-n)/2, (2^ (/a (X } ] 

where (X) is the Whittaker function, the functional 

argument being X. Since is usually very large, replacing 

infinity introduces negligibly small error. Then 


' (1-u) (2Ht) 

2^ 2 


(>^) 
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2 2 

Let M(0) be the initial magnetization per site and 
i;(n) be the relaxation time for magnetization per site in 
a cluster of size n. x (n) dept^-nds on the cluster size n and 
all sites in the same cluster have same x, because when a 
cluster relaxes, all the spins belonging to it relax together. 
In other words, relaxation time of magnetization per site in 
a clus fcer is identical with that of the cluster. The relax- 
ation time of a cluster of size n^ derived in Chapter 6 
(equation 6.7), is given by 


whe3>3 , 


t; (n) = xo exp(uK) 

E = In (2S + 1) 


(8.4) 


Because of the reilation (8.4) between n and x , a distri- 
bution iP(Q) of cluster size will always be equivalent to a 
distribution g(T^) of relaxation times. 

The magnetization per spin observed at time t 


2 magnetization of the ith spin at time t 
i 

n 

max 

2 (magnetization per spin at time t in a 
n=1 cluster of size n) nfi?(n) 
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where 


'to N 


i £ infinite 


M. = M(0) P 


cluster 

because the lifetiioe of the infinite cluster is infinitely 
large . 

Therefore , 


n 


M(t) = ^ 


max 


1 


n5’(n) exj£ -t/tCn)] dn + I'a(O) P 


n 


M(0) IT / 
o 


max 


CO 


exp(-Xn) s-^ . 

dn — b' . ) - -!-../ .. ^ e 

1^“^) p=1 P- 


1 n 

+ M(0) P 


S e"P-^ } 


—1 

where II = II /il and a = t ~ 
^0 o 


Using (8.2), 


M(t) = M(0) + I 


where , 


I = H(0) II 


00 


P=1 


p'. 


n 


max 


s exp[-(X+ pK)n 




m(0) II, 


CO 


P=1 


(- at) 

pi 


P 


(X+ pK) 


( 1 - 2 ) 

(X+ pK) 


-u 


u 




where , 


Hence , 


u = (h + pK)n 


I = M(0) II 


00 


p=1 




pi 
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[ (A+ pK) 


.fcll 


2 

2 ’ 2 


_ (ti-1) _ (^+ pE) n 

- {(?v+ pK)n^a^} ^ e ^ 


max 




(1-H) (2-|i) “-max’ ] 


Since is usually rery large we can, as before, drop the 

second term within the brackets. Using the asymptotic form 
of the Whittaker function for large ( h + pK) we get 

00 

I .t! M(0) it' ^ (^+ pK)-"* 

° p=1 P* 

We further assume that pK » X ,so that 

^ M(0) {-d:/(2S+1)}^ ^ 

I ^ __ _ _ 

M(0) n' ' 

= — K ^[Ei {-a t /(2S+1)} -0.577 

- In { a t/(2S+1)}] 

where Ei(z) is the exponential-integral function and 0.577 
is the Euler's constant. Pinally, substituting the series 
representation of the function Bi(-r z) we get 

M(t) = ]Vi(0) + A + B exp{- at/(2S+l) ? -S(l) In t 
where , 


(8.5) 
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m(0) N 

^ ^ _ o [0.577 + lnfo/(2S+1) 1 ] 


aiid 


S(T) = 


_ 


are temperature-dependent but independent of tlij^e, and 


B 


M(0)B 


K 


B (- 1 g 

k {at/{2S+'l)r 


Therefore, the relaxation is quite fast (almost exponential) 
initially because of the dominance of the third term in (8.5). 
However, this third term decays very fast. Therefore, on 
lotiger time scales the last term in (8*5) dominates and hence 
the logarithmic relaxation. Bor large t (8,5) reduces to 

I«I(t) :: M(0) - S(T) In t (8.6) 


In addition to the logarithmic relaxation few more results 
emerge ; 

(i) We know that u < 3 in general. Then 


S(T) = 


M(Q) A' 


(3 - v) i'-/2 S 



S(T) = 0 at T = T where ^ = 0, and S(T) at T=0 wtere 

S 

T 0. Hence, S(T) is expected to have a maximum in 

between T=0 and T=T„. Such behaviour of S(T), has, 

g 

indeed, been observed experimentally [76], 
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(ii) Since finite clusters exist for T >1^? one would 

O 

expect logarithmic relaxation not only for ’I <T , but 
also at higher temperatures. Such logarithmic (or 
almost logarithmic) relaxation even at T ™6 T has been 
observed in monte carlo simulation of canonical spin 
glasses by B'ernandez and Streit [374]. i'or still 
higher temperatures the relaxation should be exponen- 
tial because only single spins (and no bigger cluster) 
vTith single relaxation time exist at such high temper- 
ature s . 

(iii) One may raise objection against the form (8.1) of iP(n) 
because it is monotonically decreasing function of n 
whereas one would expect {P(n) to exhibit a maximum at 
an intermediate value of n. Best fit with computer 
experiments is obtained with an expression of the form 

exp (d- '/‘n - d^n) 

IP(n) = H 3 2_ 

° nP 


where d. ( T - T ) and ( T - T ^d is a 

positive constant. Sinco at 1 , jP(n)<^ n~ ^and 

O 

^ itP(n) = p must be at least of the order 0(1). 
Por convenience we choose p = 1 , so that 


{P(n) = h. 


. r f 

exp ( a'/^ n - b^ n) 


( 8 . 7 ) 


! 2 

where, 4 = (T - T ) , and 

& 

= If b 4 exp (-aV 8 b)/D_ 2 (“®’/'^)» 



D]^(z) being the parabolic cylinder function. Proceeding 
exactly as before we again get p67 ] (8.5). This is not 
surprising in view of the approximation pk >> X adopted 
in both cases. 

(iv) Por magnetic systems Ohikazumi [ 375 ] considered an 
example to show how a distribution of relaxation times 
gi'c ) leads to slower- than-exponential relaxation. 
Ohikazumi 's distribution 

-~r- for rg >x>t^ 

= 0 otherwise (8.7) 

can be derived using (8.1) and (8.4) assuming X= 0, 
p, = 1 and that 

= ln( ) = const f or ^ < ■^< 

= 0 otherwise. 

(v) Let us now estimate the time-interval, in terms of 
within which (8.6) holds. 

Notice that the basic characteristic time scale in the 
problem is because the latter is the relaxation 

time of a ‘loose' single spin in the spin glass. One 
must remember that 2 is the renoimalized relaxation 
time of a 'loose* spin (a cluster of size unity) 
weakly interacting with other spins via damped EKKY 
interaction rather than that of an isolated 'free* 
spin. Since the exponential term in (8.5) is negligible 
for t>2^, equation (8.6) holds for 
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t > t. 


min 


2t:_ 


( 8 . 8 ) 


On the other hand, since we started with positive M(0), M(t) 
cannot go negative even after long tinB. Since the equili- 
brium value of M(i.e., M(t = oo)) is zero, (8.6) can hold 
only up to a time t =t^^^ such that 

ULcXjL 


M(0) - M(0) n) K la 0 


In other words, (8.6) holds in the ti roe-interval t< t^^^ 

where t^j_v 3 _ given by (8.8) and (8.9) respectively. 

(vi) notice that no frustration has been explicitly assumed 
a priori in this calculation. Hence, it follows that 
frustration is not essential for logarithmic relaxation 
in spin glasses. Medina et.al. [376] also arrived at 
the same conclusion from computer simulation. 

(vii) We now explore the relation with l/f noise. It was 
Kogan [27^ who first pointed out that because of the 
logarithmic relaxation of magnetization, the noise 
spectrum in spin glasses is proportional to l/f exactly 
like l/f noise in semiconductors [276], dielectric 
i-elaxation [ 277] , dislocation mobility [377], etc. 
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Historically speaking, Schottky, using exponential 
correlation function (like ours in Chapter 4) with 
a sin^e relaxation tiioe, obtained a lorentzian that 
in the large f limit yields 1/f spectrum rather than 
l/f spectrum obser-ved experimentally. Later, a class 
of models for l/f noise was proposed based on a 
distribution of relaxation times. The situation in 
spin glass physics is clearly similar to this. Besides, 
kacfarlane [ 378 ] assumed distribution of 'cluster size' 
wliich is physically similar to our approach. 

(viii) Although our approach apparently looks similar to 
that of Stauffer [ 37^ there ar<^ some fundamental 
differences. We have considered oscillating EKKY 
interaction whereas Stauffer considered nearest 
neighbour ferromagnetic interaction and hence the 
difference in the definition of 'cluster' in the two 
models. In Stauffer's model, at very low temperatures, 
all the spins belonging to the same cluster orient in 
the same direction whereas in the percolation model, 
as discussed in Chapter 6, spins got pinned along 
the spatially random local field and hence may orient 
along different directions even in the same cluster. 
Besides percolation in Stauffer's model signals the 
onset of long range order whereas that in our problem 
leads only to a random freezing of local magnetization 
without any long range order in space. 
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8.2 Eelaxation in the Presence of External Magnetic Pie Id 668] ; 

As shown in Chapter 5 mean field approximation with 
single relaxation time is inadequate to explain the slow 
relaxation in the presence of external magnetic field. Using 
the same argument as in section 8.1 we now hawe 


(t) = u'/[M (eq) + {M (0) -k (eq) nP(n) dn 


( 8 . 10 ) 


(8.10) yields, for no t-'teo- short time scales. 


f " t I 


M (t) =[ivl (eq) + K] -[K U /ln2]lnt 


( 8 . 11 ) 


with E = ^^(O) - M^Ceq). 


Eemembering that it is that is measured experimentally, the 
form (8.11) for magnetization relaxation has been observed 
by Salamon and Tholence [ 380 } . 



GHAPOER 9 


RESISTIVITY OP SPIR SLASSES IR IHB PERCQLAOROR MODEL 


In this Chapter we shall utilize the distribution 
(8.7) for relaxation tiroes to calculate the resistiTity of 
SG within the framework of the percolation model. Experi- 
mental data near T tum out to be in better agreement with 

o 

the expression to be derived in this Chapter rather than vvith 

(4.24) . However, for T »T„ (say T^ 6 T ) the expression 

o o 

(4.24) remains valid though the approximate form (4.25) 

0(T,'i;) no longer holds. 

In Chapter 4 we derived the following expression 
for resistivity b(T) with single relaxation time t ; 


P(T) = 


m*c 


■5 


\ir 


ne 




4SI|±1) . Q(T)iM^] ]j 

(9.1) 


Since in the percolation model only a fraction P(T) of 
spins have non-zero long-time (t oo ) autocorrelation, 
equation (9»1) gets modified to 

P(I) = Q(T)E(« 

ne ^ 

/ g(T;) oL t; ] J (9,2) 

TC X 

where g(v) is given by (8.7). This can also be checked 
from earlier steps (4.25) - (4.27) where we replaced 

Submitted for publication. 
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"I li. 

^ (t = oo) by cQ« The latter replacement was valid 

in the cooperative phase transition model where all the spins 
have non-zero autocorrelation for all T < T^. This is not 
true in our percolation model as stated in Chapters 6 to 8. 


Wow, 


f g(T;) 

% X 


, 2 ri.2 2 

-ji ^ 

+3) 




T 


ln[(-n: kgTT ;2 + 3)/(Tk-gT x^ + "5) 
p ln('i:2 / x^) 


(9.3) 


Substituting (9.3) into (9*2) we get 

P(T) = [1.7^ + 3 [ QP +C ■ § -L |± l I « Qp} 

ne ^ 

In { (TiCgTr^ + 3)/(7ckgPr^ + 3> 

ln( x^ /x^) 

Notice that since Q(0)P(0) = 3^3 and In 1 =0, 
P(0) 

and, 


] I (9.4) 


m* G U ( s_ ) ' 2 2 2 

] 

ne 


p(oo ) 


m*cW (e^) 0 2 

— - [7^ + S(S+1) 3 


ne 


which are the correct Yoshida limits. 

This is rather fortuitous, because the distribution (8.7) 
and hence (9.5) holds only near T . Rewriting equation (9*4) 

O 

as. 
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P{T) 


,T„2 . ,2 


Ijv + J S(S+1) 


ne 


P - In (3 +t^^Tx^)/(3 +7ik2TT.) 

+ 3j2qp[i - S-i- £— L-] 

In 


- J'^S (S+1) [1 ~ 


p(Tg) = p(oo) - :$ 


In 


(3 +7*:gTT;2)/(3 +7&2 Tt:^) 

J jjJ 

■ In 


where , 

(h;-cj s(S+ 1) [1 in ' (^2^^ 

and hence P(Tg.) '^(oo ). 

O 

Next we calculate the slope of P(T) curve^ Brom 

(9,.4) 


1 d? d_(PQ)r, lnU3 


, 1 r S(S+1) pn-JriL^l^— 

In(T^T^) 3 ” 3 +tc kglTg 


'JiICti'P 


BM 


3+7i;k2T T ^ 


] 


(9.5) 


We can. rewrite equation (9.5) as, 

. lie ' dP fir-on\ r 




■'' ~Ti kgT r ..| 


itkg 


3 


-] 


(3 + 71 k-g^T.^ ) (3 +7tkgTTr2) 


(9.6) 
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Pirst of all, notice that the factors within the brackets f J 
in the equation (9*6) are positive because I^ext 

slightly recasting equation (9.6) we get 

^ _ dfPQ) rinr -^2 s 

3J^clvT_^(6p) ^ 

"^“^3+ iikgT 

. r 5(S+1) . , PQ , r ^ 2 _ ^ '^l i 

3 ^ '-3+ Ttk^^ 'C 2 3 + Ttkgl J 

(9.7) 

Since < 0, and — ^ >PQ ,th8 first term is negative 

whereas the second term is positive. Moreover, - « 1 

as shown in Chapter 7. Remembering that T2 >>r^, it is 

straightforward to show that the factor within square bracket 

[ ] in the first term is much less than that in the second 

term. As a result the second term in (9. 7) (which is positive )» 

dP 

the first term (which is negative), and hence ^ >0. 

dP 

Notice that ^ does not exhibit any maximum at any finite 
temperature, because at all temperatures T < T the second 
tt-rm is much larger than the first and above 1 first term 

S 

d P 

vanishes. Since ;^> 0 even at 1^, there is no anomaly 
in P(T) at T , as obsei-ved experimentally. 

So, we conclude that for T near T_, P(T) is given 

s 

by (9.4) where the approximate form (4.25) for 0'(!P,-^) is 
valid. On the other hand, for 1 »T_ (say for 1^6 1 ) 

(4.24) holds for 'P (1) although (4.25) can no longer be used 
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for 0(1},^ ) as shown in Chapter 4. This is so because for 
f >> 6 only single spins (and no bigger cluster) exist. 
Strongest evidence for cluster foimation below ~ 6 T comes 

from the deviation of X.(T) from Ourie-Weiss behaviour. 

Experimental data on P (T) show a maximum at a 
finite temperature T^> T^ as stated in Chapter 1. This does 
not follow from our model because of two reasons — firstly 
our distribution (8.7) which is valid near T may not be 

s 

applicable near T^, and secondly, our treatment goes only 
upto the order J®. The decrease of P (T) beyond 5L is 
"rendnisoent of the Kondo effect” [3433 for which one has to 
make an analysis upto at least the order and 


beyond . 



CHAPTER 10 


RAiTOOM WALK ARP MAGEBTI2ATIQH DISTRIBUTION 
OF EKKI SFIN CLUSTERS 


By RKKY spin clusters we mean clusters of spins 
interacting via RKKY interaction as defined in Chapter 6, 

In this chapter we shall apply the theory of random walk to 
calculate the probability distribution of the resultant 
magnetization of percolating spin clusters. This method also 
encompasses quantum spins, a generalization never attempted 
be fore . 

The basic problem of an unbiased one dimensional 

random walk can be formulated as follows: "a particle suffers 

displacements along a straight line in the form of a series 

of steps of equal length, each step being taken either in the 

forward or in the backward direction with equal probability 
1 

After taking R such steps, the particle could be at any 
of the points -I, -R+1 , . . , ,R~1 ,l\r. We ask; what is the 
probability W(M,N) that the particle arrives at the point 
a distance M steps away from the starting point after suffering 
R displacements?" [ 382] . The solution of this probleij. is 
well knowm from the consideration of combinatorics; 


* Submitted for publication. 
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w(M,Er) = , (10.1) 

which, in the limit of large N end M« F reduces to 

W(M,F) = exp(-MV2F) (10.2) 

Rivier and Taylor [ 383] 

remarked that 'spin glasses are realisations of a random walk 
which is the simplest of random processes. They, therefore, 
serve as a model for the study of random magnetic matter 
in general' . To our knowledge the idea of random walk has 
not been fully utilized for the study of spin passes. This 
chapter is an attempt, albeit approximate, in thet direction. 


10.1 I sing Spins 


Let us first compare a cluster of I sing spins in 
the absence of an external magnetic field with a one 
dimensional unbiased random walk. 


Random Walk 

i) Total number of steps=F 

ii) Bach step can be forward 
(+) or backward (-) with 
equal probability ^ 

iii) The net displacement 

can have any of the 
values -N,--F+1 , . . . ,F- 1 , 

I. 

iv) We calculate the 
probability W (M,F) that 
the particle arrives at 
the point M after 
suffering N displace- 
ments. 


Ising Spin Cluster 

i) Total number of spins = F 

ii) Bach spin can be either 
up ( + ) or down (-) 

iii) The net moment (M^) can 

have any of the values -F,..., 
F-1,F. 

iv) We calculate the probability 
P(M,F) that the total moment 
is M after counting F spins 
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However, there is a crucial difference between a random 
walk and an Ising spin cluster - two successive steps on a 
random walk are totally uncorrelated, each with probability 
^ for forward and backward move. But in an Ising cluster, 
though the probability of 'up* and 'down' orientations are 
each the successive spins counted are correlated via 
= tanh( p H^). In other woids, as if the drunkard was 
not totally drunk. In order to take this correlation properly 
into account one has to work with conditional probabilities. 

Pq. (£ ) H ) = the probability that a spin at distence R 
from the origin has orientation a when the resulting field 
at the origin is directed along H, as was originally suggested 
by Adkins and Eivier [127] (.AR) in a related context . But 
exact calculation of P^j (R , H) involves higher order 
correlation function and the hierarchy of such correlation 
functions has to be truncated somewhere judiciously for some 
practical calculation. Following IR we shall truncate the 
hierarchy at the very first step by decoupling the joint 
probability using the approximation 

(P+ - P„) (R, H) = sign H . m(T) (10.5) 


AR calculated P(H) for a macro scopically large system 
whereas in our problem the spectrum of cluster size ranges 
from single spin (H=1) to very large values. 
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i . e . equi-valently 

0+ - = M(T) (10.3) 

where and n_ are the number of up and down spins in a 
cluster possessing magnetization M. 

Prom the above consideration it is clear that under 
the approjciiration (10.3), one can write, for large values of 

P(M,]5r) = exp (-mVpN) (10.4) 

This result (10.4) was derived by Abrikosov [372] without 
invoking the idea of random walk explicitly. Note that in 
this case the probability distribution for both the total 
spin as well as the spin projection (magnetization) is govern 
ned by (10.4) since the spin dimensionality is one. 

10.2 Glassical Vector Spins 

The analogy drawn in the preceeding section can easily 
be generalized for the case of classical vector spins which 
may orient in any direction. Prom the theory of random walk 
in three dimensions [382] the probability W(r,N) that the 
position of the particle after N displacements is r is 
given by 

W(r,I) = ( ^ exp[-3rV2N] (10.5) 


Hence we must also have 
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’X / O “ 

P(M,lf) = ( 2 -^) exp [-3 mV2N] (10.6) 

Because of isotropy in direction this leads to 
P(M,];T) d^M = 47i:lA(M,lf) dM 

= 6xp(-3MV2lO<iM (10.7) 

which is again identical with Abrikosov's corresponding 
result [ 3723 . Notice that (10.7) gives us the probability 
of the total spin rather than the spin projection (magneti- 
zation), The probability of the magnetization may also 
be obtained dizectly by integrating over ¥ 1 ^ and iviy to 
obtain 

7 1 - 31i£ 

Although the contents of the sections 10.1 and 10,2 are 
restatements of earlier known results we stress that those 
results can be obtained very easily once the analogy between 
the random walk and spin clusters is recognized. 

The derivations of (10.2) and (10.6) were very 
simple because of the applicability of vector addition for 
coupling of vector spins. This was a direct consequence of 
the commutation of the various spin components. If it is 
possible to visualize the random walk as a process where the 
particle takes small steps 6 x in small time interval gt such 
that x=l 6 x and t=N 6 t then in the continuum limit ( 62 : -» 0 , 

5 t - 0 ) one gets for W(x,t,XQ) a Pokker-Pianck equation. 
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Similarly, a continuum limit is possible for classical spins 
and P(M,l'ir) satisfies the S’ olcke3>- Planck equation (10,2) and 
(10.8) are solutions of such Pokker- Planck equations. 

In the next section we shall see how these ideas 
are to be modified for quantum vector spins. 

10.3 Quantum Spins 

In this section we shall denote the magnetization 
(and spin) values of individual spins by small letters and 
those of clusters by capital letters. 

The first remarkable diffeirence between the classical 
and quantum spins is that the latter can have only (2S+1) 
orientations, where s(s+1) is the eigenvalue of s . This 
leads to two consequences; the ordinary vector addition cannot 
be applied for coupling of quantum spins and secondly the 
continuum limit cannot be taken. In other words, P(M,U) is 
now a discrete, rather than continuous, function of M, 

We shall first consider integral spins. Since we 

are interested in the probability distribution of S = S s. 

i=1 

and these projections obey simple addition law, the problem 
again reduces to that of effectively one dimensional random 
walks. However, there is a difference between this case 
and the one considered in section 10,1. Por Ising spins 
magnetization can change in each step only by +1 or -1 with 
equal probability -g-, and hence only nearest neighbour 
hopping was allowed in the corresponding random walk. In 



214 


the case of Heisenberg spin s, the magnetization can change 
each step by any of the values -s, -s+l , . . . , s-1 ,s. Therefore, 
we have to consider a generalized random walk where the 
'drunkard' can hop a 'distance* -s, -s+1 , . . , , s-1 , s in a 
■§iQ)Sile step with equal probability l/(2s+1). We now write 
for the probability that a cluster of H spins s has 
a net moment M, it satisfies a difference equation 

Pg(M,H) = P^p^ (m-s,N-1) + p^Pg(M-s+1,H-1) 

+ (I^I“SH-2,]!ir-1) + ... + P2gPg(M+s,H-1) 

which can be rewritten in a more compact notation as 

2s 

Pg(M,H) = ^^2 (M-s+l,F-1) (-,0.9) 

where p^ is the probability of a jump of length 1 in one step, 

Por an unbiased random walk p^ = l/(2s-i-1) for all 1, Phe 
ordered set (M,H) forms a two dimensional grid and our aim 
is to determine 1^(1, H) on this grid. The obvious boundary 
condition forall cases is Pg(M,0) = . (10.10) 

Before embarking to obtain a general solution of (10.9) let 
us illustrate with a few examples. 

Example 1 : 

Bet s = 1, H = 2 so that p = 1/3 and P<M,0) = 6^^ 

Besides, from the consideration of normalization we get the 

S P„(M,H) = 1 . (10.11) 

M= -S ® 


sum rule 
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Now, 

P^(M, 1 ) = “ [P^(M-1,0) +P^(M,0) +P^(M+1,0)] 

- 1 r 

" 3 ^M- 1,0 %,0 ^^M+ 1 , 0 ^ 

" 3 ^ "■ ^vi ,0 + %,-1 3 

Hence, P^ ( 1 , 1 ) = P^ (o,1) = P^ (- 1 , 1 ) = and the sum 
rule (10,11) is satisfied for S = 1. 

Next, 

P^(0,2) = 1 [P^ (- 1 ,- 1 ) + P^(0,1) + P^(1,1) ] 

= 5 [5 + 3 + 3 3 _ = 5 

Pl(l, 2 ) = I [P^ ( 0 , 1 ) + p^( 1 , 1 ) + p^ ( 2 , 1 )] 

= 5 [j 5 + 0 ] = 2/9 

Pl(-1,2) = 1 [P^ (-2,1) + P^ (-1,1) + (0,1)] 

= ^[0+^+^] = 2/9 

P-|(2,2) = j [P^( 1 , 1 ) 4. P^(2,1) + P^(3,1)] 

= ^[j + 0 + 0 ] = 1/9 

P^(-2,2) = ^ [ P^ (-3,1) + P^(-2,1) + P^(- 1 , 1 ) ] 

=j [0 + 0 + ^] = 1 / 9 . 

Again the sum rule is satisfied for S = 2. Notice that 
(i) the sum rule (10.11) is satisfied for any N, Thereby 
preserving normalizations, (ii), Pg(M,N) is symmetric with 
respect to the sign of M. 
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Example 2 ; 

let s =2, IT=2. In this case p = l/5 and the boundary 
condition is p2(M,0 ) = q. 

Proceeding exactly as in example 1 we get 


P2(0,2) = 1/5, =P2(-1,2) =4/25, 

P2(2,2) = P2(-2,2) = 3/25, 

P2(3,2) = p2(-3,2) = 2/25 
P2(4,2) = p2(-4,2) = 1/25 
and the sum rule (10,11) is again satisfied. 

Example 3 s 

Let s = 3/2, N=2. Hence p = l/4 and ^^^^2 ^>0)= ^ 

Therefore , 

P3/2(M,1) = 1/4 [P3/2 (M-3/2,0) + - -i, 0) 

+ P3/2(^'^ + ^* ,0) + 23/2^ + 3/2, 0) ] 

4 t\l,3/2 M,i ^ M,- i '^M,--3/2 ^ 


23/2(0,2) = I [P 3 / 2 (- d + P 3 / 2 (- is 1) + 

^ 3 / 2 ^ 1 ’ 5 

-Irl + 1 + 1 + 1^= 1 
“414 + 4^4 4J 4 

^3/2(1? 2) = [^3/2 2* ”l ) + E3/2('2d ) + ^3/2^2’^^ 

= |[| + 4 '^t + o5" ^ " V2 
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^ 3 / 2 ^ 2 >'*^ ^ •^ 3 / 2 ^ 2 '^^ 

■'* ^3/2^2’"’^ ^ 
~ ^3/2 2) 

P3/2(3,2) =j [F^y'2(V2,1) + P3/2^2’”’^ ^3/2^2’^^ 

+ ^ 3 / 2 ^' 2 ^ ^ 3 

^ 4-0 + 0 - f 0 ] = 1/16 = ^ 2^2 


P (M,IT) for any other half integral values of s can be calcu-- 
s 

lated sinrilarly.. 

To get a feeling for the variation of T^(M,lO with 
M as well as N for given s we show the values of P^(M,U) 
in a UjM grid in table 5 ? and figuire 26 . Since the probabi- 

lity is drawn out from the centra to the wings with increasing 
Ef symmetrically, the peak value P^(0,N) decreases with 
increasing N to preserve the sum rule. 


TABLE 5 


■-.N 

M'- 

1 

2 

3 

4 

-4 




1/81 

-3 



1/27 

4/81 

-2 


1/9 

1/9 

10/81 

-1 

1/3 

2/9 

2/9 

16/81 

0 

1/3 

1/3 

7/27 

19/81 

1 

1/3 

2/9 

2/9 

16/81 

2 


1/9 

1/9 

10/81 

3 



1/27 

4/81 

/I 




1/81 
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10.4 G-eneral Solution 

In order to solve the general difference equation, let 
us proceed in two steps. Slrst let us partition the N steps 
into a configuration where there are n^j^ steps of length m 
(m = -9,...,s), The probability of one such configuration 
is given by the multinomial expression 

^3 N) = --g (p)^ (10.12) 

m= -s 

To obtain P (M,W) we can take any configuration which satis- 
s 

fies the conditions 

s 

S = N (10.13a) 

m = -s 

s 

E m n^ = M ( 10.13b) 

m = -s 

Thus, 

P {H,H) = E P ({II ), M) (10.14) 

^ %e(10.13) ® ^ 

where the summations to be carried out over all configurations 
satisfying (10.13). It is straight forward to show that the 
general solution (10,14) satisfies the difference equation 
<10.9). 

Thus we conclude that the spin (orientational) 
configuration in a percolating cluster in a RKKJ SG can be 
thought of as the realization of a random walk in some 
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hyperspace. In addition, the spatial (positional) configur- 
ation of the magnetic atoms in a SG can also be thought of 
as a realization of a random walk. Comparing a EKKY spin 
cluster with a random ferromagnetic cluster at 1=0 we see 
the following difference; the latter represents an ordered 
walk (for spin degrees of freedom) on a randum walk (for 
random position of the magnetic atoms in space), whereas the 
former is the realization of a random walk on a random walk 
[ 384] . Eivier and Taylor [383] described a SG as a 
"drunkard's walk transfixed"; we shall call a SG a "drunkard’s 
walk transfixed on a drunkard's walk". 



GHAPIBR 11* 

t HOLS' IS Hm LQGAI. glEIiD DIS TRIBUITIOIT 

Y/e assumed in Chapter 6 that the number of sites 
where local magnetic field vanishes is vanishingly small at 
sufficiently low temperatures. In other words, we assumed 
the existence of a 'hole* in the local field distribution 
P(H). We shall prove it in this chapter. 

Theorem ; 

In any spin glass there exists a field 6 such that 
P(H) = 0 for 6 at sufficiently low temperature. 

Proof ; 

To prove this theorem one has to go beyond tbe simple 
minded weiss mean field approximation and take the reaction 
field into account. We stress that one does not necessarily 
have to work with the TAP model (or, more generally, with 
any cooperative phase transition model) in order to take 
the effect of tho reaction field into account. The physical 
origin of the reaction field is much more general and can 
be derived as follows; the local field produced at the jth 


* Submitted for publication. 
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site loy the ith spin is field induces a local momeni 

J- ^ j'th site, where X. is the local suscep- 

-LJ J X J 

tibility at the jth site. This induced moment, again, 

2 

produces a local field X at the ith site. Hence the 

J J X 

reaction field produced at the ith site by all other spins 
3(3 ^ 1 ) is S 

If the cavity field at the ith site (i.e. the local field 
at the ith site in the absence of the ith spin) is hj_, then 


h. 

1 


S J. .M. 
13 3 




or 


% 


EJi;^!. 


= h. + 


EJJ. X .M, 


( 11 . 1 ) 


Notice that so far we have not assumed anything about the 
nature of the spin glass transition, i.e., whether it is a 
cooperative phase transition or a gradual freezing process, 
let us define the spatial average of a local property by 

i = (1/N) I 

We may thus replace J^ - X . by AX ('l1.2) 

X J J 

where A is a temperature- independent constant and X ia 
average susceptibility. The right hand side of equation 
( 11 . 2 ) is also independent of i because of statistical 
homogeneity. 

Any homogeneous MHT of GPTM of EA type yields 


Xcc p (1 - Q) 


(11.3) 
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whereas in the percolation model X «: p (1-PQ). Prom Young 
and Kirkpatrick's [I25j computer simulation and Pischer 
and Hertz's [389] discussion (so far as the notation is 

c 

concerned, their q and q correspond to our Q and PQ 
respectiYely) , we get 

1 » Q = (11.4) 

and 

1 - PQ = y/ P (11 *5) 

at sufficiently low teirperaturesi a and y being two constants 
independent of temperature, (11.4) and (11.5) are consistent 
with each other because P <1 for all T> 0, 

We also know that 

= tanh (|3hj_) (11.6) 

because local magnetization is determined by the local cavity 
field hj^ rather than the total Weiss field Therefore as 

T •* 0 


•+ sign (hj^) (11,7) 

Substituting (11.5) and (1-1.7) into (11.1) we get 

% ~ \ Y ^ 

lHj_l = lh^i+ I 6 J . 


or, 

whe re , 


6 


Y 


A 


( 11 . 

( 11 . 8 ) 
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The second term on the right hand side of (11,8) is a 
constant characteristic of the system and does not depend 
on field. Therefore, the magnitude of the local field i.e, 

I H^l can never be zero; in other words, the probability that 
is zero is zero. This proves the theorem. 

In any cooperative phase transition mode 5 since 
|3(1-Q) =« /p as T -*0 and = h^ and there can be no hole 
in P(H). The existence of 'hole’ in P(H) has been observed 
in computer experiments [ 315 , 331 ] • 



CHAPIBE 12 


WHEBE DO W STAND? 


12.1 What have we Achieved ? 

12.1.1 Comparison with other Models 

Although our percolation model appears similar to 
that of Smith [353] there are some crucial differences: 

(a) Smith did not take frustration into account. We have 
shown how T vs. c curve deviates from linearity because 
of frustration, 

(b) Smith model is applicable only to Be the lattice which 
contains no frustration. Our model is applicable to 
any realistic three dimensional lattice as well as to 
amorphous SG, 

(c) Smith suggested, intuitively, that 

X(T) = f (1 - P) 

which differs from our expression (6.16). 

In other words, Smith assumes all spins in the infinite 
cluster to be perfectly stiff (i.e, Q = 1 for all spins 
in the infinite cluster) whereas we allow for thermal 
fluctuations and hence incomplete stiffness. Our expression 
(6.16) has been derived from more fundamental expression 
(6.11) of LET. 
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Our percolation model is qualitatively closer to 
Cyrot's model [389] because the latter takes frustration 
into account. However, there are some differences in the 
quantitative calculations. 

(a) Cyrot, for simplicity, divides the system into "p blocks 
of n spins where np = N (H number of spins)". On the 
other hand, in our percolation model,, there is a 
spectrum of cluster size given by P(n),. 

(b) In Oyrot' s model a SG is distinguished from random 
ferro magnets as well as from lilattis model by the fact 
that there is a distribution PC?), the probability that 
a cluster has free energy S'. This is so because the 
number of frustrated bonds varies from cluster to cluster 
even for same n, and hence the difference in their free 
energy. In our percolation model, we leave out frus- 
trated sites at the very beginning when we construct 
spheres of correlation around the spins. Hence, in 

our model P(n) is non-random at all temperatures. 

(c) Our derivation of Vogel- Fulcher law is quite different 
from that of Cyrot [389]. 

Abrikosov’ s work [372] is complimentary to ours 
as stated in Chapters 6 and 10. The relation between 
Wohlfarth's model and our percolation model has been discussed 
in Chapter 7. So far as neutron scattering experiments are 
concerned, our v/ork in Chapter 6 is a quantification of 
Murani's qualitative ideas [85] • 
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As stated in Chapter 7, our equation (11.8) resembles 
Sompolinsky' s result (l.27d) though the underlying micros- 
copic models are quite different. Sompolinsky' s result is 
valid for only SK model which assumes infinitely weak infini- 
tely long ranged interaction. Wo derived (11.8) for realistic 
three dimensional SC with damped EKKY interaction. All the 
relaxation times diverge in the thermodynamic limit in 
Sompolinsky ' s treatment whereas there is a DRT even in the 
thermodynamic limit in our percolation model. Exactly simila 2 >* 
ly, our expression (6.16) can he recast into the form (1,27c) 
v;here 

Q(x) = 0 for X < X 

= for X > X, with x('r') = i - p(t) 

But this apparent similarity with Parisi’s result (1.27c) 
does not necessarily imply any deeper connoctfon between, our 
percolation model and. SK model, , 

There seems to be a very close relation between our 
percolation model and Anderson and coworkers’ work [320,321] 
that describes SG transition as a localization-delocalization 
problem, Anderson [ 8] comments that "it is these localized 
eigenstates which so convincingly mimic the ’clusters’ of 
the Reel interprelation of SG phenomenology", Aad, we have 
already shovm in Chapter 7 that "it is the percolating spin 
clusters that mimic the ‘grains fins’ of the Wohlfarth model 
so convincingly" . Hence we identify the localized eigen- 
states of the random matrix ^ with the finite clusters and 
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the piohility edge with the percolation threshold. Such a 
relation was also suspected by Young and Stinchoombe [ 390] . 

They also conjectured that the reason why SK model exhibits 
a phase transition is associated with the corresponding 
spectrum of which, in this case, does not contain any 
localized state, which, we think, is true. A rigorous equivalence 
between our percolation model and localization-delocalization 
model is under our investigation. This rigorous treatment 
will make small quantitative changes in our model though 
the qualitative features are likely to remain unaltered. 

Yor example, in equation (6.7 ) Tq will turn out to be temper- 
ature as well as size dependent and the constant factor 
ln(2S+1) will get modified. 

One serious objection against Smith model was that 
clusters become ill-defined for infinite range EKEY inters 
action. This difficulty can be minimized assuming damped 
EKKY interaction, as we have done in Chapter 6, because the 
latter interaction is effectively finite ranged. However, 
definition of clusters in our model is still ambiguous to 
some extent. Hor example, let us consider a cluster of three 
spins, connected by nearest neighbour ferromagnetic bonds J, 
forming a triangle. Such a cluster should freeze when 
the largest solution of the eigenvalue equation - J^- 21^=0 

is equal to k^T, i.e. at a temperature T^ = -7^ , But 

defining 00-) = k J X- 2J^ we find that 0(co)> 0 and 
0(J) < 0,.,so that 0 (^ ) must vanish for a such that 

J< < 00 . In other the cluster under consideration 
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can freeze at a temperature whereas our criterion 

J 

in Chapter 6 was Of course, for low concentration 

B 

alloys, existence of such clusters as considered ahoYe is veiy 
unlikely . 

The localized state picture makes the definition of 
clusters exact, and is identical with our definition for 
highly localized states. Bor weakly localized states near 
T (i.e. near the mobility edge) our clusters becoire fuzzy 

O 

as large clusters overlap. Hence the description of SG 
transition as localization-delocalization problem becomes 
very appropriate. In this picture participation ratio 
corresponds to the size of the clusters. So far as the 
dynamics is concerned, calculation of the relaxation time of 
clusters involves calculation of first passage time. 

12.1,2 Comparison with Experiments 

Our results derived on the basis of percolation model 
agree well with experimental results; 

(a) Our expression (6,15) and (6.16) exhibit a cusp in 
X(T) as observed, 

(b) Our expression (6,8) explains the frequency dependence 

of T^ (w) 

& 

(c) Equation (6.10) explains the q-dependence of 

(d) Bigure 22 agrees with observed T^ versus c data, 

(e ) Our model accounts for the smooth behaviour of specific 

heat across T_, 

a 
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(f) Figure 25 agrees well with the distribution of 
blocking temperatures calculated in Eef.[288], 

(g) logarithmic relaxation of magnetization (8.6) was also 
derived, 

(h) Equation (9*4) accounts for the smooth behaviour of 
resistivity P (T) 

(i) Equation (5.26) explains the M dependence of magneto- 
resistance Ap satisfactorily, 

(j) The value of shown in Table 3 is of the correct 
order of magnitude, 

(k) Equation (11.8) proves the existence of hole in P(H) 
observed in computer experiments, 

(l) In Chapter 6 we also showed that a finite system will 
become paramagnetic on very long time scales, in agree- 
ment with computer simulation results [24,361] 

12.2 What Else Gan be Done ? ; 

As already stated in Section 12.1.1 a rigorous relation 
between localization-delocalization problem and percolation 
model is under investigation. This will enable us to clarify 
some drawbacks of the percolation picture and make more 
rigorous quantitative predictipns wh4.ch can be compared with 
experimental data. 

We have investigated only two transport properties — 
resistivity and magnetoresistance. Other properties e.g. 
thermo power, Hall effect etc. can also be calculated similarly. 
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ITone of the percolation models discussed so far examined the 
effect of spin waves on the stability of the clusters* 

Recently SG- has been considered as a proto-type for 
many more complicated and less understood problems, e.g* 
according to Kirkpatrick, "a computer chip is a SG" and 
according to Hopfield, "the brain is a SG" [391] . Kielsen 
says, "the universe in a SG" and according to Anderson, "the 
origin of life is a SG" [ 391 ]* These comparisons show the 
potentiality of spin glass physics. We have shown in this 
thesis that "a spin glass is a glass" [391] • Howe-ver, a lot 
remains to be done. May be, this 'normal science' will end 
vath a "scientific revolution" (in Kuhnian sense )i 
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APPENDIX 1 


Q(T=0’ 


aiANTDM SPINS 


It was shown by Fischer [150] that 
Q = S(S+1) y - 7 ' A-- / d^r exp(- 4 r^)rt ^ 

X 

[ (S + ■^) Goth (S + ^ Coth('| 

wheDTC Z, = (pA )^Q 

and A is the second moment of the gaussian distribution 
of random bonds, As T - 0, p •* oo and hence 5 •* oo 

Besides, { (S + ■^) Goth (S + r- Goth ^ r) J **■ S 

as T -* 0, 

Thus, Q = S(S+1) at T = 0. 
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APEEKDIX 2 


Q(T=0) FOR CLASSICAL VEOTOR SOT S 


It was shown by Mookerjee [392] that in the spin glass 


phase 


00 Yi 

Q = Q = Q - ( ( C dh^ dh^ dh^ P(h) [Ooth(p h)- 

^yy ^zz J .ho6 xyz^n ^ pn 

At T = 0, (T=0) = Qyy (T=0) = (^=0) = Q(T=0) 




(2J-2q)3/2 


■z^2 g-hV2GJ^Q 


2^ 


(PixcJ-j Q; 


(X) tC 0,0 

C JU f* u2 . ^ n 2^ -h /2 cJ 5Q 
jyg j dh I dO h sxn6 Cos © e 1 

^ 0 h) 


2 

3 * 


00 
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^ ii2 dh 


Putting 


h‘ 


2cJ^ Q 


r= X, 


Q_(T=0) = Q,„,(1=0) = Q,,(T=0) = 


2^2 


•XX"- “yy" '' 'zz'“ ' ^(2^)172 


00 


dx 


X I . p(3/2) = i . for 


= 1 


In the general case, 


XX 


^^yy "■ ^ZZ“ 5 


at T = 0 



